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All-Band 2G+3G Radial Disc-Cone Antennas: 
Design, Construction and Measurements 

N.I. Yannopoulou, P.E. Zimourtopoulos, E.T. Sarris * 

Antennas Research Group, Austria - Hellas [l, 2] 

EECE Dept, Democritus University of Thrace, Hellas [2 and 3] 

Abstract 

We define as "All-Band 2G+3G" any band that includes all 
frequencies allocated to both 2G and 3G services. We define 
as "Radial Disc-Cone Antenna RDCA" any discone antenna with a 
structure of radial wires. The RDCA was theoretically ana- 
lyzed and software simulated with the purpose of computation- 
ally design a broadband model of it. As an application, a 
broadband RDCA for operation from 800 to 3,000 MHz, which in- 
clude all 2G and 3G frequencies, was designed and an experi- 
mental model was constructed and tested. In order to evaluate 
the agreement between theory and practice, mathematically ex- 
pressed measurement error bounds were computed. 

Introduction Research 

In 1945, Kandoian invented The RDCA was theoretically 

the well-known discone an- analyzed as a group of iden- 
tenna, that is a dipole made tical filamentary V-dipoles 
of a disc above a cone [1] . with unequal arms connected 
In 1953, Nail gave experimen- in parallel. The dipoles re- 
tally two naive relations for cline on equiangular vertical 
the discone dimensions [2]. phi-planes around z-axis to 

In 1987, Rappaport de- form a disconical array, 
signed discones using an Fig.lA shows two such copla- 
N-type connector feed [3]. In nar dipoles conformed with 
1993, Cooke studied a discone the apex angle a. Each 
with a structure of radial V-dipole has a total length L 
wires [4]. In 2005, Kim et equal to the sum of arm 
al. presented a double radial lengths r and s plus the gap 
discone antenna for Ultra 9 between its terminals. 
Wide-Band applications [5] . The simulation was based 

In this short paper we on a suite of developed vis- 
present an All-Band 2G+3G RDCA ual tools which are supported 
fed by an N-type/Female/50- by a fully analyzed, cor- 
Ohm connector. rected and redeveloped edi- 
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tion of the original thin- 
wire computer program by J.H. 
Richmond [6] . 

Two arithmetic criteria 
were adopted for the broad- 
band characterization of a 
model : 

(1) 50-Ohm VSWR lower than 2 

(2) Normalized radiation 
intensity U/Umax lower 
than 3 dB on the 
horizontal plane. 

A visual application pro- 
gram was specifically devel- 
oped to design a broadband 
radial discone with bare 
wires of diameter d embedded 
in free space when the wire 
conductivity, the type of 
feeding connector and the 
frequency band are given. 

The program uses the model 
of a radial discone fed by an 
N-type connector shown in 
Fig. 2. Starting with an ap- 
propriate combination of the 
relations given by [2] -[4] 
the program computes by iter- 
ation in terms of wavelength 
A, the geometric characteris- 
tics r, s, g, a, of the 
broadband model, just when 
the criteria are satisfied. 

Fig. IB shows a Ground 
Plane Antenna GPA that was 
designed for reference and 
consists of equal number of 
cone radials s and a vertical 
monopole with height r. 

As a practical application 
of the broadband design, the 
2G+3G band from 800 to 2,500 



MHz was selected to begin 
with and an experimental ra- 
dial discone of copper wire 
fed by N-type connector was 
built, as shown in Fig. 3. 




Fig . 1 : A - RDCA, B - GPA 



z 




Fig. 2: RDCA - Designed Model 
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ALL -BAND 2G+3G RADIAL DISC-CONEANTENNAS 



In order to demonstrate 
the particular behavior of 
the experimental model, the 
2G+3G band was divided as 
follows : 



2G+3G Sub-Bands 
80Q MHz - 2,500 MHz 



Sub-Band 


MHz 


I 


806 - 960 


II 


1,429 - 1,513 


III 


1,710 - 1,900 


IV 


1,910 - 2,025 


V 


2,110 - 2,170 


VI 


2,400 - 2,499 



Our measurement system 
consists of an EM anechoic 
chamber, a network analyzer, 
a number of support instru- 
ments, a set of standard 
loads of factory accuracy and 
a constructed antenna rota- 
tion mechanism with a built 
hardware control unit of its 
step motor. The combined 
characteristics of system 
parts specify a measurement 
band from 600 to 1300 MHz, 
which overlaps with the 2G/3G 
band. Developed control soft- 
ware synchronizes the system 
and collects data using the 
IEEE-488 protocol. 

A developed general mathe- 
matical method expresses the 
measurement error bounds. An- 
other set of developed soft- 
ware applications processes 
the collected data and com- 
putes the error bounds. 



Results 

The consideration of ra- 
dial discone as an array of 
at least eight 8 V-dipoles 
produces a theta-polarized 
vector radiation pattern with 
magnitude a surface almost by 
revolution around z-axis. So 
the radial discone has indeed 
on the horizontal plane xOy 
the basic properties of a 
vertically polarized almost 
omni-directional antenna, 
that is a fact that encour- 
aged the design of a broad- 
band model by using simula- 
tion . 

The application of the 
broadband criteria to 2G/3G 
band resulted to the design 
of a RDCA with the following 
geometrical characteristics: 



All-Band 2G+3G RDCA 
800 MHz - 3,000 MHz 



Geometry 


Units 


d 


1.5 


[mm] 


r 


44 


[mm] 


g 


6 


[mm] 


s 


125 


[mm] 


a 


60 


[°] 



The RDCA operates from 800 
to 3,000 MHz, which exceeds 
that of 2G+3G band. The ac- 
cordingly constructed experi- 
mental radial discone of 
Fig. 3 should be implied with 
a constructional tolerance of 
±0.5 mm and ±0.5°. 
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The broadband model has a 
directivity from about -0.5 
to 2.9 dBd with slope angle 
between -65° and +58°, but 
the directivity gain on hori- 
zontal plane stays very close 
to the desirable value of 0 
dBi, since it changes from -1 
to +1.7 dBi only. Fig. 4 shows 
that the predicted horizontal 
normalized radiation inten- 
sity remains below 3 dB in- 
deed, while it stays above 0 
dB relative to the reference 
antenna in all 2G+3G sub- 
bands indicated by the verti- 
cal gray strips, when both 
are fed by the same 50-Ohm 
Fig. 3: RDCA experimental model source. 

L/A 

co 0.35 0.70 1.05 1.40 1.75^ 

T3 CD 




Frequency, GHz 



Fig. 4: Predicted radiation intensity on horizontal plane. 
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Fig.5A shows the predicted 
normalized radiation patterns 
in dB at the center of each 
sub-band, which confirms the 
horizontal omni-directional 
radiation properties of the 
broadband model. 

At the center frequency of 
950 MHz of the measurement 
band, the predicted and mea- 
sured radiation intensity on 
the three main cuts of the 
radiation pattern are in good 



agreement, as shown in Fig. 
5B . 

This is made clearer by 
the measurement error bounds 
on a vertical plane as shown 
in Fig. 6. 

Fig. 7 shows that the 50- 
Ohm VSWR predicted results 
for the broadband discone are 
below 2 indeed and almost 
covered by the error bounds 
in the measurement band. 




Fig. 6: Measurement error bounds on a vertical plane 
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Fig. 5: A (Up) Predicted normalized radiation intensity patterns 
at the center of each 2G+3G sub-band - B (Down) Normalized 
radiation intensity pattern at the center of measurements band 
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L/A 




0.6 0.8 1 1.3 2 2.5 3 3.2 

Frequency, GHz 



Fig. 7: Standing wave ratio against frequency or ratio of total 

length to wavelength 



Conclusion 

Prediction and experimen- 
tation in the measurement 
band 600 MHz to 1,300 MHz 
proposes a successfully de- 



signed, constructed, and mea- 
sured Radial Disk Cone An- 
tenna RDCA capable to serve 
All-Band 2G+3G applications 
from 800 MHz to 3,000 MHz. 
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Abstract 



An analytical method was developed to estimate errors in 
quantities depended on full one-port vector network analyzer 
(VNA) measurements using differentials and a complex differ- 
ential error region (DER) was defined. To evaluate the 
method, differences instead of differentials were placed over 
a DER which was then analyzed and compared with another com- 
monly used estimated error. Two real differential error in- 
tervals (DEIs) were defined by the greatest lower and least 
upper bounds of DER projections. To demonstrate the method, a 
typical device under test (DUT) was built and tested against 
frequency. Practically, a DER and its DEIs are solely based 
on manufacturer's data for standard loads and their uncer- 
tainties, measured values and their inaccuracies. 



Introduction 

In full one-port measure- 
ments with a VNA of real charac- 
teristic impedance Z G , a DUT 
with impedance Z has a reflec- 
tion coefficient p defined by 

p= (Z-Z 0 )/(Z + Z 0 ) 

and related to its measured 
value m by the bilinear trans- 
formation 

p = (m-D)/[M(m-D) + R] 

in terms of errors D, M and R 
[1] . This transformation can be 
uniquely determined from given 
distinct p n , n = 1, 2, 3 and re- 
spectively known m k , k = n [2] . 



Research 

We considered p n , m k as the 
elements of given ordered 
triples (A, B, C), (a, b, c), 

solved the resulting system 
and appropriately expressed 
its solution by 

F = X cC(B - A) 

D = X abC(A - B)/F 

M = X c(B - A)/F 

R= [l"l (A - B) (a - b) ]/F 2 

where X and |"l produce two 
more terms from the one 
shown, by rotation of the or- 
dered triple elements. These 
errors were then considered 
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as depended on the indepen- 
dent variables p n , m k . There- 
fore, their differentials 
were expressed in the same 
manner by 

dD = [n (a - b) X (B - C)BC dA 
+ X (b - c) 2 (B - A)(C - A) 

BC da] /F 2 

dM = [X (a - b)(c - a)(B - C) 2 dA 

- n (A - B) X (b - c)da]/F 2 
dR = {X [F + 2(a - b)B(A - C)] 

[(B - C) 2 dA n (a - b) 

- (b - c) 2 da n (A - B) ] }/F 3 

After that, the differential 
of p was expressed by 

dp = [-RdD - (m - D) 2 dM - (m - D)dR 
+ Rdm]/[M(m - D) + R] 2 

and was considered depended, 
through dD, dM and dR, on L = 
7 independent variables and 
their independent differen- 
tials: p n , n = 1, 2, 3 and m k , 
k = n or k = 0 with m 0 = m. 

The developed expressions 
were mechanically verified 
using a developed software 
program for symbolic computa- 
tions . 

Manufacturer's data for 
standard loads used in full- 
one port VNA measurements are 
substituted in p n , and for 
their uncertainties in dp n . 
Since Z 0 is real, the domain 
of each p n is the closed unit 
circle [3]. For |p n | = 0 or 1, 
care must be exercised to re- 
strict its differential value 
onto its domain. The VNA mea- 



surements have specified 
bounded ranges for their mo- 
dulus and argument, so that 
the domain of each m k is a 
bounded circular annular with 
its center at the origin O of 
the complex plane. Measure- 
ment data are substituted in 
m k and manufacturer's data 
for measurement inaccuracy in 
dm k . Uncertainty and inaccu- 
racy data outline domains for 

dp n and dm k . If z = | r | e- 1 ^, stands 
for any of the independent 
variables and dz for its dif- 
ferential then the contribu- 
tion of dz to dp is a summa- 
tion term of the form Wdz, 

with W = | U | e^ v , so that 
Wdz = |U|e j(v + <p) d|r| 

+ |u|e j(v + (p+ TT/2) | r | dtp 

where W is in fact a known 
value of the respective par- 
tial derivative and d|r|, dtp 
are the independent real dif- 
ferentials of the complex dz 
in polar form. Each expres- 
sion Wdz outlines a contour 
for a partial DER around O. 
If z ^ 0, the partial DER is a 
parallelogram with perpendic- 
ular sides d | r | and | r | dcp, 
stretched or contracted by |U| 
and rotated by (V + tp) around 
O. If z = p n = 0, the partial 
DER is a circle with radius 
| U | d | r | . Accordingly, a DER 
is the sum of either L paral- 
lelograms or ( L - 1 ) parallel- 
ograms and 1 circle. DER is 
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then a convex set with con- 
tour either a polygonal line 
with 4L vertices at most, or 
a piecewise curve composed of 
4( L - 1) line segments and 4(L 
- 1) circular arcs at most. 
The greatest lower and least 
upper differential error 
bounds are the end-points of 
DEIs for the real and imagi- 
nary parts of dp and result 
from the projections of DER 
for p on the coordinate axes. 
These conclusions can be ge- 
neralized for any other quan- 
tity directly or indirectly 
depended on all, some or just 
one of the above independent 
variables and their differen- 
tials. Thus, the quantity has 
an L-term DER, where 7>L>1. 
For example, the impedance Z 
of a DUT has the 7-term DER: 

dZ = 2Z 0 dp/(l - p) 2 

Results 

All of the following data 
are specified by manufactur- 
ers of the parts for our mea- 
surement system. This system 
operates from 1 to 1300 MHz 
with 100 Hz PLL stability and 
consists of a type-N Z 0 = 50 Q 
network analyzer, a number of 
support instruments and a set 
of standard loads. The stan- 
dards are: a short circuit A, 
a matching load B with re- 
flection coefficient 0.029 
and an open circuit C with 
reflection coefficient 0.99 
and phase accuracy ±2°. In 
the absence of manufacturer's 



data for A we considered its 
uncertainty equal to that of 
C. So, the following values 
were substituted in the de- 
veloped expressions: 

A = -1, 0 < d | A | <0.01, -180° < 
dtp A < -178° or 178° < dtp A < 180°, 

B = 0, | dB | = 0.029, 

C = 1, -0.01 < d | C | <0, -2° < dtpc 
< + 2 ° 

The annular domain for m k 
of VNA is specified from 0 to 
-70 db in modulus and ±180 
degrees in argument. Measure- 
ments m k result with a deci- 
mal floating point precision 
of 4 digits, for both modulus 
and argument. We consider the 
modulus and argument of dm k 
equal to ±1/2 of the unit in 
the last place of the corre- 
sponding mantissa in modulus 
and argument of m k Conse- 
quently, our system produces 
a DER, either for p or Z, as 
a sum of ( L - 1 ) = 6 parallelo- 
grams and 1 circle, with a 
contour of (4L + 4L) = 48 ver- 
tices at most. 

A suite of developed soft- 
ware applications: (i) con- 

trols the system and collects 
the data in terms of fre- 
quency using the IEEE-488 
protocol, (ii) processes the 
collected data and computes 
the vertices of DER and the 
end-points of its DEIs (iii) 
sketches pictures of DER for 
p and Z in terms of the fre- 
quency steps and make a film 
using them as frames. 
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77.8 86.7 95.5 

Input Resistance R, Ohm 



Fig.l: A typical DER for the impedance Z 




0.2 

Fig. 2: DER for the reflection coefficient p and for its 
associated impedance Z against frequency 
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Fig. 3: Greatest lower and least upper differential error 
bounds for resistance R and reactance X against frequency 



A typical resistor with a 
nominal DC impedance of 50 Q 
±20% tolerance was soldered 
on a type-N base connector 
and enclosed in an aluminium 
box to serve as a simple DUT 
for testing its Z from 2 to 
1289 MHz in 13 MHz steps. The 
center frequency f c = 639 MHz 
was chosen to reveal the de- 
tails of the proposed method 
in Fig.l, where the contour 
of a typical DER for Z is 
outlined with small circles 
as its vertices. This contour 
surrounds that of the 4-terms 
DER due to inaccuracy of mea- 
surements (1) and that of 3- 



terms DER for the uncertainty 
of loads (2). A properly cir- 
cumscribed rectangle of DER 
shows graphically how the 
DEIs for R and X result. The 
commonly used error from the 
matching load only is shown 
as a dotted circle. This is 
in fact a 1-term DER which is 
surrounded from the contour 
of the DER by a factor of 
about 125% to 185% in all di- 
rections. Finally, in the 
same figure, 2 7x2 differences 
AZ resulting from the same 
dp n and dm k , dense enough to 
appear as stripes, were 
placed over DER to compare 
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them with differential dZ 
values. Notably, almost all 
of AZ values are belong to 
DER while the computation 
time for these AZ exceeds 
that for DER by more than one 
order of magnitude. To demon- 
strate the method, a set of 
selected DER frames for p and 
Z are shown in Fig. 2, as 
beads on space curved fila- 
ments against frequency. 



Finally, the computed DEIs 
for R and X are shown in 
Fig. 3 against frequency. 

Conclusion 

The proposed method may be 
efficiently used in the same 
way, to successfully estimate 
errors in any quantity de- 
pended on full one-port vec- 
tor network analyzer measure- 
ments . 
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Abstract 



Since S-parameter measurements without uncertainty cannot 
claim any credibility, the uncertainties in full two-port 
Vector Network Analyzer (VNA) measurements were estimated us- 
ing total complex differentials (Total Differential Errors). 
To express precisely a comparison relation between complex 
differential errors, their differential error regions (DERs) 
were used. To demonstrate the method in the most accurate 
case of a direct zero-length thru, practical results are pre- 
sented for commonly used Z-parameters of a simple, two-port, 
DC resistive T-network, which was built and tested against 
frequency with a VNA measurement system extended by two 
lengthy transmission lines. 



Introduction 

It is well known that in 
full two-port VNA measure- 
ments the S-parameters for a 
two-port Device Under Test 
(DUT) are given in terms of 
their 4 measurements m^, i = 1, 

2, j = 1, 2 by 

~~ £[(^11 - D ) / R] [1 + (m 22 
-D' )M'/R' ] - L(m 21 - X)(m 12 
- X' )/(TT' )}/H (1) 

S 21 = {[l+(m 22 -.D')(M' - L ) /R ' ] 
(m 21 - X)/T}/H (2) 

H = [1 + (m^ - D)M/R] [1 + (m 22 

- D 1 )M ' /R ' ] - LL' (m 21 - X)(m 12 

- X ' ) / ( TT ' ) (3) 



S 22 , S 12 have expressions 

that result from (l)-(2) by 
substituting i, j with j, i 
and D, M, R, L, T, X with D' , 
M', R 1 , T', L', X' and vice- 
versa [1] . These 12 quanti- 
ties have been defined as 
system errors [2] . Stumper 
gave non-generalized expres- 
sions for the partial devia- 
tions of S-parameters due to 
calibration standard uncer- 
tainties, in 2003 [3] . Fur- 

thermore, the developed total 
differential errors for full 
one-port VNA measurements [4] 
are also not generalized in 
the two-port case. To the 
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best of the authors' knowl- 
edge, there are no analytical 
expressions for total differ- 
ential errors in full two- 
port VNA measurements. 



Research 

Since S-parameters are 
functions of 16 complex vari- 
ables, their total differen- 
tial errors were initially 
expressed as 



ds ii = {T T ' (i - MS 11 ) [ R ' + M' (m 22 - D' )] (dm 1:L - dD) 

- RR 1 L(1 - L 1 S 1;L ) [ (m 21 - X) (dm 12 - dX' ) + (m 12 - X' ) (dm 21 - dX) ] 

+ M ' T T 1 [ ( m lrL - D ) ( 1 - MS 1:l ) - RS 1:L ] (dm 22 - dD' ) 

- TT'S 1 i(m 11 - D)[R' + M ' (m 22 - D ' ) ] dM 

+ TT'(m 22 -D')[(m 11 -D)(l - MS 1:L ) - RS 1:L ]dM' 

- ( R ' L( 1 - L'S 1:l ) (m 12 - X' ) (m 21 - X) + 

+ TT'S 11 [R' + M ' (m 22 - D ' ) ] )dR 

- (RL(1 - L'S 11 )(m 12 - X' )(m 21 - X) 

-TT'[(m n -D)(l - MS lrL ) - RS lrL ] )dR' 

- RR' (m 12 - X' ) (m 21 - X) [(1 - L'S 1;L )dL - LS^dL' ] 

+ [(m lrL - D)(1-MS 11 ) - RS 1± ][R' +M'(m 22 - D')] 

(T ' dT+TdT ' ) }/P (4) 

dS 21 = {- MT T'S 21 [R' + M' (m 22 - D' )] (dm 1:L - dD) 

+ RR' LL'S 21 (m 21 - X) (dm 12 - dX' ) 

+ R{T ' [R' + (m 22 -D')(M'-L)] + R ' LL ' S 21 (m 12 - X' )}(dm 21 - dX) 

+ T' (R(m 21 - X) (M' - L) - M'TS 21 [R + M(m 1;L - D)] ) (dm 22 - dD' ) 
-TT'S 21 (m 11 -D)[R' + M ' (m 22 - D ' ) ] dM 
+ T' (m 22 - D' ) (R(m 21 - X) - TS 21 [R + M(m 1;L - D) ] )dM ' 

+ { (m 21 - X) (T ' (m 22 - D ' ) (M ' - L) + R ' [T ' + LL 1 S 21 (m 12 - X ' ) ] ) 
-TT'S 21 [R' + M' (m 22 - D' )]}dR 
+ (R(m 21 - X) [T ' + LL'S 21 (m 12 - X' )] 

- T T'S 21 [R + M(m 1;L - D)] )dR' 

+ R(m 21 - X) [R ' L'S 21 (m 12 - X') - T' (m 22 - D' ) ]dL 
+ RR' LS 21 (m 12 - X' ) (m 21 - X)dL' 

- T'S 21 [R + M(m 1;L - D)] [R* + M' (m 22 - D' )]dT 
+ (R(m 21 - X) [R 1 + (m 22 - D ' ) (M ' - L)] 

- TS 21 [R + M(m 11 - D) ] [R ' + M'(m 22 - D' )])dT'}/P (5) 
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P = TT'[R , +M , (m 22 -D , )][R+H(m 11 -D)] 

- RR'LL' (m 12 - X' )(m 21 - X) (6) 



dS 22 and dS 12 resulted from 

(4), (5) with the mentioned 

substitutions. X, X' errors 
stand for crosstalk measure- 
ments. D, M, R (D 1 , M', R 1 ) er- 
rors are uniquely determined 
in terms of 3 standard loads 
A, B, C (A 1 , B 1 , C') and their 
3 measurements a, b, c (a', b', 
c'), by full one-port VNA mea- 
surements, so the number of 
independent complex variables 
increases from 16 to 22. L, T 
(L', T') errors are accurately 
determined after the replace- 
ment of DUT with a direct 
thru (or approximately, if an 
adapter is used instead) in 
terms of new measurements 
tn, ^21 ( t 22 ' ^- 12 ) anc ^ °f previ- 
ously found quantities. Their 
expressions were appropri- 
ately stated as 



L= [I(ab + ct 11 )C(B-A)]/E (7) 

T= (t 21 -X)[n (A - B) (a - b) ]/ 

(E [IcC(B- A)]) (8) 

E = X(ab + Ct 11 )(B-A) (9) 

where X and |"| produce two more 
terms, from the given one, by 
cyclic rotation of the let- 
ters a, b, c (a 1 , b', c') or A, 
B, C (A ' , B ' , C ' ) . In this way, 
each S-parameter has as total 
differential error dS, a sum 
of 22 differential terms: 

16 due to measurement inaccu- 
racies dm^, dX, dX', d t , da, 

db, dc, da', db', dc' and 6 due 
to standard uncertainties 
given by their manufacturer 
dA, dB, dC, dA', dB 1 , dC'. The 
expressions for dD, dM, dR 
(dD 1 , dM', dR') are known [4]. 
The expressions for the rest 
of differential errors were 
developed as 



dL = {X (B - C)(b - t 1:L )(c - t 1;L ) [ (B - C)(b - a)(c - a)dA 

- (b - c)(B - A)(C - A)da] + [n (A - B)(a - b)] dt^j/E 2 (10) 

dT = {X (t 21 -X)(b- c)(B-C)[(t 11 - c)(b- a) 2 B(A 2 +C 2 ) 

+ (b - t lrL )(c - a) 2 C(A 2 + B 2 ) - 2ABC( b - c)(t 1:L (b + c - 2a) 

- be + a 2 )] [(B - C)(b - a)(c - a)dA 

- (b - c)(B - A)(C - A)da]}/(E 2 [X cC(B - A)] 2 ) 

+ [n (A - B)(a- b)]{[(t 21 -X) Xa(B- C)/E]dt 11 

+ dt 21 - dX}/( E [X cC( B - A) ] ) (11) 
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Each complex differential 
error defines a Differential 
Error Region (DER) on the 
complex plane with projec- 
tions to coordinate axes the 
Differential Error Intervals 
(DEIs) [4]. Obviously, any 



quantity differentiably de- 
pendent on the above vari- 
ables has also a DER. For ex- 
ample, after another correc- 
tion to the given S to Z-pa- 
rameters relations [5], the 
Z-DERs are resulted from 



dZu _ 2Z 0 [ ( 1 S 22 ) 2 dS 11 + (1 S 22 )S 21 dS 12 + (1 S 22 )S 12 dS 21 

+ S 12 S 21 dS 22 ] / [(l - S 1:L )(1 - S 22 ) - S 12 S 21 ] ( 12 ) 

dZ 21 - 2Zg[(l - S 22 ) S 21 dS 11 + S 21 dS 12 + (1 - S 11 )(l - S 22 )dS 21 

+ (l-S 11 )S 21 dS 22 ]/[(l-S 11 )(l-S 22 ) -S 12 S 21 ] 2 (13) 



while dZ 22 , dZ 12 result from 

(12), (13) by application of 

the mentioned substitutions. 

Results 

Six calibration standards, 
in pairs of opposite sex, 
were used and their manufac- 
turers' data were substituted 
in the developed expressions: 

A = -1 = A 1 , 0 < d | A | = d | A 1 | < 
0.01, -180° < dcp A = dcp A , < -178° 

or 178° < dcp A = dcp A , < 180°, 

B = 0 = B' , | dB | = 0.029 = | dB ' | , 
C = 1 = C ' , -0 . 01 < d | C | = d | C ' | 
<0 and -2° < dcp c = dcp c , < +2° . 

The inaccuracy of any VNA 
measurement was conserva- 
tively considered as a sym- 
metric interval defined by 
just 1 unit in the last place 
of the corresponding man- 
tissa, both in modulus and 
argument. Consequently, each 



S-DER is a sum of 20 paral- 
lelograms and 2 circles, with 
a contour of 160 vertices at 
most [4] . 

To demonstrate the method, 
a typical T-network of common 
resistors with nominal DC 
values Z 1 = 24.2 Q, Z 2 = 120 Q 

for the horizontal arms and 
Z 12 =l. 1 £2 for the vertical arm, 

were soldered on type-N base 
connectors of opposite sex 
and enclosed in an aluminium 
box, to form a two-port DUT. 

The VNA measurement system 
was extended by two transmis- 
sion lines of 3.66m and 14m, 
respectively, up to the DUT. 
The DUT was tested from 2 to 
1289 MHz in 13 MHz steps. The 
frequency 1003 MHz was se- 
lected to illustrate the pro- 
posed method for S-DERs shown 
in Fig. 1. 
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To study the total differ- 
ential error, dS was expressed 
as dll + dl, where dll is due to 
the uncertainty of 6 stan- 
dards and dl to the inaccu- 
racy of 16 measurements. The 
contribution of these, con- 
servatively considered mea- 
surement inaccuracies to the 
total differential error is 
as much significant as the 
uncertainties of standard 
loads are. For example, com- 
putations for S 12 over the 

whole measurement band show 
that max|dU| and max|dl| con- 
tribute ~35%-80% and ~25%-70% 
to max|dS 12 |, respectively. 

In addition. Fig. 1 shows how 
the projections of each S-DER 
result its real and imaginary 
DEI. To display the variation 
of S-DER against frequency, a 
number of selected S-DER 
frames are shown in Fig. 2 as 
beads on a space-curved fila- 
ment. It is worth mentioning 
that S 1:l -DER (S 22 -DER) was 

greater than it resulted from 
appropriately organized full 
one-port measurements, as it 
was expected. Finally, the 
computed Z-DEIs are shown in 
Fig. 3, along with their LF 
Z-values . 






Fig . 1 : S-DERs at 1003 MHz 
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1 400 800 1300 1 400 800 1300 

Fig.3:Z-DEIs against frequency 



Conclusion uncertainties in any case 

The proposed method may be where the process equations 
efficiently used to estimate (1), ( 2 ) and (4), (5) can find 

application. 
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Direct and Related Far-Field Inverse 
Scattering Problems for Spherical 
Electromagnetic Waves in Chiral Media 

N.M. Berketis * 



Independent Researcher, Athens, Greece 

Abstract 



This paper studies the direct and inverse scattering 
problem when the incident electromagnetic field is a time 
harmonic point - generated wave in a chiral media and the 
scatterer is a perfectly conducting sphere. The exact 
Green's function and the electric far-field patterns of the 
scattering problem are constructed. For a small sphere, a 
closed-form approximation of the scattered wave field at the 
source of the incident spherical wave is obtained. Also 
treats the same inverse problem using far-field results via 
the leading order term in the low-frequency asymptotic ex- 



pansion of the scattering cros 

Introduction 

In a homogeneous isotropic 
chiral media the electromag- 
netic fields are composed of 
left - circularly polarized 
( LCP ) and right - circularly 
polarized (RCP) components, 
which have different wave 
numbers and independent di- 
rections of propagation. 

The LCP and RCP components 
are assumed to be spherical 
Beltrami fields since in 
practice such wave fields are 
more readily realized. 

In this work, the author 
has studied the electromag- 
netic waves in chiral media 
produced by a point source in 
the vicinity of the scat- 



-section . 

terer. In particular, in [1], 
[2], reciprocity, optical and 
general scattering theorems 
for stimulation of point- 
source asymmetric media have 
been demonstrated. This paper 
studies the inverse problem 
of far field [2] . Specifi- 
cally, we measure the scat- 
tering cross-section for a 
five-point source area. 

In the second Section, 
considering Bohren decomposi- 
tion into suitable Beltrami 
fields, we formulate the di- 
rect scattering problem of a 
spherical electromagnetic wa- 
ve by a perfectly conducting 
obstacle. This problem is 
well posed, the existence and 
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uniqueness has been proved in 
[1], [3]. 

In the third Section, af- 
ter expanding the incident 
field in terms of spherical 
wave functions, we obtain the 
exact solution of the scat- 
tering problem as well as an 
expansion for the electric 
far-field pattern [2]. 

Finally in the fourth Sec- 
tion, we consider either LCP 
or RCP incidence and we ob- 
tain an approximation of the 
scattering cross-section. For 
the far-field experiments, we 
measure the scattering cross- 
section for various point 
source locations [2] . 



vector r 0 with respect to an 
origin 0 in the vicinity of 
the scatterer. In order to 
define spherical electromag- 
netic fields E r , H r , we make 
use of the Bohren decomposi- 
tion into Beltrami fields 
Q/_,r 0 and Q RI . , as follows 

E r 0 = Q(.,r e + Qr, r 6 

1 ( 1 ) 
H r e ~ ^ (Qt, r 0 " QR,r e ) 



Where 


q = (p/e) 1/2 is 


the 


in- 


trinsic 


impedance 


of 


the 


chiral 


medium. The 


Beltrami 



fields satisfy the equations 
[4], [5], 



Statement of the problem 

Our goal is to study the 
direct and inverse scattering 
problems when the incident 
electromagnetic field is a 
time harmonic point - gener- 
ated wave in a chiral medium 
and the scatterer is a per- 
fectly conducting sphere of 
radius a centered at the 
origin. The exterior space 
(r = |r| > a) is an infinite ho- 
mogeneous isotropic chiral 
medium with chirality measure 
(3, electric permittivity e 
and magnetic permeability p. 

We consider a time har- 
monic spherical electromag- 
netic wave due to a point 
source at P g with position 



Vx Ql-,r 0 - Y lQl^q 

x q r, r 0 = — y rQ-l,^ 



where y L andy R are wave num- 
bers given by, 



Yl = 



k 

1 - /cp ' ” 



k 

1 + /cp 



(3) 



with k = cj(ep) 1/2 , w being the 
angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form 
[ 1 ], [ 2 ], 
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DIRECT AND RELATED FAR- FIELD INVERSE SCATTERING PROBLEMS 



1 1 



QZ£(r \Pl) = -(i + — vv + — Vxi) ( 

2 Yl Yl 



hjYi |r ~ r B |) 
MVV 0 ) 



) ' P L = 



2 nr g 1 ^ Lr ° ~ ~ 

^ [/ t G fs (r, r 0 ) + VxG fs (r, r 0 ) ] • p t 



1 ~ 1 1 

-(i + 4vv- — 

2 Yl V R 



A(Yr |r-r B |) 

MW 0 ) 



2tt r p _1 ' v ° 

^ [y R G fs (r, r 0 ) - VxG fs (r, r 0 ) ] • p R 



( 4 ) 



(5) 



where b 0 (x) = b 0 (x) = e lx /( ix) is 
the zeroth-order spherical 
Hankel function of first 

kind, I = xx + ijjijj + zz is the 
identity dyadic, r 0 =|r o | and 

G fs (r, r 0 ) is the free space 
dyadic Green function [4] . 
The constant unit vectors 
p L and p R satisfy the rela- 
tions 

r B P L = r e p R = 0 

r 0 xp L = ip L (6) 

r 0 x p R = - ip R 

We note that when r 0 -»oo, the 
incident electric field 

Ef; c (r|p t ,p H ) = 

(7) 

= Q4 ( r I Pl) + Q“ro ( r I P«) 



reduces to plane electric 
wave with direction of propa- 
gation -r a and polarizations 
p L , p R , since 

limQ4(r|pJ = e 1 ^ r p,= 

r 0 ^» ( 8 ) 

= Qf c (r;-f 0 ,pJ 

limQ^(r|p R ) = e- ly ^ r p R = 
r ^“ (9) 

= Qf c (r; -r 0 , p R ) 

We consider E^" c is incident 

r o 

upon a perfectly conducting 
sphere of radius a . Then, we 
want to calculate the scat- 
tered electric field E* c , 

which is the unique solution 
of the following exterior 
boundary value problem 
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( 10 ) 
( 11 ) 

The Silver-Muller radiation condition is modified as follows 
rxVxE= c (r)-P/ 2 rxE= c (r) + ^E= c (r) = o(i), r -> ® (12) 



uniformly in all directions 
r g S 2 , where S 2 is the unit 

sphere in M 3 , n is the out- 
ward normal unit vector on 

the scatterer and Y 2 = Y l Yr- 
The scattered electric field 
will be depended on the po- 
larizations p L , p R and will 
have the decomposition 

Er 0 c (r | P„ p R ) = Q£ t (r | p„ p R ) 

, . (13) 

+ Q"r 0 (r I P R ) 

where Qi, c r e O' I P^ p r ) and 



The functions g Lr and g R 
are the LCP and RCP far-field 
patterns respectively [4], 
[ 6 ], 

If either a LCP or a RCP 
spherical electric wave 
Er" c ( r I Pa) , A = L,R, is inci- 
dent upon the scatterer, then 
the scattering cross-section, 
is given by [6], 



A,r„ 



f [4 |g 4 r e (f I Pj | 2 + 

S 2 Y L 

+ \ |g R ,r 0 (r I P jf] ds(r) (16) 



Q£. ( r I Pl, Pr) are the corre- 
sponding scattered Beltrami 
fields which have the follow- 
ing behavior, when r oo 

Qi, c r 0 ( r \Pl'Pr) = h g (y L r)- 

. . . ( 14 ) 

gi,r 0 (r I Pi, P R ) + 0 (p) 

Q^r e (r|Pi, P R ) = h 9 (y R r)- 

... -i ( 15 ) 

•g R ,r 0 (r|Pi,P R ) + O(^) 



Exact Green's function 



We take spherical coordi- 
nates ( r,6,cp ) where 6 e [0, n] 
and e [0, 2n), with the ori- 
gin at the center of the 
spherical scatterer, so that 
the point source is at 

r = r e , 0 = 0. 



Thus, r 0 = r 0 z , p L = 



12 ( * - ^ 



and 


Pr J2 


(x 


+ lip) , 


where 


x, ip 


and 


z are 


unit 


vectors 


in 
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the x, ip and z directions, 
respectively. Using spherical 
vector wave functions, [4], 
[5], [7], and taking into ac- 
count (4), and (5), we obtain 
[ 2 ], 

oo 

( r I Pj - X B n {^clniyj) + 

(17) 

+ iL ( e! n (y,r)} 



or 



Qi?"r e ( r I Pr) = Z B n {RLiUKr 1 ") - 



n= 1 (18) 

-IRLiUKrO} 



where, for r < r g 

1 2n + 1 



2V2 h B {y A r B ) n(n + 1) 



W„(VV e ) 



and 



(19) 



H n i.Y A r B ) = h n (y A r B ) - i h n (y A r B ) 

with A = L, R . The h n is a 
spherical Hankel function of 
first kind of order n, 

h{x) = x~%{x) + h' n (x ) , 

and R^ n , with s = e, or 
o (even or odd) are the 
spherical functions [4], [5], 

L si„(W r ) = M si„(W r ) + N si n (W r ) # N 

( 20 ) 

RsiU^) = M^UKrT) - Ni? n (y R r) 

where p = 1, 3 , the and 

N^ n are known spherical vec- 
tor function [7]. The scat- 
tered electric field that 
comes from a LCP incident 
field or a RCP incident field 
has a similar expansion to 
(17) or to (18) 



E r 0 C ( r I Pj = Qm- 0 (r I Pi.) + Qr^ s (r I p L ) = 

00 00 
= Z B n a n i^olniYj) + j-O/jO } + £ {R® (KrO + 1R®, (K r 0 } 



( 21 ) 



or 



(T I Pr) = Q^r a (r I Pr) + Q“r. ( r I Pr) + 

00 00 

= Z BX{L^ n (y L r) - 1 L< 3 > n (y L r) } + Z « {R®, (y R r) - iR^for)} 

n= 1 n = 1 



( 22 ) 



Using the boundary condition (11) on r = a, we obtain [2], 



MONDAY 18 NOVEMBER 2013 



V2-11 



FUNKTECHNIKPLUS # JOURNAL 




N.M. BERKETIS 



a L Jn (Y L a) K (Ml + J'n (/l 3 ) h n (V R S ) 

h n (y L a) h n ( y R a ) + h n ( y,a ) h n {y R a) 

and 

a * J n (Y L a) K (V L a) ~ I (/^) b n (y L a) 
h n (y L a) h n (y R a) + h n ( y,a ) h n (y R a) 

or 

b L Jn (ygg) /i„ (/ R a) - Jn fog) K foa) 

K foa) h n foa) + h n foa) h n (y R a) 

and 

b R Jn (y R a) K foa) + j„ foa) /?„ fog) 
h n (y L a) h n (y R a) + h n (y L a) h n (y R a) 

Using the asymptotic forms [4], 



(23) 



(24) 



(25) 



(26) 

[7], 



L ( s 1„ for) ~ >(n + l)(-i) n b 0 (y L r) ffo (r ) ( 27 ) 

R ( si n for) ~ Vn(n + l)(-i) n b 0 (y R r) f s R ln (r) ( 28 ) 

where let us introduce LCP Beltrami angular ffo(f) , and RCP 
Beltrami angular f R ln (r) [4], satisfy by relations. 



fsin(r) = C sln (r) + iB sln (r) , f s R ln (r) = C sln (r) - iB sln (r) 
we calculate the electric far-field patterns [2], 



g"r 0 (r I Pl) = £ 



(2n + l)(-i) n 1 H n (y L r B ) 
ti 2^/2 n(n + 1) h B (y L r g ) 



afoffo(f) -iffo(f)} 



(29) 



(30) 



or 



9"r 0 (r I P«) = Z 



(2n + l)(-i) n 1 
2^2n(n + 1) 



H n (VVa) 

h B (y R r e ) 



bU-fLir) 



iffo(f)} 



(31) 
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A far-field inverse problem 

So far, all of our formu- 
las are exact. In the asymp- 
totic results to follow, 

there are three parameters 
y A a , with A = L, R and r = a/r B . 
We note that the geometrical 
parameter r must satisfy 

0 < r < 1 because the point 
source is outside of the 

sphere . 

We assume that |y^a| <k 1, as 
well; that is we make the so- 
called low-frequency assump- 
tion. From (23), (24), (26), 

(27), we obtain [2], 



ai = - 1 (V L a ) 3 + 0{{y L af) 

a 2 = - i 1 g 0 P/c ( y L a ? + 0((y L aY) 
y L a^>0 (34) 



i(l - p/c) 3 
2(1 + p/c) 2 



(y L a) 3 + 0((y L a) 5 ) 



i(l - p/c) 4 
36(1 + p/c) 3 



(y,a) 5 +0((y,a) 7 ) 



y L a -> 0 



(35) 



L _ 1 + / 1 / a \2 

" 2i£ 2 (2n + 1) ^ 



R i (1 - p/c) n + 2 
a! ^ ^^(y.a) 2 



2nC (l + p/c) n 



y L a 0 



(32) 



or 






i (l + p/c) n + 2 
2n<C 2 (1 - p/c) n + 1 



(' V„ a ) 



2/7 + 1 






i(i-P^) 
2^(2n + 1) 



(/«a) 



2n + 1 



or 






i(l + p/c) 3 
2(1 - p/c) 2 



(y R a) 3 +0((y L a) 5 ) 



i(l + p/c) 4 
36(1 - p/c) 3 



(K R a) 5 + 0((y L a) 7 ) 



y R a -> 0 



(36) 



R= _i(l-P/0 3 + )5) 

6 

< 

b " = - l(1 9 0 P/C) ^ a ) 5 +O((w0 7 ) 

y R a^0 (37) 



y R a^0 (33) 

where 

<■„= l-3-5-(2n - 1) = (2n)!/(2”n!) . 

In particular, 



In order to calculate 

with an error of 0((y„a) 4 ) we 
only need the following [4] 
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C 0 u(r) = —r= (cos <p0 - cose sin<p<p) 
v2 



B n11 (r) = — r= (cos 6 sin<p0 + 

V 2 (38) 

+ cos<p<p) 



B ell (r) = —j= (COS0 sin<p0 - sin<p<p) 
V2 



C eii(f) = (~ sincpO - 

- COS0 cos<p<p) 



(39) 



C 012 (r) = (3/2) 1/2 - 
(cos0 cos <p0 - cos2 0 sin<p<p) 

B el2 (r) = (3/2) 1/2 - 
(cos20 cos <p0 - cos0 sin<p<p) 

B 012 (r) = (3/2) 1/2 - 
(cos2 0 sin<p0 + cos0 costpcp) 

C el2 (f) = (3/2) 1/2 - 

•(- cos0sin<p0 - cos20cos<pcp) 



So from (30) and (31), we finally obtain [2], 

g“r e (r I pj = *A — — [- (V^a) T^f^f) + a^if^ (r) ] + 

O 

+ (Y A a) 2 { 2ir[f e /, 11 (r) + s^if^ (r ) ] + r 3 [f^ 12 (f) + m A if 0 A 12 (r) ] } + 

+ (Y A a) 3 { 2[f e /, 11 (r) + ^if^ (r ) ] + r 2 [f,f 12 (f) + a^if^ (r ) ] } ] + 

+ o((/,a) 4 ) 

and 



g"r e (r I P„c) = 



3(1 - tn,(3/c) 2 



a (y /1 a)r 2 [f e A 11 (f) - ra.iOr) ] 



8(1 + OT^p/c) 

(W*) 2 {- ro, 31(1 - ^k) T[f ; ii( p) _ ra ,if^(f) ] - 



5i (1 TH^P/C) 3 p x -x ■ rA /~\ -| -1 

f— T L * el2 (O * ol2 v ) J J 



16V3 1 + mfik 

+ (Y A a) 3 {-™ A 



3 r 3( 1 + TU^p/c) r ■ f A /~\ -| 

. L’ellvO 



(40) 

(41) 

(42) 

(43) 

(44) 



(45) 
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- w 



5 (1 w&k) ^ [f . z(f) _ ^if^f) ] } + 0((y,a) 4 ) 



where m,= 



12^3 
-1, A = L 



1, A = R 



and if A = L, R at that case A c =R r L. 



Now the scattering cross-section, by LCP or RCP spherical 
Beltrami fields, from (16), is given by the relations [2], 



<r. = | l~2 te (r I Pjf + A fe (r I Pjf] ds(f) 



Yl 



Yr 



-(1 + P k ) 2 r 16 



= (na 2 ) C‘ r 4 + (y,a) 2 (M r 2 + 16r B ) + (y,a) 4 (^ + ^ r 4 ) ] 



\ 4/64 , 16 _4\ 



64 



(1 - (3/c) 4 3 



(1 + (3/c) 2 l 4 
+ (Vi a ) 4 (3 + 



[- r + (y L a) (3r + 



5(1 - (3/c) 2 



r 6 ) 



4/0 , 5 ( 1 P^) T 4 nl , „//., „\6 

9(1 + (3/c) 2 



16(1 + (3/c) 2 

r 4 )]} + 0((y L af), y L a^O 



(46) 



or 



= J [A- fe I P*)f + |gj?r 0 (r I P R )f] ds(r) = 



Yl 



Y R 



= (na 2 ){ (1 6 ^ /C) r 4 + (y R a) 2 (^ r 2 + I6r 6 ) + (y R a) 4 (^f + ^ r 4 ) ] + 



\ 4/64 . 16 t-4-i 



(1 + P/C) 4 r 3 



(1 - (3/c) 2 l 4 
+ (W*) 4 ( 3 + 



[- r + (y R a) ( 3r + 



5(1 + (3/c) 2 



r 6 ) + 



4/0 , 5(1 + (3/c) \ m , ,-,//., „\6 

9(1 - (3/c) 2 



16(1 - (3/c) 2 

r 4 )]} + 0((y R a) 6 ), y R a 0 



(47) 



In the special case r G ^- oo (r 0) , by the relations (46), 
(47) we obtain [2], 



cr = (na 2 ) { 



(1 + (3/c) 2 3(1 - (3/c) 4 



(1 + (3/c) 2 



} (Vi a ) 4 + 0((y L a ) 5 ) , y L a ^ 0 



(48) 
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or 

a = (TT a 2 ) { (rL ~ 3 P/C) + 3 ^p ^2 > (V « a ) 4 + °((Y R af) , y R a -> 0 (49) 

likewise in the case y L a^>0 and y R a^Q, by the relations 
(46), (47) we obtain [2], 



<C. = 4 K) ( na2 ) ( a / r 0 ) 4 



where A = L, R , with f„(( 3, /c) = 



(1 - ^p/c) 2 [ 3(1 + ^[3/c) 4 
3 (1 - cr A p/c) 2 



Choose a Cartesian coordi- 
nate system Oxipz, and five 
point-source locations, name- 
ly (0,0,0), (1,0,0), (0,1,0) 

(0,0,1) and (0,0,21), which 
are at (unknown) distances 
r 0 , r 1 , r 2 , r 3 and r 4 , respec- 
tively from the sphere's cen- 
ter. The parameter 1 is a 
chosen fixed length. For each 
location, measure the leading 
order term in the low- 
frequency expansion of the 
scattering cross-section. 

Thus, our five measure- 
ments are [2], 

mj = ±f A V, k) naHf-y (52) 

J = 0, 1, 2, 3, 4 

Dimensionless quantities 
related to m ] are 



Y = -J— = 2 I i 1 (^-) 2 

S f A^k)na ( a> (5 , 
j = 0, 1, 2, 3, 4 

equivalently, we obtain [2], 



r j 2 = | *0 n T 

J = 0, 1, 2, 3, 4 



(54) 



There are six unknowns 
namely r 0 , r 1 , r 2 , r 3 , r 4 and a. 
Furthermore, r 0 , r 3 and r 4 are 
related using the cosine 
rule, r 4 2 + r 0 2 = 2r 3 2 + 21 2 . So, we 
can find the six unknowns. 
The center of the spherical 
scatterer is obtained from 
the intersection of the four 
spheres centers at (0,0,0) 
(1,0,0), (0,1,0) and (0,0,1), 
with corresponding radius r 0 , 
r 1 , r 2 , r 3 respectively. 



FUNKTECHNIKPLUS # JOURNAL 



V2-16 



ISSUE 2 - YEAR 1 




DIRECT AND RELATED FAR- FIELD INVERSE SCATTERING PROBLEMS 



References 

[1] Athanasiadis C., Berketis N . , "Scattering relations for 
point - source excitation in chiral media". Mathematical 
Methods in the Applied Sciences, Vol. 29, 2006, pp. 27-48 

[2] Berketis M. Nikolaos, "Scattering spherical Electro- 
magnetic waves in chiral media", PhD thesis. University 
of Athens - Department of Mathematics, December 2007 

[3] Arnaoudov Y, Dassios G, Had j inicolaou M., "The resistive 
coated sphere in the presence of a point generated wave 
field". Mathematical Methods in the Applied Sciences, 
Vol. 22, 1999, pp. 73-90 

[4] Lakhtakia A., "Beltrami Fields in Chiral Media", World 
Scientific, Singapore, 1994, §2-1.3 (pp. 42-44), §2-2.3 
(pp. 49-52), §4-1 (p. 130 Chapter 4), §5-1 (p. 164 Chap- 
ter 5), §6-1.1 (pp. 214-397), §8-2.1 (pp. 391-397), §9- 
2.3 (pp. 468- 473). 

[5] Lakhtakia A, Varadan V.K., Varadan V.V., "Time-Harmonic 
Electromagnetic Fields in Chiral Media", Lecture Notes 
in Physics, Vol. 335, Springer, Berlin, 1989, §3-§4-§7- 
§11-§12-§13 

[6] Athanasiadis C., "On the far field patterns for electro- 
magnetic scattering by a chiral obstacle in a chiral en- 
vironment", Journal of Mathematical Analysis and Appli- 
cations, Vol. 309, 2005, pp. 517-533 

[7] Morse P.M, Feshbach H., "Methods of Theoretical Physics, 
Part II", Me Graw-Hill, 1953, pp. 1864-1866, 1870, 1871, 
1887 



Preprint Versions 

"Direct and inverse scattering problems for spherical 
Electromagnetic waves in chiral media" 
Christodoulos Athanasiadis and Nikolaos Berketis 
"http : //arxiv . org/abs/0812 . 2169" 



MONDAY 18 NOVEMBER 2013 



V2—1 7 



FUNKTECHNIKPLUS # JOURNAL 




N.M. BERKETIS 



* About the Author 

Nikolaos Berketis, was born in Athens, Greece, 1952. 
Mathematics degree, Aristotle University of Thessaloniki, 
1976. MSc in Applied Mathematics-Numerical Functional Analy- 
sis, Department of Mathematics, National and Kapodistrian Uni- 
versity of Athens (UoA), 2003. Ph.D. in Applied Mathematics, 
Department of Mathematics, UoA, 2007. Professor of Secondary 
Education, 1980-2006. High school Principal, 2007-2012. Tea- 
cher at the UoA: Department of Informatics and Telecommuni- 
cations, 2002-2003. Postgraduate program of the Department of 
Geology and Paleontology, 2003-2004. Research experience: 
Participation in the following research projects: 1) "The 

far-field operator in solving inverse scattering problem". 
Special Research Account 70/4/6412, UoA, 2002. 2) "Scatter- 
ing of spherical electromagnetic waves by chiral materials". 
Special Research Account 70/4/6412, UoA, 2004. 3) "Mathe- 
matical analysis of wave propagation in chiral electromag- 
netic and elastic media." Pythagoras II: Reinforcing of Uni- 
versity research groups ( EPEAEK II), Department of Mathemat- 
ics, UoA, 2005. Member of the Question Bank department of 
Evaluation and Training Research Center, Hellenic Mathemati- 
cal Society ( HMS) , 2001-2004. Editorial Board experience: 

Member of the Scientific Committee of Euclid C' journal, 
Hellenic Mathematical Society (HMS). Member of the Scien- 
tific Committee of the project "Digital Training Textbooks 
for Pan-Hellenic examinations" -for Mathematics, Science and 
Technology Educational Path, 3rd grade, High School and 
Mathematics & Elementary Statistics, General Education, 3rd 
grade, High School- implemented by Computer Technology In- 
stitute University of Patra, for the Ministry of Education. 
Eighteen published peer-reviewed papers, more than 150 ref- 
erences, in peer-reviewed journals, utilizing both Hellenic 
and international bibliography. Participation in the pro- 
ject: "Change of speed of pulse wave after surgical narrow- 
ing of the descending thoracic aorta, "1 st prize, Best Pro- 
ject award, 28th Annual Pan-Hellenic Medical Conference in 
the field of basic research, oral presentation, Athens, May 
21-25, 2002. Personal work on the scattering of spherical 

electromagnetic waves in chiral media is posted at "The 
Smithsonian/NASA Astrophysics Data System" website. Author- 
ing experience: Scientific research articles and books. 

nberketis@gmail . com 

This paper is licensed under a Creative Commons Attribution 4.0 
International License - https : //creativecommons . org/licenses/by/4 . 0/ 

FUNkTECHNIkPLUS # JOURNAL 



V2-18 



ISSUE 2 - YEAR 1 



ELECTRICAL ENGINEERING 



EXPERIMENTATION 



Experimental Results on the Behavior of 
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Democritus University of Thrace, Xanthi, Greece [l, 2] 
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Abstract 



This paper investigates the influence of various parame- 
ters on the behavior of water droplets on polymeric surfaces 
under electric fields. An inclined plane test was carried 
out to understand the droplet behavior in strong electric 
field. Parameters such as, water droplet conductivity, 
droplet volume, polymeric surface roughness and droplet po- 



sitioning with respect to the 
flashover voltage is affected 
ters. The droplet positioning 
than the droplet volume. 

Introduction 

Water droplets on a poly- 
meric surface may cause 
corona under the influence of 
an electric field and can 
cause deterioration to insu- 
lation surface even in condi- 
tions of low pollution level. 
Water droplets on a polymeric 
surface increase locally the 
applied electric field. Local 
field intensifications lead 
to partial discharges (PD) 
and/or localized arcs, which 
may render possible the dry 
bands on the polymeric sur- 



electrodes were studied. The 
by all aforementioned parame- 
is in some cases more vital 



face. Local arcing will even- 
tually bridge the dry bands 
and a complete flashover will 
finally ensue. This mechanism 
is valid to a greater or 
lesser extent for both out- 
door and indoor insulation, 
although each of the afore- 
mentioned categories have 
their own particular charac- 
teristics, namely that indoor 
insulation is stressed more 
and is subjected to a differ- 
ent type of environmental in- 
fluences than outdoor insula- 
tion [1], [2]. A combination 
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of water droplets and dust- 
like impurities on the sur- 
face of a polymeric surface 
may lead to a conducting con- 
tamination layer, which may 
cause a reduction of the 
flashover voltage. The design 
of high voltage insulators, 
they can be for indoor or 
outdoor use, one should take 
into account not only the 
pollution level, the insula- 
tor material and the appro- 
priate voltage level, but 
also the influence of water 
droplets on the flashover 
voltage. Previous work, car- 
ried out in this laboratory, 
tackled the behavior of water 
droplets for a wide range of 
water conductivities (1.7 
pS/cm - 10000 pS/cm) [3], [4]. 
In both publications, it was 
shown that, among the factors 
influencing the behavior of 
the water droplets, were the 
water conductivity, polymer 
surface roughness, droplet 
volume and droplet position- 
ing with respect to the elec- 
trodes . 

In the present work, a 
study of the aforementioned 
parameters on the water 
droplet behavior under the 
influence of a uniform elec- 
tric field in the range of 
1.7 pS/cm - 2000 pS/cm was 

carried out. All tests were 
performed with an inclined 
test arrangement, in order to 
simulate the behavior of wa- 
ter droplets on the surface 



of a real insulator. The an- 
gle used with respect to the 
horizontal was 10°. Such an 
angle was chosen because of 
its immediate relevance to 
industrial insulators. 

Force balance at the droplet/ 
polymer surface interface 

A modeling of a wet con- 
taminated surface was given 
in other publications and 
only a brief outline is pro- 
vided here [5] . Condensation 
of droplets on the surface of 
a high voltage insulator can 
come about from droplet 
germs. In Fig. 1, the forces 
exercised on the droplet are 
shown in case where no elec- 
trical field is applied. Such 
forces are the surface ten- 
sion of the liquid (t l ), the 
surface tension of the solid 
(t s ) and the interfacial ten- 
sion between liquid and solid 
(5 sl ) ■ When an electric field 
is applied, the droplet de- 
forms because of an addi- 
tional force. The tangential 
electric field on the surface 
of the insulator creates a 
force on the surface of the 
droplet which causes its de- 
formation. The deformation of 
the droplet affects the field 
distribution. Local field in- 
tensifications may result, 
which will cause micro-dis- 
charges between the droplets. 
This is the beginning of the 
chemical deterioration of the 
insulator surface. 
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Contact angle Water drop 




Fig. 1: Force balance at the interface solid/liquid at a 
water droplet on an insulating Surface (after [5]) 



Hydrophobicity may locally 
be lost. The voltage differ- 
ence across the droplet will 
be diminished and micro-dis- 
charges will follow. Solvable 
nitrates, which are the re- 
sult of the electrochemical 
deterioration, cause a higher 
conductivity of the water 
droplets. Dry zones may fol- 
low. It is important to bear 
in mind that not only the in- 
fluence of the applied elec- 
tric field on the shape of 
the droplet is of great sig- 
nificance, but also the in- 
fluence of the disintegrated 
droplet on the electric field 
distribution [5], [6]. 

Hydrophobic polymeric sur- 
faces are characterized by a 
low surface conductivity 
which in turn gives a low 
discharge activity and a 
higher flashover voltage. 
This holds also for polluted 
environments. Reduced hy- 



drophobicity implies a higher 
risk for flashover of the in- 
sulator. Hydrophilic materi- 
als, on the other hand, are 
very sensitive to polluted 
environments, and are charac- 
terized by a significant ac- 
tivity of local discharges 
[7] . It is well known the 
classification of the Swedish 
Transmission Research Insti- 
tute ( STRI ) regarding the hy- 
drophobic and hydrophilic 
surfaces. STRI classifies the 
various surfaces according to 
their hydrophobicity from 
class 1 (most hydrophobic, 
with only discrete droplets 
on the surface with contact 
angle larger than 80°) to 
class 7 (most hydrophilic, 
with continuous water films 
forming on the surface). How- 
ever, the truth is that no 
matter if the insulator has 
some sort of humidity and it 
is only slightly polluted or 
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it is heavily polluted, in 
both cases surface discharges 
play a most important role, 
and such discharges may start 
from water droplets. 

Experimental arrangement and 
preparation of the samples 

The aim of this paper is 
to study the behavior of wa- 
ter droplets under the influ- 
ence of an electric field. 
The voltage supplied was from 
a 20 kV transformer (in prac- 
tice the transformer may de- 
liver voltages up to 1.2 
times of its nominal voltage 
without loss of the accuracy 
of the measurement. Conse- 
quently, the applied voltages 
were accurate up to 24 kV) . 
The electrodes used were of 
copper. A top view as well as 



a cross section of an elec- 
trode is shown in Fig. 2. The 
electrodes were half cylin- 
drical in shape. Attention 
was paid to the smoothness of 
the electrode surfaces, so 
that no unnecessary field en- 
hancements could be noticed. 

The water droplets were 
positioned on the polymeric 
material surface with the aid 
of a special arrangement con- 
sisting of a metallic frame 
and three rules, one of which 
had two laser indicators. The 
water droplets were put on 
the surface with a syringe. 
Detailed information on the 
way the droplets were posi- 
tioned on the polymeric sur- 
face is given in [3] . The 
photograph of the inclined 
plane test is shown in Fig. 3. 





Fig. 2: Top view (above) and cross section (bottom) of the 
electrodes used (all dimensions in mm) 
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Fig. 3: The inclined Plane test setup (side view) 



The polymeric materials 
used were PVC, rubber and 
silicone rubber. Surface 
roughness and resistivity of 
the material were measured. 
Surface roughness were mea- 
sured using perthometer (Type 
Perthometer M4P) . They gave a 
roughness of 0.25 pm for PVC, 
0.79 pm for silicone rubber 
and 1.10 pm for rubber. Re- 
sistivity of the material 
were performed with a Megger 
( BM25 type) and they gave a 
resistivity of 206 GQ for 
PVC, a resistivity of 3100 GQ 
silicone rubber and a resis- 
tivity of 2660 GQ for rubber. 
The above values of surface 
roughness and surface resis- 
tivity were not isolated val- 
ues, but each of them was the 
mean of three measurements 
[8], [9], 

In the present work, by 



mixing known quantity of NaCl 
in distilled water forming 
solutions with conductivity 
in the range 1.7 pS/cm to 
2000 pS/cm, were used as 
droplet. The range of conduc- 
tivity were chosen based on 
the conductivity of natural 
rain and its values lie in 
the range 50 - 150 pS/cm, 
whereas the tests with porce- 
lain and glass insulators are 
performed with conductivities 
of 2500 pS/cm [10] . 

Experimental procedure 

The materials used were 
PVC, silicone rubber and rub- 
ber. Various droplet arrange- 
ments were studied. These ar- 
rangements are given in Fig. 
4. Each droplet had a volume 
of 0.2 ml. The electrodes 
were positioned at a distance 
of 4 cm from each other. 
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Fig. 4: Top view showing the droplet arrangements. Starting 
from top left, the arrangements were named as 
(1) arrangement with 1 droplet, 

(2a) arrangement with 2 droplets, 14-12-14, 

(2b) arrangement with 2 droplets, 12-16-12, 

(2c) arrangement with 2 droplets, 8-24-8, 

(3) arrangement with 3 droplets, 

(5) arrangement with 5 droplets and 

(9) arrangement with 9 droplets 

All dimensions given are in mm and they symbolize 
the distances of the droplets from the respective 
electrodes and the distances between them 



The parameters investi- 
gated were the water conduc- 
tivity, the roughness of the 
insulating surface, the posi- 
tioning of the droplets and 
their volume. The insulating 
surfaces were used as they 
were received from the manu- 
facturer without any further 
treatment. After putting the 



droplets on the polymeric sur- 
face, the voltage was slowly 
raised until flashover oc- 
curred. After that and after 
cleaning the surface and put- 
ting new droplets on it, the 
voltage was raised again up 
to the previous flashover 
value minus 1.2 kV, so that 
no new flashover would occur. 
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At this voltage the arrange- 
ment would stay for 1 min. If 
no flashover occurred, the 
voltage was raised by 0.4 kV 
and the procedure was re- 
peated until flashover oc- 
curred. The reason we left 
every time the voltage on for 
1 min, was in order to give 
necessary time interval for 
the droplet(s) to deform and 
for the partial discharge to 
initiate . 

It should be noted that it 
is observed a tendency for 
the droplets to slide, espe- 
cially for PVC because of its 
smooth surface. The droplet 
slide was minimal in the case 
of rubber, which was the 
rougher of the three materi- 
als used. An elongation of 
the droplets was observed, as 
the applied voltage was 
larger. A more evident oscil- 
lation of the droplet was ob- 
served with silicone rubber. 
The reason for that was be- 
cause the aforementioned ma- 
terial is more hydrophobic 
than the other two. Conse- 
quently, the droplet, for a 
defined droplet volume, has a 
smaller contact area with 
silicone rubber, and for this 
reason it oscillates more 
[8] . In some cases, such as 
with PVC with a droplet con- 
ductivity of 1.7 pS/cm and 
with the arrangement (1) of 
Fig. 4, ejection of minute 
charged droplets was observed 
just before flashover [11]. 



Experimental results 

At first, experiments were 
performed without any 
droplets between the elec- 
trodes. This was done in or- 
der to have reference values 
of the flashover voltage and 
also to understand influence 
of number of droplets between 
the electrodes that would re- 
sult in a reduction of the 
flashover voltage. The flash- 
over voltages without any 
droplets measured were 23 kV 
(± 0.5) for PVC, 25 kV (± 
0.5) for silicone rubber and 
24 kV (± 0.5) for rubber. The 
flash-over voltages of the 
three materials used were 
very similar. 

In Figs. 5 - 11 the varia- 
tion of flashover voltage 
with respect to the droplet 
conductivity for different 
droplet arrangements is 
shown . 

It is evident that sili- 
cone rubber presents a higher 
flashover voltage than the 
other two materials. It 
should be noted, however, 
that in the case of droplet 
arrangements (5) and (9) 
where rubber seems to be as 
good as silicone rubber. A 
possible explanation might be 
that in such a case, the 
droplets cover a significant 
part of the polymeric surface 
and hence they play an even 
more important role than the 
polymer itself. This in com- 
bination with the fact that 
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the rubber has a rougher sur- suit the lesser oscillation 
face compared to the other in the case of rubber. 




0 500 1000 1500 2000 2500 

Conductivity (|jS/cm) 

Fig. 5: Flashover voltage for droplet arrangement (1) 




Fig. 6: Flashover voltage for droplet arrangement (2a) 
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Flashover voltage for droplet arrangement (2b) 




Fig. 8: Flashover voltage for droplet arrangement (2c) 
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Flashover voltage for droplet arrangement (3) 




Fig. 10: Flashover voltage for droplet arrangement (5) 
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Fig. 11: Flashover voltage for droplet arrangement (9) 



The better performance of 
silicone rubber is due to its 
hydrophobicity [5], [12]. The 

larger contact angle the 
droplets have minimum contact 
with the insulation material 
as in the case of silicone 
rubber. Figs. 12 - 14 show 

the influence of droplet vo- 
lume on the flashover volt- 
age. It is clear that the 
number of droplets affects 
the flashover voltage, i.e. 
the larger the number, the 
lesser the flashover voltage. 
An exception to that we have 
with the arrangement of 3 and 
5 droplets. It is evident 
that larger flashover volt- 
ages were observed with 5 
droplets than with 3 drop- 



lets. A possible explanation 
of that is that in the case 
of 3 droplets, the distance 
between electrode and droplet 
is smaller than in the case 
of 5 droplets. Consequently, 
one might say that there are 
occasions where the position- 
ing of the droplets with re- 
spect to the electrodes which 
plays a more vital role than 
the whole droplet volume. A 
further validation of the 
above consists of the compar- 
ison of the flashover volt- 
ages in the cases of 3 and 9 
droplets. It is observed that 
the flashover voltages for 
both these arrangements are 
not that different although 
the droplet volume triples. 
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An interesting case con- 
sists also of the droplets 
arrangements 2a, 2b and 2c. 
Higher values for flashover 
voltage were observed for ar- 
rangement 2a, then for ar- 
rangement 2b and the lower 
flashover voltage was ob- 
served for droplet arrange- 
ment 2c. This fact reinforces 
the above observations, name- 
ly that the positioning of 
the droplets play a crucial 
role, i.e. the closer the 
droplets to the electrodes, 
the lower the flashover volt- 
age. It is to be noted that 
similar observations were 
made also in [3], [4], where 

not an inclined arrangement 
was used but a horizontal 
one. What is presented in 
this paper is an approach of 
the behavior of water drop- 
lets on polymeric surfaces 
with an inclined electrode 
arrangement. The results were 
reproducible but not that 
many tests were carried out 
which would allow a statisti- 
cal study of the collected 
data. The main interest of 
this paper concentrates on 
the study of the behavior of 
the droplets. In the present 
context, no emphasis was gi- 
ven to the quantification of 
the studied parameters. 



A comment should be made 
on the results with the con- 
ductivity of 1.7 pS/cm: in 
such a case, with such a low 
conductivity, the water path 
behaves like a load, i.e. 
like a resistance connecting 
the two electrodes. It is for 
this reason that we observed 
a lowering of the voltage at 
the output of the power 
source. The flow of current 
through water of low conduc- 
tivity (i.e. of large resis- 
tance) means practically an 
increase of water temperature 
because of the power loss in 
the resistance of the water 
path. From the relation P = 
I 2 R, we can conclude that as 
the resistance is larger, as 
in our case the water path of 
very low conductivity, the 
power loss at this resistance 
is larger. Consequently, the 
temperature developed in such 
a resistance is enough for 
the boiling of the water. A 
quantity of water evaporates 
and the water path becomes 
narrower. Dry zones are 
formed, micro-discharges en- 
sue and finally the flashover 
follows. Such a phenomenon 
was observed in the inclined 
arrangement experiments, as 
they were observed before 
with non-inclined test ar- 
rangements [8] . 
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Conductivity (pS/crn) 



Fig. 12: Flashover voltage for various conductivities, 
positioning and volume of the droplets. PVC used 




Conductivity (pS/cm) 



Fig. 13: 
positioning 



Flashover voltage for various conductivities, 
and volume of the droplets. Silicone rubber used 
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Conductivity (pS/crn) 



Fig. 14: Flashover voltage for various conductivities, 

positioning and volume of the droplets. Rubber used 



Discussion and thoughts for 

further research 

In the present paper, some 
parameters influencing the 
droplet behavior on polymeric 
surfaces were investigated, 
such as water conductivity, 
droplet volume, polymeric 
surface roughness and droplet 
positioning. An increase of 
conductivity causes a de- 
crease of flashover voltage. 
This is a statement valid ir- 
respective of the polymer 
used. The surface roughness 
affects in a positive way the 
flashover voltage, when the 
number of droplets is large. 
The surface roughness func- 
tions as a hindrance to the 
movement of the droplets, and 
consequently renders their 
oscillation more difficult. 



An increase of droplet volume 
causes a decrease of flash- 
over voltage. This is in 
agreement with experimental 
observations published before 
with either ac or dc electric 
fields [13] . The position of 
the droplets with respect to 
the electrodes is of vital 
importance. With the droplets 
nearer the electrodes, the 
flashover voltage decreases. 
This is a phenomenon ob- 
served, albeit in different 
circumstances and conditions, 
also with enclosed cavities 
in solid dielectrics, where 
discharges become much more 
intense when one of the en- 
closing walls is an electrode 
[14]. 

The above show clearly 
that the polymeric material 
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plays a predominant role in 
determining the flashover 
voltage and the behavior of 
water droplets. Hydrophobic 
materials, such as silicone 
rubber, perform better than 
PVC or rubber. With this in 
mind, one should also note 
that most polymeric materials 
for outdoor applications pre- 
sent some sort of hydropho- 
bicity. However, the advan- 
tage of silicone rubber con- 
sists in the fact that it 
does not only have this prop- 
erty, it can also regenerate 
it [15]. 

The formation of water 
paths, between the droplets 
as well as between the drop- 
lets and the electrodes, gen- 
erally follow the direction 
of the applied electric 
field. The general activity 
in the form of discharges and 
droplet movement with rougher 
surfaces, sets in at higher 
voltages. In the case of just 
one droplet, with the appli- 
cation of the field, a defor- 
mation starts turning later 
to instability. Such behavior 
was observed with the in- 
clined arrangement as well as 
with previous horizontal ar- 
rangements [3], [4]. Also in 
the case of the inclined ar- 
rangement, the role of the 
'triple points' (i.e. the 
points where air, polymeric 
surface and droplet meet each 
other) is vital. The forces 
exercised on the droplets. 



because of the applied elec- 
tric field, are quite strong, 
and therefore, the 'triple 
points' move towards the 
electrodes. Experimental data 
published recently, validate 
what is reported here [16] . 
Such movement of 'triple 
points' causes the spread of 
the droplets. The spread of 
droplets is perhaps the most 
characteristic phenomenon ob- 
served with the inclined 
electrode arrangement. It is 
not, however, the only one 
observed. Droplet oscilla- 
tion, formation of water 
paths, collapsing of two 
droplets into a larger one, 
ejection of small charged 
droplets from a larger one, 
were also noted during the 
experiments. In this respect, 
the present work offers simi- 
lar conclusions with those in 
[3], [4], [8], [9], [16], [17], 
The importance of the triple 
points should be emphasized. 
In other works, it has been 
reported that partial dis- 
charge (PD) activity is some- 
times marginal, not easily 
detectable by a conventional 
electrical PD system but by a 
photomultiplier. Smaller wa- 
ter droplets offer higher 
flashover voltages. This can 
also be explained by consid- 
ering - in the case of smaller 
droplets - that the inter- 
molecular forces are in equi- 
librium with the surface ten- 
sion and, consequently, the 
electrical forces required to 
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disturb the equilibrium are 
higher. On the other hand, 
for larger droplets, the in- 
ner forces may be lower and, 
if the surface tension stays 
the same, the equilibrium can 
more easily be disturbed. In 
such a case, the PD activity 
will start earlier and at 
lower voltages, and there- 
fore, the flashover voltage 
will be lower [18] . The fact 
that the contact angle in- 
creases with the decrease of 
water droplet volume, is 
something that cannot be un- 
derestimated [18] . Generally 
speaking, although the obser- 
vations of [18] refer to sil- 
icone rubber samples only, 
the relevant conclusions are 
not different from the ones 
presented in this paper. 

It is to be noted that the 
inclined electrode arrange- 
ment, used here, should not 
be compared by any means with 
the well known arrangement of 
the inclined plane test [19] . 
In the latter, a film of 
electrolyte is arranged to 
trickle down the back surface 
of a sheet and the samples 
are rated in terms of the 
voltage which causes a track 
to form in one hour [20] . In 
other words, the inclined 
plane test is a means of 
evaluating resistance to 
tracking and erosion of insu- 
lating materials for outdoor 
use, whereas the inclined 
electrode arrangement used in 



this work is a setup to study 
some parameters affecting the 
droplet movement on polymeric 
materials. The inclined plane 
test is an accelerating test 
[19] . The angle which was 
used in our experiments, i.e. 
the angle of 10°, was taken 
from real insulators. The 
purpose was to see the drop- 
let behavior under an elec- 
tric field in, as much as 
possible, real conditions. 
The present work confirms 
some general tendencies noted 
in [3], [4], 

It would be interesting to 
try experiments, in which the 
contact angle will be accu- 
rately measured for both 
smaller and larger water 
droplets w.r.t. time. The 
purpose would be to see 
whether the contact angle di- 
minishes with the passing of 
time and which diminution is 
more dramatic, that of the 
contact angle of smaller 
droplets or that of the con- 
tact angle of the larger 
droplets. Recent research in- 
dicated that the contact an- 
gle of droplets decreases 
with time, without, however, 
precising whether the rate of 
decrease is larger for smal- 
ler or for larger droplets 
[21] . Moreover, as noted be- 
fore, an interesting point 
can be to study different 
modes of droplet deformation 
in terms of the four parame- 
ters investigated in this pa- 
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per [18], [22], 

One last point should be 
raised: the research men- 
tioned above was carried out 
with conventional polymers. 
It is remarkable that, most 
of the aforementioned points, 
i.e. the importance of the 
positioning of the droplets, 
the significance of the drop- 
let volume and/or of the 
droplet number as well as the 
importance of the water con- 
ductivity, are points which 
are also important for non- 
conventional polymers, i.e 
for nanocomposite polymers. 
Work done recently in this 
laboratory emphasizes that 
parameters such as those in- 
vestigated here, are also 
significant for nanocomposite 
polymers. In fact, the nature 
of phenomena such as surface 
discharges or flashovers, must 
have a common underlying de- 
nominator for both conven- 
tional polymers and nanocom- 
posite polymers. Future re- 
search must also be directed 
in order to find the common 
grounds for surface discharge 
phenomena in both conven- 
tional and nanocomposite poly- 
mers [23] , [24] . 

Conclusion 

Water droplet conductiv- 
ity, polymer surface rough- 
ness, droplet volume and the 
positioning of droplets with 
respect to the electrodes 
constitute important parame- 



ters affecting the behavior 
of droplets under the influ- 
ence of an electric field 
with an inclined plane elec- 
trode arrangement. Increased 
conductivity, smoother poly- 
mer surfaces and increased 
droplet volume cause a reduc- 
tion of the flashover volt- 
age. The droplet positioning 
with respect to the elec- 
trodes plays a vital role in 
reducing the flashover volt- 
age and, on occasions, is 
more important than the drop- 
let volume. 

Appendix - Elementary 

modeling 

As mentioned above, the 
behavior of a droplet was 
modeled in [25], where the 
electric field E N developed 
at one edge of the droplet is 
given by 

E N = U h/[a(h - a) ] (1) 

where, U is the applied volt- 
age, a is the droplet radius 
and h is the distance of the 
center of the droplet from 
one of the electrodes. The 
electric field E„ on the 
opposite edge of the droplet 
is given by 

E m = U h(L - h)/[a(L - h - a)] (2) 

where, L is the distance be- 
tween the electrodes and the 
other symbols as in Eq . (1). 
As a droplet is positioned in 
the middle of the electrodes, 
the ratio E N /E M is expected to 
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be unity, and this was what 
was exactly obtained with the 
above equations. Considering 
Eq . (1), as h^a, E N tends to 

infinity. This is what we 
observed in the context of 
this work, namely that the 
closer a droplet is in one of 
the electrodes, the larger 
the electric field is get- 
ting, and consequently the 
more deleterious the conse- 
quences are. The above simple 
modeling is due to [25] . It 
was elaborated in [26], 
albeit with a non-inclined 
test arrangement. It must be 
mentioned that phenomena of 
very similar nature were 
observed also with the in- 
clined test arrangement. 



Remark 

Main aspects of this paper 
were published in M.G. 
Danikas, P. Ramnalis and R. 
Sarathi, "A study of the be- 
havior of water droplets on 
polymeric surfaces under the 
influence of electric fields 
in an inclined test arrange- 
ment", Journal of Electrical 
Engineering, Vol. 60, No. 2, 
2009, pp. 94-99. The present 
version, however, contains 
additional comments on some 
aspects of the investigated 
polymeric materials as well 
as on a possible relation of 
the observed phenomena in 
nanocomposite polymers. 
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Abstract 



In order to demonstrate the usefulness of the only one 
existing method for systematic error estimation in VNA (Vec- 
tor Network Analyzer) measurements by using complex DERs 
(Differential Error Regions), we compare one-port VNA mea- 
surements after the two well-known calibration techniques: 
the quick reflection response, that uses only a single S 
(Short circuit) standard, and the time-consuming full one- 
port, that uses a triple of SLO standards (Short circuit, 
matching Load, Open circuit). For both calibration tech- 
niques, the comparison concerns: (a) a 3D geometric represen- 
tation of the difference between VNA readings and measure- 
ments, and (b) a number of presentation figures for the DERs 
and their polar DEIs (Differential Error Intervals) of the 
reflection coefficient, as well as, the DERs and their rec- 
tangular DEIs of the corresponding input impedance. In this 
paper, we present the application of this method to an AUT 
(Antenna Under Test) selected to highlight the existence of 
practical cases in which the time consuming calibration tech- 
nique results a systematic error estimation stripe including 
almost all of that of quick calibration. 



Introduction 

The systematic error in a 
full one-port calibrated VNA 
measurement p of a given one- 
port DUT (Device Under Test) 
is already estimated by its 
DER [1] - [2] : 

p = (m - D)/[M(m - D) + R] (1) 



dp = [- RdD - (m - D) 2 dM - 
- (m - D)dR + Rdm] / 

/ [M(m - D) + R] 2 (2) 

where m is the VNA complex 
reading and D, M and R are 
the complex system errors of 
Fig. 1. 

The relations holding be- 
tween this complex reflection 
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coefficient p and its respec- 
tive impedance Z, as well as, 
between their DERs are [1] - [2] : 

Z = Z 0 (l + p)/(l-p) (3) 

dZ = 2Z G dp/(l - p) 2 (4) 

In this paper, we express 
the DERs for systematic error 
estimation in VNA measure- 
ments calibrated by the much 
simpler and quicker reflec- 
tion response technique, in 
order to be in place to make 
some practical decisions from 
the different calibration tech- 
niques comparison. 



1 




R 



Fig. 1: Full one-port error 
model 

Response Calibration 

The reflection response ca- 
libration technique can be 
accomplished with the mea- 
surement of only one standard 
load, instead of three in 
full one-port, usually of a S 
short circuit [3] -[4]. This 
means that the flow graph of 
Fig. 1 is simplified a lot, 
since the two surrounded by 
dashed boxes system error 
branches of directivity D and 
source match M do not exist, 
equivalently D = 0 and M = 0 
and (1) results to: 



R = m/ps = s/S (5) 

where s is the VNA complex 
reading of the S short cir- 
cuit standard with a nominal 
value of S = - 1, m is the com- 
plex reading of a given DUT 
and ps is its complex reflec- 
tion coefficient as it is 
measured after this response 
calibration : 

p s = (m/s)S (6) 

which, from (2), has the dif- 
ferential error: 

dps = (S/s)dm - (Sm/s 2 )ds + 

+ (m/s)dS (7) 

The corresponding total DER 
is then the sum of L = 3 par- 
allelograms. Therefore, this 
DER contour is a polygonal 
line with 4L = 12 line seg- 
ments and vertices at most, 
in contrast with the DER of 
the measurement after a SLO 
full one-port calibration, 
which is a piecewise curve com- 
posed of 4(L - 1) = 24 line 
segments, 4(L - 1) = 24 circu- 
lar arcs and 8(L-1)=48 ver- 
tices, at most [1] - [2] . 

Application Results 

By following the error es- 
timation process, we already 
detailed in [l]-[2], we take 
as dS the considered manufac- 
turers' standard S uncer- 
tainty data: 

-0.01 < d | S | < 0, -2° < dzS < +2° 
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and as dm and ds the VNA in- 
accuracy of ±1 digit in LSD 
of their corresponding read- 
ings, for either the ampli- 
tude in decibels or the phase 
in degrees. Moreover, the one- 
port DUT that was considered 
is the same typical UHF ground- 
plane antenna (that is: AUT) 
mentioned in [1] . 

The difference between the 
3 nominal values (-1, 0, 1) 
of the 3 full one-port cali- 
bration standards (S, L, 0), 
respectively, and their 3 cor- 
responding VNA readings (s, 

1, o), can be estimated by 
the extent of the surfaces 
shown in the triptych of Fig. 

2, where the vertical axis 
segment represents the range 
of the distinct stepped fre- 
quencies. Each surface is for- 
med by parallel to horizontal 
plane lines. Each such line 
expresses the complex differ- 
ence between the standard no- 
minal value and its corre- 
sponding VNA reading, in each 
stepped frequency. 

In the triptych of Fig. 3, 
and from left to right we 
have the difference surfaces 
made by distance lines be- 
tween : 

(a) the measured reflection 
coefficient p after a full 
SLO one-port calibration 
(black solid points) and the 
corresponding VNA readings m 
for the AUT measurement (col- 
ored magenta points), 



(b) the measured reflection 
coefficient ps after S re- 
sponse calibration (black 
ring points) and the corre- 
sponding VNA readings m for 
the AUT measurement (colored 
magenta points), and 

(c) the two measurements (p, 
Ps) . 

All the involved, previ- 
ously shown, quantities are 
projected on the horizontal 
complex plane of Fig. 4. The 
magenta colored spiral repre- 
sents m, while, the black 
curves the reflection coeffi- 
cient: solid points, for p, 

and ring points for ps. All 
of 1 VNA readings are close 
enough to complex origin (co- 
lored green points). It is 
rather difficult to distin- 
guish the two curves for s 
and o VNA readings, which are 
close enough to the unit cir- 
cle circumference (colored 
red solid points and colored 
blue ring point, respec- 
tively) . 

The p-DERs and ps-DERs, 
for all 4 MHz stepped fre- 
quencies covering the range 
of [600, 1000] MHz, are over- 
lapped on the complex plane 
of Fig. 5, forming a light 
and a dark gray stripes, re- 
spectively. From each stripe 
we selected 11 DERs out of 
101, drawn with dark gray and 
white colors respectively, to 
illustrate their outline de- 
pendence on frequency. 
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1 V 1 .08 .08 1 1 

Fig. 2: Difference between s and S, 1 and L, o and 0 




Fig. 3: Difference between m and p, m and ps, p and ps. 
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Fig. 4: VNA s, 1, o, m readings and p, ps measurements. 



Moreover, we selected to 
magnify a part of this figure 
in the sub-range of [892, 
1000] MHz, to further illus- 
trate the DER outlines and 
their overlapping in Fig. 6, 
where the clearly shown rip- 
ple of the simple response 
calibration stripe over the 
relatively smooth full one- 
port calibration stripe re- 
veals the superiority of the 
latter in the production of 
more accurate measurements. 

The comparison between the 
AUT measurements based on 



these two calibration tech- 
niques is extended to the 
comparison against the fre- 
quency : 

(a) of the computed polar 
DEIs of the reflection coef- 
ficient magnitude and argu- 
ment stripes in Fig. 7, 

(b) of the rectangular DEIs 
for the corresponding R input 
resistance and X input reac- 
tance stripes, in Fig. 8 and 

(c) of the Z-DERs, and Zs- 
DERs stripes in Fig. 9. 



SUNDAY 15 DECEMBER 2013 



V2-45 



FUNKTECHNIKPLUS # JOURNAL 



A/. I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 




Fig. 5: Complex p-DERs and ps-DERs in [600, 1000] MHz 




Fig. 6: Complex p-DERs and ps-DERs in [892, 1000] MHz 
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Frequency, MHz 




Frequency, MHz 



Fig. 8: Rectangular DEIs of input impedance 
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Fig. 9: Complex Z-DERs and Zs-DERs in [600, 1000] MHz 



Conclusion 

From all that, it must be 
clear now that in this inten- 
tionally selected for presen- 
tation particular AUT case 
there was no advantage at all 



in selection of full one-port 
calibration over the reflec- 
tion response one, due to their 
remarkable in all aspects co- 
incidence. Of course this is 
just another one conclusion 
a-posteriori . 
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Antennas Research Group, Austria - Hellas [l, 2] 

EECE Dept, Democritus University of Thrace, Hellas [2] 

Abstract 

Based on the experience we have gained so far, as indepen- 
dent reviewers, we thought that may be proved useful to pub- 
licly share with the interested author, especially the young 
one, some practical implementations of our ideas for the in- 
teractive representation of data using 3D/4D movement and ani- 
mation, in an attempt to motivate and support him in the de- 
velopment of similar dynamic presentations, when he is look- 
ing for a way to locate the stronger aspects of his research 
results in order to prepare a clear, most appropriate for 
publication, static presentation figure. For this purpose, we 
selected to demonstrate a number of presentations, from the 
simplest to the most complicated, concerning well-known an- 
tenna issues with rather hard to imagine details, as it hap- 
pens perhaps in cases involving Spherical Coordinates and Po- 
larization, which we created to enrich the very first ever 
made Virtual Laboratories of Antennas, that we distribute 
over the Open Internet through our website Virtual Antennas. 
These presentations were developed in a general way, without 
using antenna simulators, to handle output text and image 
data from third-party CAS Computer Algebra Systems, such as 
the Mathematica commercial software we use or the Maxima 
FLOSS we track its evolution. 

Introduction 

The very first ever made 
Virtual Laboratories of An- 
tennas, which are available 
on the Internet through Vir- 
tual Antennas website [1], 
are exclusively based on our 
alternative form of Antenna 
Theory [2], and are founded, 
by applying the learning by 



teaching method to Antennas 
education [3], in order to 
support the quick comprehen- 
sion of 3D space matters un- 
avoidably related to anten- 
nas. The initial material of 
Virtual Antennas was prepared 
during the years 1996-1997 
and presented to Antenna stu- 
dents during the spring seme- 



WEDNESDAY 1 JANUARY 2014 



Vl—7 



FUNKTECHNIKPLUS # JOURNAL 




A/. I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 



ster of 1998 [4] . On 1999, 

the existence of this visual 
material was announced to the 
first ever appeared on the 
web EMLIB Electromagnetics 
Library, maintained those 
years by Jet Propulsion Labo- 
ratory of NASA [5] . On the 
same year, all of the mate- 
rial was approved for inclu- 
sion to MathSource repository 
of Mathematica [6] . 

Since 2000 and until to- 
day, the Virtual Antennas ma- 
terial is constantly in- 
creased and improved. More- 
over, a number of websites 
either use this material or 
suggest its use [7] . On 2009, 
a sample of its current de- 
velopment state was approved 
for publication by the Wol- 
fram Demonstrations Project 
[8], an event that was an- 
nounced in our Creative Com- 
mons Network web pages [9] . 

During the last years, our 
voluntary reviewing work re- 
vealed some needs of poten- 
tial authors -especially the 
young ones- related to the 
presentation of their re- 
search results, which perhaps 
may be confronted by similar 
to our interactive presenta- 
tion techniques, and thus we 
decided to present in this 
paper our continuing work for 
the Virtual Laboratories of 
Antennas . 



Presentation Development 

and Use 

The implementation of a 
presentation idea demands the 
expression of the theoretical 
idea formulation in the CAS 
language we select to produce 
output data in the appropri- 
ate file format for movement 
or/and animation (we would 
like to emphasize that in the 
case of Mathematica we use, 
while we abandoned the pack- 
age written by Novak [10], we 
still run the one written by 
Donley [11]). These files 
are: (a) WRL (WoRLd: Open by 

Web3D Consortium), a plain 
text file with a known struc- 
ture described in the VRML 
Virtual Reality Modeling Lan- 
guage [12] for a visual ob- 
ject, (b) a bunch of images 
in a selected non-destructive 
format; usually BMP, and (c) 
NBP text file (Note Book 
Player: Mathematica propri- 

etary). After that, the de- 
velopment of the presentation 
requires: (a) either a text 

editor to correct bugs or to 
add special features into WRL 
output file, (b) a graphics 
editor to align images or 
correct blemishes, as well 
as, a video editor to handle 
these images as frames of the 
final movie file in a se- 
lected format, usually AVI 
(Audio Visual Interleave: 
Microsoft proprietary) . It is 
worth to notice that accord- 
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ing to our experience, the 
non-computing time needed: 
(a) to design the idea imple- 
mentation, has been reduced 
by the years, from a couple 
of weeks at the beginning to 
just a couple of days now, 
and (b) to develop the pre- 
sentation, after a tedious, 
routine work of editing each 
frame, varies considerably 
according to the frame theme, 
e.g. from 1.5 minutes for 
that in Fig. 7, to 10 minutes 
of Fig. 10. 

The software needed for 
presentations to work under 
MS Windows is: a web browser 
(MS- IE 3 . 02r+) , a VRML viewer 
add-on to that browser (World- 
View 2.1), an AVI player 
(Mplayer2), and the Mathema- 
tica 7 Notebook Player. 

3D/4D Presentation Samples 

In order to abbreviate the 
figure descriptions in the 
following, we have to notice 
that: (a) for every animation 
presentation, a typical full- 
window screen capture of the 
recommended AVI player is 
shown after its pause button 
was pressed, (b) for every 
movement presentation, a ty- 
pical full-window screen ca- 
pture of the recommended on- 
line WRL viewer, with a cur- 
sor perhaps to indicate the 
existence of an additional 
feature, and (c) for both 
presentation types, the basic 
colors (R, G, B) are used in 



that order to correspond ei- 
ther to the CCS Cartesian Co- 
ordinate System (x, y, z) 
axes, (xoy, yoz, zox) planes, 
and its unit vectors, or to 
the SCS Spherical Coordinate 
System (r-radius, theta-semi- 
circle, phi-circle) curves, 
(r-sphere, theta-semi-cone, 
phi-semi-plane) surfaces, and 
its unit vectors, while the 
same color correspondence 
holds for the radiation pat- 
tern cuts by the mentioned 
CCS planes or SCS surfaces. 

After that, short descrip- 
tions of the samples, are 
following, while all of these 
presentations will be always 
available in authors' group 
repositories in Virtual An- 
tennas [13] and Google Code 
Repository [14] . 

Fig. 1 shows 1 frame out 
of 12 of an AVI repeated-for- 
ever animation for the 3D 
plane-time presentation of 
the considered as time-har- 
monic sinusoidal current wave 
amplitude, on a 3 wavelengths 
portion of a long thin wire 
loop antenna terminated on 
some complex impedance. The 
propagated current wave p is 
decomposed in the dual cou- 
ples: (incident wave i, re- 
flected wave r) and (standing 
wave s, transmitted wave t), 
while the imposed letters and 
arrows on the figure indicate 
these current waves as well 
as the direction of their mo- 
tion . 
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Fig. 2 shows a screen-shot 
of a WRL presentation with 5 
additional predefined view 
points of 12 frames anima- 
tion, as the cursor indica- 
tes, of a 4D space-time asyn- 
chronous movement, from a 
view-point in the lst-octant. 
In essence, this presentation 
is an alternative geometric 
meaning of the considered as 
standing wave current on an 
adequate length of a thin 
wire dipole antenna. 

Fig. 3 shows a screen-shot 
of a WRL presentation of 4D 
space-time asynchronous move- 
ment from a view point in the 
lst-octant, of the 3D differ- 
ence volume element formed by 
adjacent SCS coordinate sur- 
faces and lying between the 
1st and 2nd octant. 

Fig. 4 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment from a view-point in the 
lst-octant, for the triples 
of SCS unit vectors in the 
shown directions, that is 13 
triples on the CCS main- 
planes, as well as, 2 unde- 
fined triples on the irregu- 
lar z-axis of the CCS-to-SCS 
transformation . 

Fig. 5 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment from a view-point in the 
lst-octant, for all 3 SCS 
surfaces, curves, and unit 
vectors in a direction of the 
1st- octant. 



Fig. 6 shows 1 frame out 
of 37 of an AVI animation for 
the 4D plane-time presenta- 
tion of a SCS phi-circle cur- 
ve as the intersection bet- 
ween 2 SCS surfaces, that is 
of a r-sphere and a theta- 
semi-cone, while theta angle 
-which changes from 0° to 
180°- has the value of 50°. 
Different colors illustrate 
the inner and outer surfaces 
of the r-sphere and a cut off 
spherical section between ad- 
jacent meridians permit us to 
see the interior of that 
sphere . 

Fig. 7 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment for three different po- 
larizations at three differ- 
ent points: linear, circular 
and elliptical on each of the 
coordinate axes x, y and z 
respectively to emphasize 
that, in general, the pola- 
rization of an antenna is not 
constant . 

Fig. 8 shows 1 frame out 
of 20 of an AVI animation for 
the 4D space-time presenta- 
tion of the non-linearly CCW 
polarized time-harmonic real 
electric field E, from an an- 
tenna, as well as, its decom- 
position to two linearly po- 
larized time-harmonic real 
electromagnetic fields, as 
they are uniquely defined by 
the constant of time linearly 
independent space vectors E c 
and E s . 
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Fig. 1: Current waves 
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Fig. 2: Standing current- 
voltage waves relation 
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Fig. 3: SCS: Volume element 
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Fig. 9 shows the last fra- 
me out of 23 of an AVI anima- 
tion for the 4D space-time 
presentation of the path 
traced by the electric far- 
field arrow tip, which is CW 
elliptically polarized in the 
direction of y-axis. 




Fig. 7: Polarization: 
Directional dependence 




Fig. 8: Polarization: 
Non-linear CCW 



Fig. 10 shows 1 frame out 
of 19 of an AVI for the 4D 
space-time presentation of a 
normalized CCW elliptically 
polarized electromagnetic far 
-field in the shown direction 
of propagation, from a point 
of view in the 1st octant of 
an observer mirrored to (r, 
theta) plane. 




Fig. 9: Polarization: 
CW elliptical 
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Fig. 10: Polarization: 
Far-Field CCW elliptical. 
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Fig. 11 shows 3 frames out 
of 73 of an AVI for the 4D 
space-time presentation of 
the direction dependency of 
antenna polarization from a 
center fed short crossed 
dipoles on y and z axes, 
which have a current phase 
difference of 90°. The three 
different polarizations are 
shown with their orientation 
for an observer on the xOy 
plane (theta = 90° and in 

three different phi angles of 
0°, 240°and 270° . 

Fig. 12 shows 4 frames out 
38 of an AVI animation for 
the 4D space-time presenta- 
tion of the space radiation 
pattern and its three main- 
planes cuts, for the specific 
case of an antenna consisting 
of 2 dipoles, with the same 
direction of the unit vector 
(0.2, 0.4 0.894), each of 

2.4A length long, with their 

centers placed 0.25A apart on 
an axis with a unit direc- 
tional vector (0.3 0.5 

0.812), fed with a 30° cur- 
rent phase difference. 

Fig. 13 shows a screen- 
shot of a WRL presentation 
with 1 additional predefined 
view point of 12 frames ani- 
mation possibility, as the 
cursor indicates, of a 4D 

space-time asynchronous move- 
ment, from a view-point in 
the lst-octant. This presen- 
tation concerns the interior 
of an anechoic chamber and 

imitates a discone antenna 



rotation in order to visual- 
ize the way by which the mea- 
surement of a plane-cut of 
its radiation pattern is ac- 
complished . 
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Fig. 11: Polarization: CW 
Circular - CCW Elliptical - 
Linear 
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Fig. 14 shows a sophisti- 
cated combination of a move- 
ment with an animation pre- 
sentation, which in fact is 
an application that was de- 
veloped using the version 7 
of Mathematica, in a way that 
permits the definition of any 
(theta, phi) direction in 
space for a dipole of vari- 
able length. Notably, this 
application has all the avai- 
lable VRML movement features, 
with 6 predefined view 
points, while 3 different an- 
imations can be run simulta- 
neously. In addition, the 



mesh, opacity, rendering goal, 
and the evaluation step for 
the radiation pattern can be 
also defined. 

Fig. 15 shows the last 
frame out of 100 of an AVI 
for the 4D space-time presen- 
tation of a detailed study 
for some of the linear dipole 
characteristics regarding the 
plane and space radiation 
patterns, the direction from 
the dipole axis of their 
first common maximum, the di- 
rectivity, and the input ra- 
diation resistance; this re- 
ally is a 4-tuple of frames. 




Fig. 12: Radiation Pattern: 3D and its 2D main-plane cuts 
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Fig. 13: Anechoic chamber - discone antenna 




Fig. 14 Radiation Pattern: A composite 3D/4D presentation 
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Centre -Fed Linear Dipole with Standing Wave of Current 




normalized E Radiation Pattern 



ft . Plane Pattern 



B . Space Pattern 
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Fig. 15: Antenna characteristics 



Conclusion 

The Antennas Research Group 
(ARG) carries on for the en- 
hancement of this low-cost 
Virtual Laboratory with new 
material placed to web site 
www.antennas.gr [13], as well 
as to Google Code Repository 
[14], and Fig. 13 and Fig. 14 
are only two of many more 
possible examples. 

The evolution on computer 
systems in recent years has 
increased the demand for new 
teaching, learning and re- 



searching methodologies, which 
this paper supports by de- 
scribing available material 
for a web-based Virtual Labo- 
ratory for Antennas that pro- 
vides interactive computer 
aided visualization in vari- 
ous Antenna Theory topics. 

The main goal is con- 
stantly attained since this 
continuously developed mate- 
rial cannot be produced by 
the means of the conventional 
literature . 
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Abstract 



We have developed two methods to study an inverse 
scattering problem of electromagnetic waves in chiral media, 
for a spherical scatterer perfect conductor. The first 
called Far-field inverse problem was described analytical in 
our previous work. The second called Near-field method is 
given in detail in the present paper. It is a geometrical 
method based on the scattered fields. Using Near-field 
experiments, in which the scattered field is measured at the 
source, we solve the corresponding inverse scattering 
problem that is to determine the coordinates of the center 
and the radius of the spherical scatterer. 



Introduction 

We study an inverse scat- 
tering problem of electromag- 
netic waves with time harmo- 
nic dependence applying the 
method of Near-field [ 1 ] - [ 3 ] . 
Spherical electromagnetic wa- 
ves generated by a point 
source incident on a spheri- 
cal scatterer perfect condu- 
ctor in a chiral media. 

Knowing the incident and 
scattered wave fields in the 
inverse problem we are 
looking for the coordinates 
of the sphere center and its 
radius . 

The corresponding problem 
is solved by geometrical me- 
thods using either energy 
scattering cross-section, or 



the scattered fields. In the 
first case we refer to the 
inverse Far-field problem [4] 
and in the second case to a 
Near-field inverse problem. 

The chiral materials exhi- 
bit the phenomenon of optical 
activity i.e., the phenomenon 
that in which the plane of 
polarizations of linearly po- 
larized light is rotated as 
the light passes through an 
optically active medium. 

The chirality is a proper- 
ty which is often found in 
nature and reflects the asym- 
metry in the spatial inver- 
sion. An electromagnetic wave 
into a chiral media is analy- 
zed in a counter-clockwise 
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( LCP ) and a clockwise (RCP) 
Beltrami field. Using vector 
spherical harmonic functions 
and appropriate expansions of 
the Beltrami fields and by 
extension the development of 
spherical electromagnetic wa- 
ves in spherical wave fun- 
ctions, we calculate the 
exact solution of the scat- 
tering problem for a spheri- 
cal perfect conductor in chi- 
ral media, and the correspon- 
ding back scattering [3] -[4]. 

In the inverse problem we 
calculate the scattered Near 
-field, when the source is in 
position r = r 0 . A similar pro- 
blem for achiral materials 
has been studied in the works 
[ 1 ], [ 2 ]. 

Statement of the problem 

Consider a point source at 
the position r 0 that produces 
spherical electromagnetic wa- 
ves in a chiral media near a 
scatterer Q” , perfect con- 
ductor, i.e. the surface of 
the boundary condition is sa- 
tisfied 

n x (E inc (r) + E sc (r)) = 0 

( 1 ) 

r = a 

We consider that the exte- 
rior Q + = R 3 /Q of the scat- 
terer is homogeneous chiral 
media with fixed chirality 
(3, dielectric constant s and 
magnetic permeability p . 



A spherical incident ele- 
ctromagnetic wave (E“ c (r), 

H^ nc (r)) with time harmonic de- 
pendence in accordance with 
the Bohren transformation, 

analyzed in spherical Beltra- 
mi fields Q^(r) and Q^(r) 
as follows : 

E;; c (r) = Q^(r) + Qj^(r) 

, H ^(r) = i-(Q^(r)-Q^(r)) (2) 

where q = (p/s) 1/2 is the in- 
trinsic impedance of the chi- 
ral medium. The Beltrami 

fields satisfy the equations, 

[5], [6], 



JVxQ 4re (r) = y,Q 4r0 (r) 
[ Vx Q/?,r B ( r ) = -y fi Qi., re (r) 


(3) 


where y L , y R are wave 


numbers 


for Beltrami fields 
given by, 


and are 


k _ k 


(4) 


Yl ~ l - fcp ' Vr ~ l + fcp 



With /c = «(ep) 1/2 , co being the 
angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form, 
as defined in the following 
works: [4] (issue 2, p. 9, 

relations (4), (5)), and [3]. 
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If Q^(r) and Q^r(r) are 

Beltrami fields corresponding 
via transformation to Bohren 
in Ef c and H“ c , then the 

scattering problem for the 
perfect conductor can be 
formulated in the following 
way: be found Q™(r), Q™(r) , 
which belong to the space 
C 1 (Q + ) n C(Q + ) , such that: 



(E* c (r) , H* c (r)) . The scattered 
electric field will be depen- 
ded on the polarizations p L , 
p R , (see [4] (issue 2, p. 
9, relation (6)) and will have 
the decomposition 

Er B c (r|p L ,P R ) = Q^ e (r|p L ,p R ) + 

sc . . ( 6 ) 

+ Q R ,r 0 (r I P„ P R ) 




(iii) : 



r x Q“,(r) + iQ£. (r) = o( i) 
f X Q“.(r) - iQ~,(r) = o(i) 

r -a oo 



The limits on radiation 
conditions (5(iii)), are ta- 
ken uniformly in all dire- 
ctions reS 2 , where S 2 is 

the unit sphere in M 3 and n 
is the outward normal unit 
vector perpendicular to the 
surface on the scatterer. 

The incident electromagne- 
tic wave (E“ c (r), H“ c (r)) on the 

scatterer generates the 

corresponding scattered field 



Where Q“ e (r|p t , p„) and 

Q“ r . O' I P^ P R ) are the corre- 
sponding scattered Beltrami 
fields which have the fol- 
lowing behavior, when r oo , 
[5], [7], 

Q“r.(Mp t , P R ) = h 9 (Y A ry 

. . . 1 ( ? ) 
•g*,r e (r|P^P R ) + 0(pr) 

with A = L, R . The functions 
g L and g R are the LCP and 
RCP far-field patterns respe- 
ctively, which are defined by 
the following relationship 
[7], 

g*,r 0 (f |P„P R ) = 

Jnx[y^VxE; c (r'|p L , p R ) - (8) 

S 

- TO„y 2 E* c (r' | p„ p R )e-^ r r ds(n 
f-1, A = L 

with m A =\ , y 2 = y L y R . 

{ 1 ,A = R 
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Near-field Inverse problem 

Using the exact solution 
found in the paper [4] (issue 
2, pp. 10-12) and considering 
appropriate asymptotic forms 
of the Bessel and Hankel 
functions in low frequencies, 

i.e. \y A a\ « 1 , calculate asym- 
ptotic expressions of the 
coefficients a L n , a* , b L n and 
b R . Specifically we have [3], 



a L i + P k 

n 2i^(2n + 1) 
y L a -+ 0 

p i 



(y L a ) z 



( 9 ) 



(1 - (3/c) n + 2 , , 2 

E7i(YL a ) 



2 < (i + m 



or 



bt 



i (1 + (3/c) n + 2 
2< (1 - p/c) n + 1 



(y R a) 2 " 



1 



y R a -+ 0 



(10) 



b R n 



i(l-p/c) 
2^(2 n + 1) 



(WO 



2/7 + 1 



Where 

<■„ = l-3-5 - (2n - 1) = (2n)!/(2"n!) 



In the inverse problem we 
calculate the scattered Near- 
field in the source r = r 0 . By 
the relation (21), in work 
[4], with r = r 0 and for LCP 
incidence, we obtain [3], 



E r„ c ( r 0 l PJ = 



=z 



n ( n + !) d l. 



B n a n {(h n {y L r g )x + h n {y L r B ) ip) + i{h n {y L r B )x - h n {y L r B ) ip)} + (11) 



00 n( n + 1) - 

+ Z 2 e n a n{(^n(yA)x - h n {y R r B ) ip) + i (~h n {y R r B )x - h n (y R r B ) ip) } 



Using the asymptotic rela- with y A r B -> 0 , A = L, R and the 



tions : 




series. 


ft n (>V a ) 


Zn 




00 


h e (Y A r B ) 
fy,( Ya O 


(y A r B y 

~nZ n 


< 

(12) 


^.2/7 + 1 

n-1 

CO 

V 1 2 2/7 + 


h B (y A r Q ) 


(Y A rJ + 1 


2, n r 

n=l 


h n (y A a ) ~ 


"Zn 




oo 


i(y,a) n + 2 


X (2n + 1 



(1 - r 2 ) 3 

, + 1 t(3t 2 - r 4 ) 

(1 - r 2 ) 2 



(13) 



(14) 
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n{2n + 1 )t 2 



hl t(3t 2 + r 4 ) 

= (1 - T 2 ) 3 



y 2n + 1 _^ 2 n + 1 
n 



i-2\ , o-i-2 



T - 1 

and because of the relation- 
ships : 



(15) 


E "r e ( r o 1 Pa) ~ 










(1 - gj 4 (3/c) t 5 


(20) 




2(y,a) 2 (1 - r 2 ) 3 


(16) 


y A a 0, A = L, R 




■T 2 )] 


Choose a Cartesian coordi- 


nate system Oxipz, and 


five 



1 - V k 
Y l 1 + P*c 



, 0 < r < 1 , 



where r = a/r 0 , we obtain from 
(11) [3], 



Er 0 c (r 0 I Pj ~ 

(1 + P/C) r B 



2 (y L ay (l-iy 



Pl 



(17) 



point-source locations, name- 
ly 0(0, 0, 0) , A 1 (l, 0, 0) , A 2 (0, 

1, 0) , and /A 4 (0, 0, 21) , which 
are at (unknown) distances 
r 8 , r 4 , r 2 , r 3 and r 4 , respe- 
ctively from the sphere's 
center K. The parameter 1 is 
a chosen fixed length. The 
sizes of the resulting five 

measurements |E^ r (rv | p 4 )| are 



M < 21 > 



Therefore the measure of the 
scattered field E* c (r 0 | p L ) is 

|e sc (r | 6 )| ~ ^ + ll 

| r ^° IPJ| 2 (y L a) 2 (1-t 2 ) 3 (18) 

y L a 0 



Similarly for the RCP inci- 
dence, we have that 



Er> 0 I P R )| 



(1 - (3/c) 



2(y*a) ( 1 -t) ( 19 ) 
y R a -a 0 



Finally, based on the rela- 
tionships (18) and (19), we 
conclude that 



and 




Where j = 0, 1, 2, 3, 4 . Therefo- 
re we have the following five 
measurements 






b 5 Pj 

(p- - by 



(23) 



with p J >h 2 >0. Applying the 

law of cosines in triangle 
K0A 4 and using that 

rj=Pjl (24) 
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we obtain 



P 4 = 2 + 2 P3-Pg ( 25 ) 

Because of (24), we take 

p] - 3 p)b A + 3 p)b s - 

p ( 26 ) 

-b 12 -b 5 ^- = Q 
Mj 



a is too small so that a «: 1 , 
then the relationship (23) 
shows [3], 




so we have the following sy- 
stem of equations 



Thus in equation (26) iden- 
tified the p , with j = 0, 1, 2, 

3, 4 . If the radius a of the 
sphere is known, solve the where j = 0, 1, 2, 3, 4 . 
system of six algebraic equa- 
tions (25) and (26). Also if 
you consider that the radius 



p 4 = 2 + 2p 3 - p G 
p 5 1 M 1 - b 5 = 0 



(28) 
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Interfaces in High Voltage Engineering: 

A Most Important Question for Conventional 
Solid Insulating Materials as well as for 
Nanocomposite Polymers 
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Abstract 

Interfaces consist a most important part of conventional 
insulating systems at high voltages. They are considered to 
be problem areas which have to be dealt with. Numerous pub- 
lications have contributed in rendering the mechanisms of 
interfaces understandable. On the other hand, interfaces in 
nanocomposite polymers seem to function in an entirely dif- 
ferent manner from that in conventional insulating systems. 
The present paper reviews past work on both the conventional 
insulations and in nanocomposites. Differences regarding the 
interfaces are mentioned and discussed. Whereas interfaces 
in conventional insulating systems are to be avoided, inter- 
faces in nanocomposite polymers seem to be desirable - at 
least - up to a certain percentage of nanoparticles in the 
base polymer. Although things are better understood in con- 
ventional insulating materials, more work has to be per- 
formed in order to clarify several aspects, such as space 
charges and electrical trees emanating from enclosed cavi- 
ties. Needless to say that much more work has to be done in 
nanocomposites w.r.t. their modeling and possible explana- 
tions of the surprising performance of interfaces, a perfor- 
mance that deviates strongly from the performance in classi- 
cal insulating materials. 
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Breakdown, breakdown strength, conventional insulating 
materials/systems, nanocomposite polymers, partial dischar- 
ges, electrical treeing/trees 
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Introduction 

With some sort of exag- 
geration, a well known pro- 
fessor said once that "the 
problems of high voltage in- 
sulations are problems of in- 
terfaces" [1] . Interfaces re- 
sult when there are two dif- 
ferent insulating materials 
next to each other or when an 
insulating material meets a 
conductor. An interface may 
become, e.g., the source of 
partial discharges or even the 
cause for a complete break- 
down of an insulating system. 
When two insulating materials 
of different dielectric con- 
stants have a common inter- 
face, then the material with 
the lowest dielectric con- 
stant will undergo the more 
intense electric stressing 
[2] . Depending also on other 
parameters, such stressing 
may result in a gradual dete- 
rioration of the insulating 
system and consequently in a 
complete breakdown [3] . 

In this paper, the ques- 
tions raised by the existence 
of interfaces in high-voltage 
insulating systems are dis- 
cussed. Interfaces play an 
important role in determining 
the robustness of an insulat- 
ing system, when conventional 
insulating materials are 
used. Interfaces play also a 
vital role in insulating sys- 
tems with nanocomposite poly- 
mers. Their functioning, how- 
ever, is of another nature. 



This short review is by no 
means exhaustive since the 
topic of interfaces is a vast 
one. It is the aim of the au- 
thors to give the gist of the 
problems and questions the 
researchers may face and to 
offer some comments. In the 
context of this paper, the 
terms "insulating material" 
and "dielectric material" or 
simply "dielectric" are some- 
times interchanged meaning 
the same thing. 

Interfaces in Conventional 

Insulating Materials 

Dielectric breakdown in 
insulating materials depends 
on electrode configuration, 
insulating material thick- 
ness, electrode materials, 
presence of microcavities, 
temperature, pressure, nature 
and morphology of material 
under test, type of applied 
voltage, damage path (surface 
or volume) [4] . Various di- 
electric breakdown theories 
have been put forward [5]- 
[13] . No matter whether the 
proposed theory was based on 
cumulative impact ionization 
by electrons - creating thus 
positive space charges which 
distort the field distribu- 
tion and weaken the dielec- 
tric [5] -, on the notion of 
"intrinsic breakdown" - ac- 
cording to which a large num- 
ber of electrons trapped in 
energy levels due to lattice 
imperfections can transfer 
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energy to the lattice vibra- 
tions [6] - [8] , on the "40 ge- 
nerations avalanche theory" 

[9] , on the importance of spa- 
ce charges which modify the 
local electric field value 

[10] , on the right assumption 
that the breakdown is a prop- 
erty of the dielectric mate- 
rial plus its electrode sys- 
tem [11] or on the theory 
based on the ionizing elec- 
trons and the hole traps 
[12], [13], the fact remains 
that all the above mentioned 
phenomena result from elec- 
tric field intensifications, 
i.e. from either electrode 
imperfections or mismatch of 
dielectric materials. This 
brings us to the point men- 
tioned in the Introduction of 
the present paper: that in- 
terfaces may create the con- 
ditions which may cause elec- 
tric field intensifications. 

The subject of partial 
discharges (PD) which may en- 
sue because of electric field 
intensifications and/or be- 
cause of gas (or foreign par- 
ticle) inclusions in a solid 
dielectric material, has been 
studied by Mason in his fun- 
damental publications [14]- 
[18] . Having in mind all the 
above, it is fitting to say 
that interfaces - created ei- 
ther by a mismatch of dielec- 
tric materials or because of 
intrusion of foreign parti- 
cles and/or air cavities in 
the insulating material under 



question -, are the problem 
areas of an insulating sys- 
tem . 

As was pointed out quite 
early [19], interfaces play a 
most significant role in the 
discharge and breakdown pro- 
cesses: even if an insulation 
does not contain any cavi- 
ties, at a sufficient stress 
"some event" releases gas 
with the subsequent formation 
of a cavity. The cavity is 
occupied by a gas discharge 
which increases the rate of 
gas formation with the subse- 
quent growth of both the cav- 
ity size and the discharge 
intensity. The importance of 
the differing nature of in- 
terfaces was also stressed in 
another publication, where it 
was pointed out that damage 
in internal cavities in poly- 
ethylene is little compared 
to the electrode adjacent ca- 
vities of the same dimensions 
and tested under the same ex- 
perimental conditions [20] . 
It is evident that in [20], 
interfaces between polyethy- 
lene and gas were compared 
with interfaces between poly- 
ethylene and electrodes, and 
the latter were found to be 
more dangerous and deleteri- 
ous to the insulating mate- 
rial. On the same lines, 
Kreuger showed that with PVC- 
insulated cables, the number 
of discharges increased with 
increasing electric stress in 
the dielectric [21] . In yet 
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another paper, it was indi- 
cated that the nature of in- 
ternal discharges was greatly 
affected by the assembly of 
the electrode system and the 
adhesion of the insulating 
tapes [22] . Discharges always 
start in the electrically 
weaker insulating material, 
as was commented in previous 
works [23] - [25] . 

Needless to say that phe- 
nomena related to PD, such as 
electrical treeing, are also 
closely connected to the mis- 
match of the dielectric con- 
stants of insulating materi- 
als and/or to the existence 
of gas cavities in their vo- 
lume. Earlier papers indi- 
cated that the treeing phe- 
nomenon in polyethylene ca- 
bles started from both inner 
and outer surfaces and also 
from solid particles and fi- 
bres [26] . Pioneering work 
performed with 15 kV and 22 
kV polyethylene insulated ca- 
bles reported that trees ori- 
ginated from contaminants and 
cavities, the tendency for 
tree initiation from a con- 
taminant being probably more 
affected by the contaminant 
material than by the size, 
location or shape of the con- 
taminant [27] . The importance 
of enclosed cavities in the 
initiation of trees was also 
reported more recently [28] . 
According to Ieda [29], tree 
propagation can be induced by 
internal gas discharge in the 



tree. It is to be bore in 
mind that in numerous publi- 
cations dealing with experi- 
mental work, the electrode 
arrangement that was used was 
a needle-plane electrode ar- 
rangement, indicating again 
that an electrode arrangement 
was chosen, with pronounced 
interfaces, in order to study 
the treeing phenomenon [30]- 
[33], 

Interfaces which may play 
a role in determining the 
breakdown strength of an in- 
sulating material need not be 
only interfaces between insu- 
lating material and metal or 
between insulating material 
and gas cavity or contami- 
nant. Differing phases may 
play also a role, as was no- 
ted in [34], where the inter- 
facial domain of crystalline 
and amorphous phases may de- 
termine the various proper- 
ties of semi-crystalline po- 
lymers, such as biaxially ori- 
entated polypropylene (BOPP). 
Stressing the importance of 
interfaces and experimenting 
with cross-linked polyethy- 
lene (XLPE) , McKean showed 
that a considerable improve- 
ment in cable breakdown can 
be achieved by impregnation 
with silicone oil or diethy- 
leneglycol. Such liquids can 
impregnate gas microspaces in 
the main insulation and thus 
result in an increase of the 
breakdown strength [35] . Si- 
milar observations were also 
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reported with polypropylene 
and polyethylene impregnated 
with suitable dipolar liquids 
[36]. 

Interfacial breakdown was 
studied with various elec- 
trode geometries and insulat- 
ing systems consisting of pa- 
per typical for transformers 
and transformer oil [37] . It 
was reported that interfacial 
breakdown will occur if the 
paper is not carefully dried 
or if many gaseous micro- 
porosities are left in or on 
the paper. In [37], however, 
it was also noticed that us- 
ing a carefully prepared pa- 
per-oil interface structure, 
the breakdown does not neces- 
sarily take place at the in- 
terface. Similar observations 
were made more recently by 
using silicone rubber inter- 
faces, where both perpendicu- 
lar and parallel to the ap- 
plied electric field were in- 
vestigated [38] . 

Conventional paper-oil ca- 
ble insulation was studied 
quite early and the problems 
of interfaces were noted 
[39] . Alternative insulating 
systems, based mainly on po- 
lymeric materials, were pro- 
posed with considerable com- 
mercial success [40] - [43] . 
Modern cables with solid 
polymeric insulation did not 
avoid the problems of inter- 
faces, namely those of extru- 
sions of semi-conducting 
sheaths with the main insula- 



tion or the inclusion of mi- 
crocavities and/or impurities 
[44], [45]. Extruded cable in- 
sulation exposed to wet con- 
ditions suffered from elec- 
trochemical treeing and impu- 
rities greatly deteriorated 
its electrical performance 
[46] . Moreover, operating elec- 
trical stresses may also cau- 
se premature insulation fai- 
lure in 15 kV polyethylene 
cables, if combined with un- 
favorable interface profiles 
and moisture [47] . Interfaces 
between polyethylene and 
small contaminants or micro- 
cavities may cause bow-tie 
trees in polyethylene cables 
[48]. 

Interfaces either perpen- 
dicular to the applied elec- 
tric field or parallel to it 
or at an angle with it were 
dealt with in [49], where it 
was noted that such a variety 
of interfaces may be encoun- 
tered in applications, such 
as capacitors, cables and in 
transformer windings. Compo- 
site insulating systems must 
preserve low dielectric los- 
ses. Higher dielectric losses 
may imply high ionic concen- 
tration in a solid/liquid in- 
sulating system, i.e. high 
ionic concentration in the 
liquid component of such a 
system [50] . 

Before concluding this 
section, it is fitting to men- 
tion the composite insulating 
systems of electrical machi- 
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nes, which consist mainly of 
epoxy resin and mica sheets. 
Previous work done in this 
context indicated that 
electrical treeing propagates 
through the epoxy resin and 
generally stops at the mica 
sheets, as mica is harder and 
electrically stronger than 
epoxy resin. The importance 
of such interfaces was re- 
ported before using a needle- 
plane electrode arrangement 
[51], where experiments were 
carried out without and with 
a mica sheet inserted in 
epoxy resin (Figs. 1 and 2). 
Evidently electrical trees 
were propagating more easily 
in the case of absence of the 
mica sheet and with much more 
difficulty with mica sheet. 

Simulation work done re- 
cently showed that mica she- 
ets prevent electrical trees 
from reaching the opposite 
electrode [52] . The purpose 
of mentioning experimental 
results regarding the compo- 
site system of epoxy resin/ 
mica sheets is to show that 
the electrical trees propa- 
gate through the electrically 
weaker insulating medium. The 
simulation results indicate 
that even the slightest vari- 
ations of dielectric constant 
may cause the electrical tree 
growth. In other words, the 
simulation data indicate that 



local fluctuations of dielec- 
tric constant imply local - 
even microscopically minute - 
formations of interfaces, 
which in turn may mean local 
field intensifications, en- 
couraging thus the growth of 
electrical trees (Figs. 3 and 
4). Such observations w.r.t. 
the local variations of die- 
lectric constant have also 
been reported for polyethy- 
lene [53] , [54] . 

It is evident from all the 
above that interfaces in 
classical insulating systems 
seem to cause problems (pos- 
sible dielectric constant 
mismatch, PD, treeing pheno- 
mena and ultimately risk of 
ultimate insulation failure). 
Due attention should be paid 
in choosing the insulating 
materials for specific appli- 
cations and to the construc- 
tion of the composite insu- 
lating system. Too many 
things depend on the quality 
of the construction of the 
interfaces [55]-[58], too ma- 
ny things that cannot be ig- 
nored. Interfaces in tradi- 
tional insulating systems are 
considered as the weak as- 
pects of such systems. Keep- 
ing this in mind and without 
exaggerating, it is not far 
from the truth if we state 
that an insulating system is 
as good as its weakest inter- 
faces . 
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Fig. 1: Electrical tree propagation without the presence of 
mica sheet (applied voltage 28 kVrms, 50 Hz) (after [51]) 




Fig. 2: Electrical tree propagation with the presence of 
mica sheet. The mica sheet increases the propagation time of 
the tree (applied voltage 28 kVrms, 50 Hz) (after [51]) 
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Fig. 3: Simulated electrical tree propagation with one mica 
sheet. Needle-plane electrode arrangement used 




Fig. 4: Further expansion of electrical trees. 
Electrical trees stop at mica sheet. 
Needle-plane electrode arrangement used 
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It goes without say that 
this short review regarding 
classical interfaces does not 
by any means cover the whole 
subject and variety of solid 
insulating materials (for ex- 
ample, no mention in this pa- 
per was offered about the in- 
terfaces in outdoor polymeric 
insulators [59], [60] or in 
indoor polymeric insulators 
[61]). Both in the libraries 
and in the Internet, the in- 
terested reader may find pra- 
ctically tens of thousands of 
publications referring to the 
questions and problem areas 
of the solid insulating mate- 
rials and insulating systems. 
What this short review tried 
to do is to show that inter- 
faces, electric field inten- 
sifications, pre - breakdown 
phenomena (such as PD and 
electrical trees) and break- 
down mechanisms are all in- 
terwoven and interrelated. 
Having said that, the next 
question related to inter- 
faces, is whether they play 
the same detrimental role in 
the new generation of insu- 
lating materials, the nano- 
composite polymers. This is 
to be examined in the follow- 
ing section. 

Interfaces in Nanocomposite 

Polymers 

More than twenty years 
ago, nanocomposite polymers 
came to our lives [62] . The 
first application of nanocom- 



posites appeared in 1990, 
when Toyota Motor Corporation 
introduced nanocomposite ny- 
lon in their car industry 
[63] . Since that year, many 
car industries introduced the 
use of nanocomposite poly- 
mers. Use of nanocomposite po- 
lymers was noted in other in- 
dustries, such as in the op- 
tics and the electronics in- 
dustries as well as in the 
food industry. A seminal pa- 
per by Lewis gave the impetus 
for research also in the in- 
sulation branch [64] . 

For the electrical insula- 
tion, nanocomposite polymers 
are defined as conventional 
polymers in which particles 
smaller than 100 nm are added 
and dispersed in such a way 
that at the end one gets a 
homogeneous mixture [63] . The 
addition of such nanoparti- 
cles (the term "nanofillers" 
is also widely used) is being 
done in very small quanti- 
ties, usually less than 10 wt%. 
Nanocomposite polymers con- 
sist of three components: 

a) the base polymer (or poly- 
mer matrix), 

b) the nanoparticles (or na- 
nofillers) and 

c) the interaction zone (or 
interface zone) between the 
matrix and the nanoparticles 
[63], 

Regarding the polymers used, 
these may be either ther- 
moplastics, thermosettings or 
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elastomers. Nanoparticles may 
be classified w.r.t. their 
dimensions, and they can be 
distinguished as 

1) mono-dimensional (i.e. ex- 
tremely thin), 

2) two-dimensional (nanotu- 
bes) and 

3) three-dimensional (inorga- 
nic oxides) . 

The most usual nanoparticles 
for the purposes of electri- 
cal insulation are 

1) silica nanoparticles SiC> 2 , 

2) montmorillonite nanopar- 
ticles (layered silica), 

3) metallic oxides such as 
AI 2 O 3 , Ti0 2 , MgO and ZnO and 

4) carbon nanotubes. 

Nanocomposite polymers can be 
obtained in two types of 
structures, namely, 

(i) intercalated nanocompo- 
sites (formed when there is 
limited inclusion of polymer 
chain between the clay layers 
with a corresponding small 
increase in the interlayer 
spacing of a few nanometers 
and 

(ii) exfoliated nanocomposi- 
tes (formed when the clay 
layers are well separated from 
one another and individually 
dispersed in the continuous 
polymer matrix [65], [66]. 

As mentioned above, nano- 
particles are added and dis- 
persed in relatively small 



quantities in the base poly- 
mer (usually no more than 10 
wt%) . Since nanoparticles are 
smaller than microparticles 
(smaller by three orders of 
magnitude), their interaction 
with the surrounding polymer 
matrix is much greater [62] . 
The so-called interaction zo- 
ne is the main factor con- 
tributing in the improvement 
of the properties of the base 
polymer [67] . In the case of 
addition of nanoparticles in- 
to a polymer, the interfaces 
are far more numerous and far 
larger than in the case of 
microparticles. As the size 
of the added particles is re- 
duced, the interface becomes 
larger and larger. The dis- 
tance between the nanoparti- 
cles is also extremely small. 
It seems that interfaces de- 
termine to a great extent the 
properties of nanocomposite 
polymers . 

The size of nanoparticles 
and the distance between them 
is of the order of magnitude 
on nanometers. Such particles 
may interact with the polymer 
matrix both physically and 
chemically in the nanometer 
scale. This has as conse- 
quence the appearance of pro- 
perties that are somehow dif- 
ferent from those we already 
know in a more macroscopic 
scale [68] . In contradistinc- 
tion to the interfaces in 
classical insulating materi- 
als, and also to what we know 
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from classical high voltage 
textbooks, the improved insu- 
lating properties of the na- 
nocomposite polymers are due 
to 

a) the large surface area of 
nanoparticles, which creates 
a large interaction zone, 

b) the changes in polymer 
morphology because of the 
large interaction zone, 

c) the changes in the space 
charge distribution and 

d) a dispersion mechanism 
[69]- 

Both the size of the nanopar- 
ticles and the chemical pro- 
perties of their surface play 
an important role in deter- 
mining the electrical, ther- 
mal and mechanical properties 
of nanocomposites. Needless 
to say that the chemical com- 
patibility between the intro- 
duced nanoparticles and the 
polymer matrix is of para- 
mount importance for the ge- 
neral properties of the nano- 
composite [70] . 

One of the most signifi- 
cant characteristics of nano- 
composite polymers is the in- 
crease of their breakdown 
strength as the size of the 
added nanoparticles tends to 
extremely small values. This 
increase is not in agreement 
with the conventional wisdom, 
which suggests that as the 
number of interfaces in- 
creases, the breakdown stren- 
gth decreases dramatically 



[2], [3]. Nanocomposite poly- 
mers seem not to agree with 
what we already know for 
classical insulating materi- 
als or systems [71] . Diffe- 
rences in breakdown strength 
between conventional epoxy 
resin and epoxy resin with 
nanoparticles was reported in 
[72] . Such observations were 
also noted later, when six 
different materials based on 
epoxy resin with various with 
and/or micro- and nanoparti- 
cles of alumina/silica were 
tested. It was shown that 
epoxy resin with nanosilica 
particles was the most suit- 
able to obtain high values of 
breakdown strength [73] . 

Addition of the percentage 
of nanoparticles to epoxy 
resin up to a certain level 
favors the increase of break- 
down strength both with a.c. 
and d.c. voltage, as was no- 
ticed in [74] . Why nanoparti- 
cles act in such a favorable 
way, despite the numerous in- 
terfaces they create? The in- 
crease of the breakdown 
strength may be due to 

(i) the increase of the sur- 
face area of the interfaces, 
which somehow alters the be- 
havior of the polymer, 

(ii) the changes of space 
charge distribution inside 
the insulation structure, 

(iii) the dispersion mecha- 
nism, and 

(iv) the changing properties 
of the insulating material, 
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more specifically its volume 
resistivity, its tan5 and its 
dielectric constant. 

It is possible that the elec- 
trons moving in such a nano- 
composite polymer, loose 
their kinetic energy because 
of the nanoparticles. Since 
the distances between the na- 
noparticles are extremely 
small, the electrons cannot 
acquire enough speed so that 
they can contribute to the 
breakdown process. Consequen- 
tly, epoxy resin with nano- 
particles presents a higher 
breakdown strength than con- 
ventional epoxy resin [74] . 

Normally the introduction 
of particles in polymeric ma- 
terials has as result the in- 
troduction of defects and 
subsequently the worsening of 
its electrical properties. 
Nanocomposite polymers seem 
not to obey the above rule, 
as the mechanisms of conduc- 
tivity during the breakdown 
process are influenced from 
the applied electric field, 
the dielectric constant of 
the nanoparticles and their 
number. The combined effect 
of these factors is difficult 
to fully understand at this 
stage and we need more work 
[71] . Similar results were 
obtained with epoxy resin 
with nanoparticles of Ti 02 as 
it presented a much higher 
breakdown strength than con- 
ventional epoxy resin [75] . 



On the other hand, elec- 
trical treeing propagation 
was found to be easier in 
conventional polymers than in 
nanocomposite polymers [76], 

[77] . Even a small wt% addi- 
tion of nanoparticles affects 
in a positive way the elec- 
trical treeing resistance of 
the nanocomposite polymer. It 
seems that electrical tree 
propagation paths go through 
the base polymer and around 
the nanoparticles (experimen- 
tal evidence for this was 
presented in SEM photographs 
published in [76]). Conse- 
quently, the more the nano- 
particles in a polymer, the 
more difficult the formation 
of treeing paths. It seems 
that nanoparticles act as ex- 
tremely small barriers, thus 
preventing the easy growth of 
electrical trees. Electrical 
trees propagate through the 
base polymer (in other words 
through the polymer matrix) 
and not through the nanopar- 
ticles. In some cases the 
electrical trees stop at the 
nanoparticles and they do not 
progress any further. Such 
observations were made in 
simulation studies recently 

[78] - [80] . 

Further research showed 
that a small percentage in- 
troduction of nanoparticles 
into a conventional polymer 
may increase its resistance 
to electrical treeing [81] . 
It is interesting to note 
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that nanoparticles may func- 
tion as barriers preventing 
the tree growth even in 
minute quantities [82] . Ear- 
lier work indicated that as 
soon as the electrical tree 
tou-ches the nanoparticle, 
the physico-chemical proper- 
ties of it are such that very 
high energies are required in 
order to cause its deteriora- 
tion [83]. Although the lat- 
ter paper is an old one, it 
may give a clue as to why 
nanoparticles act as elemen- 
tary barriers and why they 
prevent (or they delay) elec- 
trical treeing. More re- 
cently, a similar argument 
was given by some Japanese 
researchers [84] . 

Loading (i.e. the percent- 
age of included nanoparticles 
expressed in wt%) plays also 
a vital role in determining 
the resistance to electrical 
treeing. More loading (i.e. 
more nanoparticles, that is 
more interfaces) implies a 
better resistance to treeing. 
This is probably because 
trees interact with many more 
nanoparticles and this delays 
their growth [85] . Another 
possible explanation was of- 
fered some years ago, where 
the authors proposed that in 
front of a tree a damage 
process zone is formed in a 
conventional polymer. Such a 
zone cannot progress easily 
when it meets nanoparticles 
[86] . Simulation data indi- 



cated that loading affects 
the tree growth. The nanopar- 
ticle size plays also a role 
in delaying tree propagation, 
smaller nanoparticles offer a 
better tree resistance than 
the larger ones [79] . In 
Figs. 5-8 simulation results 
regarding the loading of na- 
noparticles as well as the 
size of nanoparticles are 
shown. It is evident that mo- 
re loading makes tree growth 
more difficult. It goes also 
without say that smaller 
nanoparticles offer a better 
resistance to tree propaga- 
tion . 

There is no need to empha- 
size that there is also in 
the field of nanocomposite 
polymers a vast body of tech- 
nical literature, too vast to 
be mentioned here in the con- 
text of this paper. From this 
short review it is obvious 
that interfaces in nanocom- 
posites play an entirely dif- 
ferent role from the one they 
play in classical insulating 
materials. Why is this so? 
This may be because the phy- 
sics and/or chemistry somehow 
change in the nanoscopic 
world of such materials. In- 
terfaces become highly desir- 
able - at least up to a cer- 
tain percentage of added na- 
noparticles. The surface area 
of the nanoparticles is huge 
if compared with that of mi- 
croparticles for other high 
voltage applications. 
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Fig. 5: Simulation with epoxy resin filled with HO 2 nano- 
particles (loading of 2 wt%, nanoparticle diameter 100 nm) 




Fig. 6: Simulation with epoxy resin filled with Ti02 nano- 
particles (loading of 6 wt%, nanoparticle diameter 100 nm) 
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Fig. 7: Simulation with epoxy resin filled with TiC >2 nano- 
particles (loading of 6 wt%, nanoparticle diameter 200 nm) 




Fig. 8: Simulation with epoxy resin filled with Ti02 nano- 
particles (loading of 6 wt%, nanoparticle diameter 100 nm) 
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Models regarding the ex- 
planation of the functioning 
of nanocomposites were pro- 
posed as well as how the 
nanoparticles behave inside 
the base polymeric material 
[62], [87], [88]. Such models 
try to explain the higher 
breakdown strength values of 
nanocomposites and also their 
higher resistance to electri- 
cal treeing. The explanations 
seem to be plausible but more 
hard evidence is needed. Such 
evidence will be offered by 
many more photographs (SEM 
and TEM photos) showing in 
detail how the electrical 
trees circumvent the nanopar- 
ticles and how they propagate 
through the polymer matrix. 
Such photographs are very 
difficult to obtain. 

Where all this leaves us? 
How can we understand in a 
unified way interfaces in 
both classical and nanocom- 
posite insulating systems? 
Can the experimental results 
and simulations on electrical 
treeing with typical machine 
insulation ([52]) provide a 
hint also for a possible ex- 
planation in the nanoscopic 
world? A recent paper on in- 
terfaces posed some pertinent 
questions regarding classical 
and nanocomposite insulating 
materials [89] . From this pa- 
per it is obvious that al- 
though more questions are in 
need of an answer for the na- 
nocomposites, the subject of 



interfaces in classical insu- 
lating systems is by no means 
finished. For example, charg- 
ing of larger interfaces, 
such as found in cable joints 
and terminations, needs to be 
further explored regarding 
the mechanisms of space char- 
ges. Another aspect in need 
of further discussion is whe- 
ther electrical trees may em- 
anate from enclosed cavities 
in conventional polymers. The 
latter question has been par- 
tly answered in [28], [90] - 
where some experimental evi- 
dence was offered as to the 
possibility of electrical 
trees stemming from enclosed 
cavities - but further expe- 
rimental data is needed. 

One thing that should not 
be forgotten - and it is com- 
mon to conventional as well 
as nanocomposite polymers - 
is that an insulating system 
to a great extent is as good 
as its interfaces. This means 
that, no matter whether we 
deal with conventional insu- 
lating materials or nanocom- 
posites, preparation and con- 
struction in both cases has 
to be not only careful but 
meticulous . 

A last remark on the li- 
terature presented here: the 
interested reader may find 
that the authors dwell per- 
haps too much in the older 
scientific literature. This 
is not done because they tend 
to ignore the more recent re- 
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search: they simply would 
like to show that even in the 
old days, the problems were 
more or less the same. More- 
over, they would also like to 
show that fundamental ideas - 
which are with us even today 
- on the various mechanisms 
in insulating materials came 
about quite early. 

Conclusion 

This review - by no means 
exhaustive - tackled the sub- 
ject of the importance of in- 
terfaces both in conventional 
insulating materials and in 
nanocomposite polymers. Whe- 
reas interfaces are to be 
avoided in conventional mate- 
rials, they seem to be a 
blessing in nanocomposites. 



Whereas in conventional mate- 
rials they cause problems of 
compatibility and sometimes 
high field intensifications 
with all the bad consequences 
such intensifications entail 
(i.e. PD, trees), in nanocom- 
posites they seem - up to a 
certain loading - to be de- 
sirable and they prevent (or 
the delay) tree growth. Whe- 
reas in conventional insulat- 
ing systems the introduction 
of more interfaces seems to 
cause sometimes insurmount- 
able problems, the introduc- 
tion of interfaces (because 
of the nanoparticles' intro- 
duction to the polymer ma- 
trix) seems to alleviate elec- 
trical trees and to distrib- 
ute more evenly electric 
fields and space charges. 
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Abstract 

The Beltrami fields are extremely important for the desc- 
ription of time-harmonic electromagnetic fields in chiral 
media. This paper introduces the spherical Beltrami fields 
in chiral media, extending the corresponding known results 
for plane electromagnetic waves. Two key theorems are given 
for the scattering of spherical electromagnetic waves in 
chiral environment from a perfect conductor: the Reciprocity 
and the General Scattering theorems. 
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Introduction 

The name "Beltrami field" 
comes from the fact that 
these fields A satisfy the 
equation type V x A = crA , 
where a a 0 is a reciprocal 
length [1] . The word chiral 
regards the concept of "chi- 
rality", derived from the 
Greek word hand. In geometry, 
a figure is chiral if it can- 
not be mapped to its mirror 
image by rotations and tran- 
slations alone. Typical exam- 
ples are human hands, snail 
shells, spirals and spirals 
in general. 

The discoveries of natural 
optical activity in special 



materials have been known 
since the beginning of last 
century. Though optical acti- 
vity has been considered in 
optics and in quantum mecha- 
nics for many years, its 
analysis within the framework 
of classical electromagnetic 
field theory arose much la- 
ter. Recently, there has been 
a considerable interest in 
the study of scattering and 
diffraction by chiral medium. 
Chiral media are isotropic 
birefringent substances that 
respond to either electric or 
magnetic excitation with both 
electric and magnetic polari- 
zations. Such media have been 
known since the end of the 
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nineteenth century (e.g. the 
study of chirality by Pas- 
teur) and find a wide range 
of applications from many 
sciences [1] . 

In general, the electro- 
magnetic fields inside the 
chiral medium are governed by 
Maxwell equations together 
with Drude-Born-Fedorov equa- 
tions in which the electric 
and magnetic fields are coup- 
led [1] . The chiral medium is 
characterized by the electric 
permittivity e, magnetic per- 
meability p and chirality 

measure (3 . 

Scattering problems with 
incident waves have been stu- 
died: for spherical acoustic 
waves in the works [2] -[4], 
as well in the works [2], [5] 
for spherical electromagne- 
tic waves in non-chiral elec- 
tromagnetic waves and finally 
for spherical electromagne- 
tic waves in chiral media in 
the works [6] , [7] . 

In this work, we derive 
reciprocity relation and gen- 
eral scattering theorem, from 
spherical electromagnetic wa- 
ves emanating from point sour- 
ces in chiral media and scat- 
tered by a perfect conductor. 
Similar theorems are proved, 
analogously, for a chiral di- 
electric [6] , [7] . 

As it is well known [1], in 
a homogeneous isotropic chi- 
ral medium the electromagne- 
tic field is composed of left 
-circularly polarized (LCP) 



and right -circularly polari- 
zed ( RCP ) components which 
are propagated independently 
and with different phase 
speeds. When either a LCP or 
a RCP or a linear combination 
of LCP and RCP electromagnet - 
tic waves is incident upon a 
chiral scatterer then the 
scattered field is composed 
of both LCP and RCP compo- 
nents. So, using Bohren de- 
composition [1], [6], [7], 
the electromagnetic waves are 
expressed in terms of LCP and 
RCP Beltrami fields. The use 
of the Beltrami fields for 
problems of scattering in 
chiral material leads to a 
more simplified relationship. 

This is due to the fact 
that the differential equa- 
tions which satisfy Beltrami 
fields is of first order whi- 
le the electric and magnetic 
field satisfy the Helmholtz 
modified equation which is of 
second-order. Moreover, the 
Beltrami equations and the 
conditions of radiation for 
the Beltrami fields are in 
effect separately for each 
one of them but on the scat- 
terer surface, when applying 
the boundary conditions, both 
Beltrami fields are present. 
This fact makes easier to use 
these equations mainly in 
scattering problems where the 
behavior of the scattered 
field away from the scatterer 
is studied. 

In the second section, we 
formulate the problem for 
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spherical electromagnetic wa- 
ve in a chiral medium. In the 
third section we formulate 
the problem of electromagnet - 
tic wave scattering from a 
perfect conductor as a fun- 
ction of the Beltrami fields 
using the Bohren transforma- 
tion. The spherical Beltrami 
fields are defined in the 
fourth section. Finally, Re- 
ciprocity relation and the 
General Scattering theorem are 
given in fifth section. 

Problem Formulation 

Let Q be a bounded and 
closed subset of M 3 having a 
C 2 - boundary S, i.e. dQ=S. 

The set £T will be referred 
to as the scatterer. The ex- 
terior Q + = R 3 /Q~ of the scat- 
terer is an infinite isotro- 
pic homogeneous chiral medium 
with electric permittivity e, 
magnetic permeability p and 
chirality measure (3. 

The scatterer is filled 
with a isotropic chiral me- 
dium with corresponding phy- 
sical parameters £“, p and 
. All the physical parame- 
ters are assumed to be real 
positive constants. 

Let (E r lnc , H* nc ) be a time- 
' '0 ' '0 ' 

Harmonic spherical electro- 
magnetic wave, due a point 
source located at a point 
with position vector r 0 with 



respect to an origin 0 in the 
vicinity of Q . 

This wave is incident upon 
the scatterer Q and let 

(E^, Hj: c ) be the corresponding 

scattered field. Then the to- 
tal electromagnetic field 
(E, 0 , K e ) in CT is given by 

E^(r) = Ef(r) + E-(r) 
f (!) 
\\ r (r) = H“ c (r) + H* c (r) 

r 0 \/ r 0 \/ r 0 \/ 

We assume that the scat- 
tered field satisfies the 

Silver-Muller radiation con- 
dition 

iE-(r) + rxH-(r) = o(r- 1 ), ^ 

r ^ oo 

uniformly in all directions, 
reS 2 , where S 2 is the unit 

sphere in K 3 , r = |r| , r = r/r 
and rj = (p/e) 1/2 is the intrinsic 
impedance of the chiral me- 
dium in Q + . It is known [8], 
that (2) can be replaced by 

rxE-(r)- n H; c (r) = o(r- 1 ), 

r ^ oo 

uniformly in all directions, 
r e S 2 . The total electric 

field E tot satisfies the boun- 
dary condition: 

n x E^r) = 0, r e S (4) 
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where n is the unit vector 
perpendicular to the outer 
surface S . In view of the 
Drude-Born-Fedorov constitu- 
tive relations [1], the total 
exterior electromagnetic field 
satisfies in the source-free 

region Q + the modified Max- 
well equations 

VxE‘Jr) = P/ 2 E‘ o (r) + 

v (5) 

+ i«p(|-) 2 H[ 0 (r) 

VxH‘(r) = PKX(r) + 

y, t ( 6 ) 

-iws(^) 2 E^(r) 

where 

/c 2 = w 2 ep , y 2 = i _^ 2p2 (7) 



independently and with dif- 
ferent phase speeds. So, we 
consider the Bohren [9], de- 
composition of E r0 and H r 
into suitable Beltrami fields 
Q, r and Q„ r , as follows 

E r 0 = Ql,t 0 + Q/?,r 0 ( 9 ) 

(10) 

where q = (p/8) 1/2 is the impe- 
dance of the material. 

Because the transformation 
Bohren the relation (8) is 
written : 



' Vx Q t , r ; 




Yl 


0 


1 

lO 

1 


_ V X Q/?,r e _ 




0 


~Yr_ 





with |(3/c| <1, [1] . 

We note that (in contrast 
to the non-chiral case) k is 
not a wave number, but just a 
parameter without physical 
significance which has dimen- 
sions of inverse length and 
a) > 0 is the angular fre- 
quency. The system of rela- 
tions (5), (6) can be written 
as 



or equivalent 

j V X Q t ,r. - Y L %r. 

j^v x Q Rrg = — y R Q R/re 



( 12 ) 



where y L , y R are wave numbers 
for Beltrami fields Q, r , Q„ r 

^■L,r 0 r ^K,r 0 

and are given by, 



Yl = 



k 

1 - Ar(3 



Yr = 



k 

1 + /c(3 



(13) 



~V x E* 
r 0 

VxH' 
r 0 

Beltrami Fields 

As it is well known [1], 
[10], in chiral media LCP and 
RCP waves can both propagate 



respectively. Applying the 

transformation Bohren we 

split the field E[ , H[ into 
the Beltrami fields Q L , Q„ , 

so that the problem (2)-(6) 
becomes equivalent to the 

following problem: We are 



(Y* 



k 2 p iwp 
-iwe /c 2 p 



( 8 ) 
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looking for Q L , Q R in Q + 
such that 

V x Q L (r) = y L Q L (r) , r e Q + (14) 
VxQ R (r) = -y R Q R (r), reQ + (15) 

nxQ t (r) = i^|nxQ s (r), reS (16) 

r x Q t (r) + iQ t (r) = o(i) , r -> oo ( 17 ) 

r x Q R (r)-iQ R (r) = o(i), r ^ oo (18) 



Spherical Beltrami Fields 

So, we consider the Bohren 
decomposition of E^ nc and H“ c 
into suitable incident sphe- 
rical Beltrami fields 

and which have the form 

[6], [7], 

Q^(r|pJ = ^B t (r,r 0 ).p t (19) 

Q^(r|p R ) = ^ R B R (r, r B )-p R ( 20 ) 

where 



B L (r , r e ) = + iw + V x i)h{y L \r - r 0 |) (21) 

B„(r, r 0 ) = ^%y R I + fw - V x I)h(y R |r - r 0 |) (22) 

oTT/ Y R 



and h(x) = e lx /(ix) is the ze- 
roth order spherical Hankel 
function of the first kind 

and I = xx + ijjijj + zz is the 
identity dyadic. We recall 
that B L and B R are the funda- 
mental Green dyadics, for the 
Beltrami fields [1], [11]. 

The constant unit vectors p L 
and p R are assumed to satisfy 
the relations: 

•V P, = f 0 P R = 0 

r 0 xp L = ip L (23) 

r.xp B = - ip R 



The constants A L and A R 
are evaluated so that as the 
location of the point source 
goes to infinity along the 
ray in the direction r 0 , the 
spherical fields degenerate 
into plane LCP and RCP Bel- 
trami fields propagated in a 
direction from r a towards 0 , 
with polarizations p L and 
p R , respectively. Using the 
asymptotic forms: 

|r-r 0 | = r 0 -r 0 -r + 0(i), 

1 0 

p— = - f a + O(-p), (24) 

\r ~ r 0 | r e 

r 0 -> oo 
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in the relations (19) 

(20) we obtain [6], [7], 

Q^(r|p j ) = A J .2e_^K j( -r 0 )- 
'0 

.p. e -^ r + 0(i), (25) 

*i n 



Qf; r :(ripj-r,(^^)p t (3D 

Q^(r|P R ) = ^(£^)-P R (32) 

where 



r e -> w 



for j = L, R , where 

K,(-r o ) = J(I-r 0 r 0 -ir o xl) (26) 

K R (-r 0 ) = |(I-r 0 r 0 +ir 0 xI) (27) 

In view of (23) the dy- 
adics K l and K R satisfy the 
relations 



jp L = K,(-r 0 )-p t 
lp R =K R (-r 0 )-p R 

Hence, if we take 



(28) 



r L =2^(y,I + i vv + Vxl) (33) 

f ' = 4 (/Bi+ ^ VV_Vxi) (34) 



and u = |r-r 0 |. Using asympto- 
tic forms (24) for r oo we 
obtain [6], [7], 



Qj^(r | Pj) = F^(f | PjMYjr) + 

+ o(i) 

where 



(35) 



F^(r I Pj) = 



e~ lyjr ' re 

^(y/ s ) 



K,(r)Pj (36) 



A 

A 



4ny 2 

ky L 

,-iw. inz! 

ky R 



then 



(29) 



lim Q^(r | Pj ) = e 



_ _ 



= Qf (r; -r B , pj) 



(30) 



It is convenient to write 
the incident spherical Bel- 
trami fields as [6], [7], 



are the far-field patterns of 
the point source incident 
Beltrami fields, which sati- 
sfy the relations 

r-F/™(r | Pj) = 0 ( 37 ) 

The scattered electric 
field E* c will be dependent 
on the polarization p L , p R 
and will be have the decom- 
position 



E r s C ( r I Pi./ Pr) = Q4r e ( r I P L , P r ) 
+ Q"r e ( r \PuPr) 



(38) 
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whprp *ZS r I P,' Pr) = h (Yf) + 

(39 ) 

Q"r 8 ( r I P L ' Pr) and *ZS r I P Lf Pr) + O(-i-) 

are the corresponding scat- r 

tered Beltrami fields which The functions g ( f) given by 

have the following behavior, J,e 

when r -» oo [l], [10], [12] 

g ,,r 8 ( r) = ^-K t (f)-J n x (y L V x E“(r') + y 2 E~(r') )e^ f r cfs(r') (40) 

g R , ro (r) = K R (r)- f n x (y R V x E~(r') - y 2 E~(r') )e ^ r 'ds( r') (41) 

are the LCP and RCP far-field where S is the surface of 
patterns, respectively [10], the scatterer £T and n is 

and they are dependent also the outward unit vector on 

on the polarization p L and s . We , also, consider two 

p R . We note that the far- locations for the point sour- 

field patterns g 3r (f) satisfy ae ' a and b from which the 

J,re time-harmonic incident sphe- 



the relations 

fg J , re (r) = 0 
■ f xg L re (f) = -ig L ro (f) 
f xg R re (r) = ig Rir0 (r) 



rical electric waves 



Ef c (r|p„p R ) = 

= Ql”o ( r I Pj + Qr"o ( r I P R ) 



for a = a, b , emanate. E* c and 
the corresponding scattered 
electric waves E* c of form 



Reciprocity and General 

Scattering Theorems electric waves E a of form 

For two vector functions < 38 )' have the Bohren decom - 
u and v we introduce the bi- Position of (9) in terms of 
linear form r 10 1 LCP and RCP Beltrami fields. 

In particular, we have 

{U,V> S = r . 

I £inc _ Qinc _j_ Q 2nc 

= f n (u x V x v - v x V x u)dS- (43) ° 0 ° ( 45 ) 

J v ' ' ' psc qSC qSC 

S ( c o Xl_,o T XR,o 

- 2(3y 2 j n (u x v)dS and the following properties 

s for the fields Q“ c and Q“ 0 

[6], [7]: 
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{Qf"a(r; pj, Qf b (r I Qu %)}s 



2y 2 

k 



J n (Q“ c (r; pj x Q“(r | q L , q R ) )ds{ r) 



(46) 



{QlT (r; pJ , Q"„(r I q R )} s = 0 



(47) 



{Q^a(r;p R ), Q"b(r | q L , q R )} s = - n-(Q^ a c (r;p R ) x Q*'„(r | q L , q*))ds(r) (48) 



{Q R ; a c (r; P R ),QL C b(r-|q„q R )} s = 0 



(49) 



Now a "reciprocity theo- 
rem" for spherical electro- 
magnetic waves in chiral me- 
dia is formulated as follows 
[6], [7]: 

For any two incident 
spherical electric waves 

Ea” c (r I P L , p R ) and Ef c (r|p,,pJ 
of (44) and for any scatterer 
in a homogeneous isotropic 
chiral medium, we have: 



-i V L a 



Qr b ( a iq^q R )-PL + 

ae ^Q'bCalqL, q R )-p R = 



y R 

be- iY ^ 

Yl 

be- ±Y ° b 

V R 



Q“.( b I q u q R )-q L + 

Q R 'a(b I q„ q R ) q R 



(50) 



Let 

Ef c (r|p L ,p R ) and Ef c (r | p L , pj 
be two incident spherical wa- 
ves of the form (44). We de- 
fine spherical far-field pat- 
tern generators for LCP and 
RCP spherical Beltrami fields 
by 



0j, b ( a I q^ q R ) = 



ae 



i Yi a 



1 Vi 



g J .b(a|qb,q R )-K J (f)e 1 ^ a ds(r) - 

S 2 

-Q-: b (a|q„q R )} (51) 



for j = L, R . 

This terminology and defi- 
nition is appropriate because 
when both the observation 
point and the source goes to 
infinity, G Jib (a | q„ q R ) , reduce 

to the far-field patterns for 
an incident plane electric 
wave propagating in the di- 
rection -b and of polariza- 
tions q L for LCP and q R for 
RCP fields. Far-field pattern 
generators in acoustic and 
achiral electromagnetic scat- 
tering have been defined in 
[13] . When the point sources 
are transferred at infinity, 
the generators G., b (a | q L , q R ) are 

transformed into far-field 
patterns scattering spherical 
electrical wave, namely [6], 
[7], 
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limGj b (a|q L ,q R ) = g ■ b (-a|q L ,q R ) (52) 

Using this notation, the 
"general scattering theorem" 
for spherical electric waves 
in chiral media is formulated 
as follows [6] , [7] : 

Let 

E a” c O' I fV P R ) and E“ c (r|p,,pJ 
be two spherical electric wa- 
ves of (44) incident upon a 



scatterer in a chiral medium. 
Then it is valid 

G L,b( a I q L , q R )p l + G R,b( a l q t . q R ) p« + 

+ G L,a(b I P u P R ) q L + G R,a(b I P u Pr) q R = 

= h ^) F-(f)dS(f) + 

Zn s’ Yl 

+ f ^^) FR S ,b(r)t/S(f)} (53) 

? Yr 
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Self -Standing End-Fed Geometrically Uniform 
Linear Arrays: Analysis, Design, 
Construction, Measurements and FLOSS 

K.Th. Kondylis, N.I. Yannopoulou, P.E. Zimourtopoulos * 
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Antennas Research Group, Austria - Hellas [2, 3] 

Abstract 



The array factor of a both geometrically and electrically 
uniform array is the simple formula for the complex geometric 
progression sum. This fact, although results in the simplest 
of all possible analytical designs, obviously does not in the 
least simplify the complicated practical problem of feeding 
the array elements using multiple driving points. In order to 
begin the examination of uniform linear arrays with a single 
driving point, this paper presents a compact study of the 
end-fed space arrays with application to geometrically uni- 
form, self -standing linear arrays of parallel dipoles. A num- 
ber of test array models were simulated, constructed and 
their radiation pattern was then measured. The experimental 
and computational results were found to be in good agreement. 
The developed software applications are available through the 
Internet as FLOSS Free Libre Open Source Software. 

Keywords 

End-fed, single driving-point, self -standing arrays 



Introduction 

A space array of 1 < k < N 
parallel, arbitrarily shaped, 
dipoles with identical cur- 
rent density distributions, 
has a complex vector radia- 
tion pattern E = AG (by apply- 
ing the "Principle of Radia- 
tion Patterns Multiplication" , 
i.e. by the linear property 
of the volume integral), whe- 
re G is the radiation pattern 
of the first dipole (which 



"generates the array"; the 
"Generator Pattern"), and A 
is the Array Factor, 




i.e. the complex numerical ra- 
diation pattern of N invented 
"isotropic point sources", 
each of current I R and poin- 
ted by the dipole center vec- 
tor R k , with projection R k to 
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the unit direction vector r. 
Linear Arrays have dipole cen- 
ters on a straight line. Ful- 
ly Uniform Linear Arrays are 
both geometrically uniform, 
i.e. the dipoles are equidis- 
tant, and electrically uni- 
form, i.e. the consecutive di- 
pole currents are of equal 
amplitude and constant phase 
difference [1] . 

End -Fed Arrays 

Perhaps, the simplest pra- 
ctical array is the one con- 
structed from a linear two- 
wire transmission line that 
supports the arms of paral- 
lel, linear, symmetrical di- 
poles, vertical to its plane. 
In order to operate the line 
as a balanced one, it is end- 
fed through a coaxial line 
balun from a coaxial conne- 
ctor, which at the same time 
supports the weight of the 
whole, self -standing, array. 
A possible analysis of such 
an array is shown in Fig. 1. 

In this one port linear 
network, the number of cir- 
cuit voltage and current va- 
riables is equal to the sum 
of 2 variables for the input 
port plus 2N variables for 
the dipole ports plus 4(N-1) 
variables for the line seg- 



ment ports, that is a total 
number of 6N-2 variables. The 
number of the linear rela- 
tions between these variables 
is equal to the sum of 3 re- 
lations, for the generator 
dipole 



1 < N 


: V = V = .V 

i ii' 






I = I + . 1 , 

i ii 


(2) 


plus 


2 relations, for 


the 


last 


dipole 




2 < N 


: V = V , 

O N-l N' 






I„ +1=0 

O N-l N 


(3) 


plus 


3 ( N-2 ) relations, 


for 



the intermediate dipoles 



3 < N, 1 < k < N-2 : 
V 



o V k 



k+1 



A+i' 



o^k + ■'‘k+1 + i^k+l ® 



(4) 



plus N relations, for the cou- 
pling between these dipoles 

N 

1 < k < N : V k =Xz kM I, (5) 

M-l 

plus 2 ( N-l ) relations, for 
each one of the lossless 
transmission line segments of 
electrical length (31 and 

characteristic impedance Z 0k 



0 < |31 k t vtt, v = 1, 2, 3, . . . 



ii^k io^k 


_ ^ z 0k 


cospi k 


1 


oi^k oo^kj 


sinpi k 


1 


cospi k 



pl k = vn, v = 2p + 1, p = 0,l,2, . . . : 0 V k = - i V k , 



Ik' 



I, - .1, = 0 

ok l k 



(7) 
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B1 =vtt, v = 2u + 2, u = 0,l,2, . . . : V = V., L + ,1 = 0 (8) 



1 k N 



Tl 

I 






V L 


i l k-l o | 


| i k N-l o 


<r' 







Fig. 1: End-fed linear array 
of linear dipoles 



that is a total of 6N-3 li- 
near relations to solve for 
6N-2 variables, in terms of 
the source voltage V, which 
is considered as a parameter. 
Anyhow, the resulting current 
ratios are clearly indepen- 
dent of V. 

Single Driving -Point 

Self -Standing Linear Array 

The authors' group has li- 
mited available technical re- 
sources for antenna construc- 
tion. This fact restricts the 
experimentation to thin-wire 
self -standing array models of 
a low total weight that is of 
a small total transmission 
ine length and of a small 
number of dipoles. Thus, the 
practical application of the 
analysis was carried-out for 
N = 2, 3 and 4 dipoles only. 
In order to demonstrate the 
procedure in use, the smal- 
lest case of N = 2 dipoles is 



presented in some detail. In 
Fig. 2, the resulting 6N-3 = 9 
linear relations between 6N-2 
= 9 variables + 1 parameter, 
are shown in a compact form, 
for the case of (J1 *vti, v = l, 

2, ... in (6), with each cell 
value to be the coefficient 
of the variable in the first 
row of its column in an im- 
plied summation, while, if 
(31 1 =vn, v = 1, 2, ... then the 

two last rows with the gray 
background have to substituted 
by the rows of Fig. 3 or 4, 
according to the odd or even 
value of v given in (7)-(8). 



Vi 


V 2 


II 


I 2 


dYl 


oVi 
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Fig. 2: The linear system of 
9 relations for arrays with 
N = 2 and pi x t vtt, v = 1, 2, ... 

Obviously, the complexity 
of the expressions increase 
with the number of dipoles, 
from the simple, of 2 dipoles 
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in which the equality of self 
and mutual impedances have 
been taken into account, re- 
sulting from the system of 
Fig. 2, or of the simplest 
I 2 /I 1 = +l or I 2 /I 1 = -l of Fig. 

3 or Fig. 4 respectively, to 
the most complex one for 4 
dipoles, shown in Fig. 5, 
which covered about one and a 
half A4 page. 
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Fig. 3: Last rows replacement 
for v = 2p + 1, p = 0, 1, 2, ... 
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Fig. 4: Last rows replacement 
for v = 2p + 2, p = 0, 1, 2, ... 
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Fig. 5: Equivalent circuit of the 4 dipoles linear array 



Three Visual Fortran ap- 
plications were developed for 
the computation of current 
ratios. The GUI application 
form for N = 4 dipoles, is 
shown in Fig. 6 [2], [3]. In 
this form, the input data are 
the N-l distances between di- 
poles, the dipole radius and 
length, the length, the char- 
acteristic impedance and the 



velocity factor of each tran- 
smission line segment. In ad- 
dition to current ratios, 
each application produces the 
text files needed by the 
[RadPat4W] application of a 
FLOSS mini-Suite of tools, 
which plots in Virtual Real- 
ity the Generator, Array Fac- 
tor, and Array 3D radiation 
pattern, as well as, their 2D 
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main-plane cuts [4]. The de- 
duced formulas for the deter- 
mination of the 1 + 2 + 3 = 6 
current ratios (I k /I 1 ), were 

mechanically verified using 
Mathematica . 

Array Design, Construction 

and Measurements 

A number of dipole array 
designs were carried out us- 
ing the developed applica- 
tions and the [RichWire] si- 
mulation program, which is a 
fully analyzed, corrected and 
redeveloped edition of the 
original Moment Method thin- 
wire computer program [4], 

[5] . Eleven arrays were fi- 
nally constructed and tested. 
The results for the current 
ratios of four selected ar- 
rays are shown in Tab. 1. The 
arrays were designed for op- 
eration at the frequency of 
1.111 GHz. A two-wire trans- 
mission line, of Z 0 = 200 £2 and 

velocity factor vf = 1 was con- 
structed to feed the dipoles 

[6] . This balanced line was 
then connected to an unbal- 
anced 50 £2 type-N/F base con- 
nector through a 4:1 balun 
made from a segment of RG- 
174U coaxial cable (Z 0 = 50 £2, 

vf = 0.66) with total length 
A/2. 

The arrays were constru- 
cted by bare copper wire of 1 
mm (0.0037A) radius and they 
are self -standing using an 
orthogonal piece of a two- 



sided printed board (3 cm x 
4.48 cm), on which the two- 
conductor line was soldered. 
A few Teflon spacers of low 
relative dielectric constant 
(= 2) were fabricated to fix 
the distance between the two 
line wires at 5.5 mm. The 
measurement system consists 
of a 50 £2 Vector Network Ana- 
lyzer external to an anechoic 
chamber [7] . Each array under 
test was azimuthally rotated 
around its three main axes, 
by a 360° built positioner, 
under the developed software 
control of a built hardware 
controller. The stationary an- 
tenna was a UHF standard gain 
antenna [8]. Fig. 7 and Fig. 
8 show the results for the 4 
test arrays, in 3 groups of 
rows, as follows: 1st group: 
The screen captures of the 
produced Virtual Reality ra- 
diation patterns in dB for 
(a) Generator, (b) Array Fac- 
tor, and (c) Dipole Array - 
2nd group: The constructed 

model, and 3rd group: The 2D 
radiation pattern cuts by (d) 
xOy, (e) yOz, and (f) zOx 
main-planes. In Fig. 7, the 
measurements for the array of 
N = 2, dipoles with equidis- 
tance d = 0.325A, were carried 
out with and without balun. 
In Fig. 8, the array of N = 2 
dipoles with equidistance d = 
0.633A, has been designed to 
exhibit the maximum radiation 
pattern direction off the 3 
main-axes . 
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Fig. 6: GUI for the analysis of a linear array 
of N = 4 dipoles 



Tab. 1: The experimental array characteristics and 

current ratios 



Dipoles 


2 


3 


4 


2 


L/A 


0.5185 


0.5 


0.5 


0.44 


d/A 


0.325 


0.325 


0.325 


0.633 


Ratios 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


V 1 ! 


IVM 


0.534 


0.205 


0.286 


0.280 


0.115 


0.226 


0.527 


^I k /Il[°] 


-80.5 


-114.9 


173 


-115.9 


150.4 


68.1 


49.5 
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Conclusion 

Although, the unavoidable 
mechanical supporting ele- 
ments which exist in the ane- 
choic chamber in the neigh- 
borhood of the antenna under 
test may affect the radiation 
pattern measurements, the ob- 
served differences in the 4 
of the 12 patterns between 
the analysis and simulation 
results, on the one hand, and 
the measurements, on the 
other, have to be charged re- 
spectively: (1) On a cone-cut 
of the array radiation pat- 
tern, instead of the expected 
yOz main-plane cut, in Fig. 
7(e), N = 2 and Fig. 8(e), N = 
4, (2) On an inclination of 
the rotation axis relative to 
the expected linear polariza- 
tion measurement plane, in 
Fig. 8(d), N = 2, and (3) On a 



loosed connection during the 
array rotation, in Fig. 7(d), 
N = 3. These conclusions are 
amplified by the careful stu- 
dy of the corresponding Vir- 
tual Reality space radiation 
patterns. Therefore, under the 
given measurement circumstan- 
ces, the experimental and com- 
putational results were found 
to be in good agreement and 
no attempt was made to modify 
any design or repeat any mea- 
surement . 

The results for the single 
driving-point self -standing 
fully Uniform Linear Arrays, 
i.e. those including the elec- 
trical uniformity, along with 
their application to the con- 
strained pattern design will 
be presented in a future pa- 
per (Part 2) . 
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Abstract 

The aim of this paper is the study of various parameters 
affecting the ageing of transformer oil in distribution 
transformers of 20/0.4 kV. Fifty (50) samples of oil were 
taken from such transformers. The transformers function in 
the major Athens area, Greece. Parameters, such as breakdown 
strength, oil color, humidity, interfacial tension and tan5 
were taken into account. Transformer ageing and lifetime are 
strictly related to the rate of ageing of the whole insulat- 
ing system, and mainly of the oil. 

Absence of transpiring system with silica gel has as a 
result the increase of oil humidity. The lengthy use of 
transformers under heavy load, and consequently under high 
temperatures, is a main factor for oil ageing and oxidation. 
Especially in the last few years, the increased loads re- 
quired by the transformers which have to do also with the 
climatic changes, have as a result the additional stressing 
of the oil. Various arcs, resulting from short circuits in 
the network, have as a consequence the production of gases 
and sludge. Such gases and sludge influence in a negative 
way the insulating properties of the oil as well as its rate 
of ageing. The role of silica gel is stressed. 

Keywords 

Transformer oil, distribution transformer, diagnostic me- 
thods, breakdown voltage, dielectric strength 
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Introduction 

Transformer oil is a most 
important component in trans- 
formers [1] . Transformer oil 
ageing has been studied as 
well as the parameters and 
factors influencing its be- 
havior under a variety of 
stresses (electrical, ther- 
mal, etc.) [2] - [5] . Faults in 
power and distribution trans- 
formers are rare (of the or- 
der of 1% - 2% per year), but 
when they occur, they have 
very serious technical and 
economic consequences. They 
may even lead to dangerous 
situations for human life and 
the environment. Main factors 
affecting the acceleration of 
ageing of the insulation of a 
transformer are humidity, high 
temperature, oxidation and the 
acidity of its oil. The role 
of chemical byproducts of the 
transformer insulation is al- 
so very important for its 
ageing. There is no single 
measurement able to give a 
whole picture of the state of 
a transformer. It must be no- 
ted that the variety of the 
diagnostic methods used, does 
not have as an aim the pre- 
diction of the useful remain- 
ing lifetime of the insula- 
tion, but it tries to reveal 
the increasing probability of 
faults and the corresponding 
decrease of the insulation 
reliability. Needless to say, 
that reduced reliability im- 
plies also reduced remaining 



insulation lifetime. 

With this in mind, in the 
context of the present paper, 
diagnostic methods were used 
in order to study the state 
of distribution transformers 
of the major Athens area. The 
transformers investigated we- 
re 20/0.4 kV. The whole work 
was carried with the aid of 
Public Power Company (PPC) of 
Greece, and more specifically 
with the aid of PPC Trans- 
former Division. Oil was ta- 
ken from 50 distribution 
transformers . 

Diagnostic Methods 

Warning signs about the 
state of a transformer are, 
among others, a big increase 
of partial discharges (» 
2500 pC), a visible deterio- 
ration because of foreign me- 
tallic and carbon particles, 
the presence of humidity in 
the solid insulation about 3- 
4 % and the presence of 
sludge . 

Several diagnostic methods 
were used in order to see the 
quality of the transformers 
in question. The characteri- 
zation of the oil color (DIN 
51517 - ASTM D155 ) was per- 
formed through a device 
(chromometer) including stan- 
dard glass disks and two 
glass jars with lid. The con- 
trol of dielectric strength 
was measured by a typical 
Foster test cell, according 
to IEC 156/95 (Fig. 1). 



FUNKTECHNIKPLUS # JOURNAL 



Vl—54 



ISSUE 4 - YEAR 1 




PARAMETERS AFFECTING THE LIFETIME OF TRANSFORMER OIL IN DISTRIBUTION 



The control of humidity in 
the oil was measured by a 
Metrohm - 684 KF Coulometer, 
which consisted of a glass 
container with a stirrer ti- 
tration in which the reagent 
from container storage is ad- 
ded. The device is fully au- 
tomated and once the experi- 
menter gives the settings, it 



measures the moisture content 
of the oil. The measurements 
were performed according to 
IEC 814. The control of in- 
terfacial tension ( ASTM D971 
- 91) was performed via a 
tensimeter, which gives the 
value in dynes per centimeter 
in a direct reading (Fig. 2). 




Fig. 1: Device for the measurement of dielectric strength 
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The device that performed 
measurements of tan5 and of 
resistivity, is the BAUR-DTL 
fully automated device for 
measuring dielectric losses 
of oils. Such a system has a 
fully automated process for 
measuring dielectric loss, 
relative dielectric constant 
and resistivity. The measure- 
ments were performed with a 



system counting tan5 values 
with maximum accuracy from 
0.00001 to 4.0, measured ac- 
cording to IEC 247. 

It is true that no single 
diagnostic method can give 
full information as to the 
state of a transformer. The 
aforementioned methods may 
give a better picture of its 
state . 




Fig. 2: Tensimeter Cenco du Nouy 
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Experimental Results 

The sampling was performed 
with due care, and according 
to general standard practice, 
i.e. during sampling there 
must not be any dust and hu- 
midity in the nearby space, 
the sampling cells must be 
clean and dry and they must 
be washed with oil from the 
transformer which is to be 
checked, the samples must be 
protected from light and they 
must be taken while the oil 
is hot but due time should 
pass before the various par- 
ticles settle down. Sampling 
was done from the bottom of 
the transformers. The glass 
cells used for oil sampling 
were big enough so that they 
contained enough oil volume 
for additional measurements, 
if needed. No measurement was 
performed based on only one 
value or sample. All results 
must be verified with addi- 
tional samples and measure- 
ments. Color measurements of 
the oil samples were based on 
specification DIN 51517 
ASTM D155 (Fig. 3). Dielec- 
tric strength measurements 
were based on specification 
IEC 156/95. Humidity measure- 
ments were based on specifi- 
cation IEC 814. Interfacial 
tension measurements were ba- 
sed on specification ASTM 
D971-91 . Loss factor measure- 
ments were based on specifi- 
cation IEC 247. 

Transformer oil is charac- 



terized as good, if its color 
is not dark, its breakdown 
voltage more than 40 kV, its 
humidity less than 10 ppm, 
its resistivity more than 3 
GD-m, its loss factor 0.1, 
its interfacial tension more 
than 28 mN/m. It is accept- 
able if its breakdown voltage 
is between 30 to 40 kV, its 
humidity between 10 and 25 
ppm, its resistivity between 
0.2 and 3 GQ-m, its loss fac- 
tor between 0.1 and 0.5, its 
interfacial tension between 
22 and 28 mN/m. In the case 
of an acceptable oil, sam- 
plings should be carried out 
more frequently, and the con- 
tent in foreign particles and 
contained water should be 
controlled. If the breakdown 
voltage and humidity are near 
the limit values, then the 
oil should be filtered and 
cleaned. An oil is poor when 
its color is dark (the dark 
color is an indication of 
pollution or ageing), its 
breakdown voltage less than 
30 kV (in this case the oil 
should be replaced or thor- 
oughly cleaned), its humidity 
more than 25 ppm (in such a 
case the oil should be re- 
placed or thoroughly clea- 
ned), its resistivity less 
than 0.2 GQ-m, its loss fac- 
tor more than 0.5, its inter- 
facial tension less than 22 
mN/m (in this case a check 
should be carried out for the 
presence of sludge). 
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Generally it is suggested 
that in case one of the pa- 
rameters mentioned (color, di- 
electric strength, humidity, 
interfacial tension, loss fa- 
ctor, resistivity) is accept- 
able, a more frequent sam- 
pling should be done in order 
to better monitor the oil 
quality. In the case any of 
the mentioned parameters is 
poor, cleaning or replacement 
of the oil is suggested. In 
any case, it is good to con- 
trol all parameters concerned 
in order to have a clear pic- 
ture of the state of the oil. 

From the 50 oil samples 
investigated (taken from 50 
distribution transformers), 3 
of them showed breakdown vol- 
tage lower than 30 kV (i.e. 
percentage of 6%), 45 samples 
indicated humidity higher 



than 25 ppm (i.e. percentage 
90%), 28 samples presented 
interfacial tension lower 
than 15 dynes/cm (percen- 
tage 56%), whereas 44 samples 
showed a color higher than 1 
1/2 (percentage 88%) (Fig. 
4) . It was evident from the 
data - mentioned in detail in 
[6] - that the more the number 
of years in service, the more 
the number of samples with 
lower dielectric strength. 
This is due to the fact that 
with the years there is also 
an increase of humidity in 
oil as well as an increase of 
the oxidation byproducts. 
Trans formers which are in 
cities or in urban localities 
do not usually suffer light- 
ning strokes, consequently 
the only stresses come from 
switching overvoltages and/or 
some high currents. 




Fig. 3: Chromometer Hellige Comparator 
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Fig. 4: Typical glass colored disc 



As the years pass by, hu- 
midity increases because fluc- 
tuations in temperature cause 
an increase of humidity en- 
tering the main oil volume. 
In distribution transformers 
having silica gel, the latter 
looses its absorbing property 
with an increasing number of 
years in service. Consequently 
silica gel should be replaced 
at regular intervals. The in- 
creasing number of years in 
service influences also the 
oil oxidation. Oxidation by- 
products in turn affect tan5. 
The loss factor, however, is 
not influenced from the exis- 
tence of humidity in the sam- 
ples and for this reason is 
not a criterion for its exis- 



tence. Interfacial tension is 
reduced with the years in 
service. This is due to the 
increasing quantity of humid- 
ity as well as to the bypro- 
ducts of oxidation. There is, 
however, a number of samples 
with rather reduced interfa- 
cial tension although the 
number of years in service 
was not that big. This may be 
due to the quality of the oil 
used or to the fact that it 
contained only a small quan- 
tity of anti-oxidants. Fi- 
nally, the oil color changes 
with the years in service. A 
change of color may also be 
due to the overcharging of a 
transformer . 
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Discussion 

In this work, sampling from 
50 transformers of 20/0.4 kV 
was carried out. Sampling and 
measurements were performed 
as suggested by the various 
international specifications 
[7] . All transformers were in 
the major Athens area. It 
should be noted that the trans- 
formers were not all from the 
same manufacturer. There were 
transformers from various ma- 
nufacturers, both from Gree- 
ce and abroad. This is one of 
the complications of the Greek 
Electricity System and cer- 
tainly, this fact hinders any 
thought for a statistical ap- 
proach. In case of poor oil 
color, this may be a symptom 
of humidity. Poor oil dielec- 
tric strength may imply that 
the oil should be cleaned or 
replaced. If the humidity le- 
vel is high, the oil has to 
be again either cleaned or 
replaced. If the interfacial 
tension of oil is acceptable, 
this means that a control for 
the possible presence of by- 
products may be performed. 
Acceptable values of oil re- 
sistivity may be tolerated. 
In most cases, tan5 was good. 
The parameters which changed 
most were the color, the oil 
dielectric strength, the lev- 
els of humidity and the oil 
interfacial tension. From the 
above, it is evident that one 
cannot pronounce any verdict 



on the quality of a trans- 
former oil based on only one 
parameter . 

It is evident from the re- 
sults that there was absence 
of a transpiring system with 
silica gel in most investi- 
gated transformers. This ab- 
sence has as a result that 
the humidity of the atmos- 
phere is in contact with the 
oil. The functioning of a 
transformer under heavy load, 
and consequently its func- 
tioning at high temperatures, 
results in an acceleration of 
its ageing. High temperatures 
cause oxidation of the oil. 
In addition to that, in the 
last years, increased loads - 
as a consequence of climatic 
changes (such as increased 
temperature and humidity) - 
are required and they result 
in the stressing of the oil. 
Arcs, because of short cir- 
cuits, create gaseous byprod- 
ucts and sludge. The latter 
influence in a negative way 
the dielectric strength of 
the oil. Another point which 
should be considered is that, 
with the continuing effort 
from the side of the manufac- 
turers for the reduction in 
the size of the transformers, 
the oil tends to be more eas- 
ily thermally stressed. This 
results in an accelerated 
rate of rise in temperature 
and an accelerated rate in 
its ageing. 
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Generally speaking, the 
frequent oil sampling helps 
in pointing out potential 
trouble spots in the network. 
It goes without say that, as 
in this paper, a variety of 
measurements may give a bet- 
ter picture of the actual 
state of the transformers ra- 
ther than isolated measure- 
ments or measurements confi- 
ned only on a single quantity 
(e.g. dielectric strength). 
Ageing is a complex process 
and, as is well known, a va- 
riety of factors play a role. 
Ageing factors contribute not 
only to the ageing of the 
oil, but they may also inter- 
act among themselves, render- 
ing thus the prediction of 
lifetime even more difficult 
[8], [9]. Breakdown of trans- 
former oil is also a complex 
phenomenon determined mostly 
by the operating conditions 
[ 10 ], 

As a way out of the vari- 
ous problems with transformer 
ageing, the construction of 
additional substations is pro- 
posed. This will reduce the 
load of each transformer. Need- 
less to say that, the contin- 
uous checking of the whole 
distribution network and the 
requirement from the manufac- 
turers to observe all the in- 
ternational norms and speci- 
fications, are necessary con- 
ditions for the good func- 
tioning of the whole system. 



Special attention should 
be paid to the way one gets 
the oil samples. Sampling 
should be done with due care 
and care should be taken in 
order not for the sample to 
contain water. The samples 
should be taken from the 
transformers when the oil is 
still warm. Foreign particles 
should be allowed to settle 
at the bottom, so that no un- 
due influence can be recorded 
in the various measurements. 

In the case the results 
are contradictory, one has to 
repeat the sampling and the 
corresponding measurements. 
It is true that the manufac- 
turers, try to reduce the pro- 
duction cost of the trans- 
formers and their respective 
materials by reducing the si- 
ze of the transformers. This 
on one hand has as a result 
the possibly higher break-down 
strength of the oil, on the 
other hand, however, has as a 
result the faster increase of 
the oil temperature, and con- 
sequently the faster ageing 
and stressing of the mechani- 
cal parts of the transformer . 
Anyhow, a better maintenance 
of the distribution transfor- 
mers should be followed as 
well as the placing of more 
transformers of larger power 
in localities with an in- 
creasing load demand. 

The discussion as well as 
the conclusions of the pre- 
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sent paper are in line with 
previously published work 

[11] , [12]. Such work was 
carried out on transformers 
of 20/0.4 kV as well as of 
150/20 kV. A statistical ap- 
proach cannot be offered for 
the time being, as it was 
done before [13], for the 
reasons given above. 

It is to be hoped that the 
next steps in our work will 
be a classification of the 
distribution transformers ac- 
cording to manufacturer, ac- 
cording to their power (kVA) 
and according to the previous 
history (which includes years 
in service, all faults, short 
-circuits, lightning strokes, 
switching operations etc.). 
Thereafter, an effort will be 
made in order to statisti- 
cally analyze the data. 

Concluding, it is fitting 
to note that the present pa- 
per does not claim any origi- 
nality regarding any new di- 
agnostic techniques. This pa- 
per is an application of al- 
ready existing and well known 
specifications. It is a con- 
tinuation of an effort star- 
ted few years ago, as men- 
tioned in [7], [11], [12]. In 
this respect, the present pa- 
per presents new results, al- 
though the methods used are 
the same as in [6], [7], [11], 

[12] . In other words, this 
paper presents the state of 
transformer oil in some dis- 



tribution transformers of a 
given area at a given moment 
in time. Needless to say that 
continuous monitoring is 
necessary in order to check 
and confirm the transformer 
oil quality. Continuous moni- 
toring, however, is often not 
enough. If possible, it should 
go together with fault dete- 
ction, namely partial dischar- 
ge (PD) detection. PD dete- 
ction can help preventing a 
defect from developing into a 
fault, e.g. a short circuit 
[14], [15], [16]. Nevertheless, 
condition monitoring as pre- 
sented here, should be done 
on a regular basis. 

Conclusion 

Transformers of 20/0.4 kV 
from the major Athens area 
were investigated. Various pa- 
rameters (oil color, dielec- 
tric strength, humidity, in- 
terfacial tension, tan5 and 
resistivity) of the oil were 
studied with standard methods 
based on international speci- 
fications. No verdict on trans- 
former oil quality can be ba- 
sed on only one parameter. A 
multitude of parameters is 
needed in order to pronounce 
a correct verdict on oil qua- 
lity. The continuous control 
and monitoring of the distri- 
bution transformers is neces- 
sary in order to avoid prob- 
lematic situations. The role 
of silica gel is emphasized. 
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Silica gel should be used in 
all transformers, so that hu- 
midity can be absorbed. The 
quantity of humidity increases 



with the number of years in 
service, consequently, distri- 
bution transformers should be 
checked at regular intervals. 
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Abstract 



Electrical machine insulation consists basically of epoxy 
resin and mica foils. Both are good insulating materials and 
they proved to withstand large partial discharges and/or the 
combined attack from electrical, thermal and mechanical 
stresses. Nowadays alternatives to the traditional insula- 
ting materials exist, namely those of nanocomposite poly- 
mers, which present somehow improved performance regarding 
the aforementioned stresses. It is the aim of the present 
paper to investigate mechanisms of electrical treeing and/or 
breakdown in machine insulation as well as to study possible 
improvements with the aid of nanocomposite polymers. 



Keywords 

Machine insulation, nanocomposite polymers, breakdown, par- 
tial discharges, breakdown strength, pre-breakdown phenomena 



Introduction 

Machine insulation is ge- 
nerally characterized by hard 
materials, capable to with- 
stand partial discharges of 
quite large magnitudes. An 
electrical insulation in high 
voltage machines must have a 
high breakdown strength, a 
good long term functioning 
without possible problems em- 



anating from degradation ef- 
fects (partial discharges 
(PD), treeing phenomena etc.) 
and certainly rather small - 
if any - leakage currents. 
Such insulation must also ha- 
ve a satisfactory thermal per- 
formance and must withstand 
relatively high temperatures. 
The electrical and thermal 
properties of such insulation 
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must not be deteriorated be- 
cause of some extreme condi- 
tions, i.e. the highest pos- 
sible temperature to be tol- 
erated should not provoke any 
alterations to the insulation 
itself. The electrical machi- 
ne insulation should have good 
long term functioning and it 
should present only small 
changes (to be tolerated) af- 
ter many thermal cycles. Mo- 
reover, electrical machine in- 
sulation should have a satis- 
factory mechanical behavior, 
i.e. high mechanical strength 
which should stay high even 
at higher temperatures. 

Thermal stressing may cau- 
se some sort of breakdown on 
the surface of mica sheets, 
whereas mechanical stressing 
may cause also fissures in 
epoxy resin and in mica. Elec- 
trical stressing may cause 
partial discharging in possi- 
ble existing imperfections of 
the insulation and thus lead 
to electrical failure [1] . 
Needless to say that all pro- 
perties - electrical, mechan- 
ical and thermal - are vital 
for the reliable operation of 
rotating machines. The insu- 
lation should have satisfac- 
tory electrical properties, 
it should withstand the ex- 
pansion and contraction dur- 
ing temperature cycles and it 
should respond in a good way 
to mechanical stresses. 

Rotating machines can be 
divided into two categories: 



those with voltage ratings 
less than 6.6 kV and those 
with voltage ratings higher 
than 6.6 kV. Mica has been 
used for years in the elec- 
trical machine industry. Mica 
sheets with a backing of 
glass cloth and of binding 
material, like epoxy resin, 
have been the classical insu- 
lation systems for rotating 
machines. Generally speaking, 
such insulating systems have 
proved to be reliable [2]- 

[4] ■ 

Research on possible break- 
down mechanisms in such sys- 
tems revealed that PD and 
electrical treeing may lead 
to breakdown [5]. As is poin- 
ted out in published work, 
electrical trees propagate 
through the weaker material 
and they tend to reach the 
opposite electrode, resulting 
thus in total failure [2]- 

[5] , 

It is the aim of the pre- 
sent paper to explore mecha- 
nisms of failure both in con- 
ventional insulating materi- 
als used for rotating machi- 
nes as well as to investigate 
alternative insulating sys- 
tems based on a new series of 
insulating materials, namely, 
the nanocomposite polymers. 
Simulation data will be pre- 
sented for conventional insu- 
lating materials and possible 
improvements will be sugge- 
sted. It has to be noted that 
this review is by no means an 
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exhaustive one. Rotating ma- 
chine insulation is a complex 
insulating system and, in the 
context of the present paper, 
only some aspects of it will 
be treated. It should also be 
added that the present review 
does not deal with the sub- 
ject of enameled wiring, 
which will possibly consist 
the topic of another review 
paper . 

At this stage, it would be 
fitting to state that review 
papers are useful even today. 
Why? Because they give some 
insight to the newcomer as 
well as to the experienced 
scientist regarding a partic- 
ular subject. Although we are 
flooded with much information 
from the Internet with prac- 
tically thousands and thou- 
sands of references (books, 
papers, journals etc.), a re- 
view paper can always be a 
starting point for the inter- 
ested person. The purpose of 
a review paper is not simply 
to collect information on a 
particular subject, it is ra- 
ther to offer - besides a 
wealth of information and the 
relevant references - com- 
ments and interrelations be- 
tween experimental data, cri- 
ticism and approaches for pos- 
sible new theoretical models, 
it is probably the right pla- 
ce to also propose some new 
insights from the reviewer, 
it is the place to offer some 
thoughts for possibly new de- 



velopments. This is why we be- 
lieve in the validity of such 
efforts, this is why we think 
that reviews will always be 
useful, despite the flood of 
new scientific and technical 
information, almost on an eve- 
ry second-basis. 

Conventional Insulating 

Materials for Rotating 

Machines 

Mica sheets are used in 
rotating machine insulation 
as traditional insulating ma- 
terial. Mica is a natural mi- 
neral. Its crystalline nature 
gives very strong bonds in 
one plane and very weak Van- 
der-Waals's forces in the pla- 
ne normal to this. The conse- 
quence of that is that mica 
can be split easily into fla- 
kes [2] . Mica has excellent 
tracking strength, high break- 
down field strength, very good 
resistance to PD, high volume 
resistance as well as good 
thermal stability up to 6000 
C [2]. For high voltage ap- 
plications, mica sheets use 
as bonding material epoxy re- 
sin, a thermosetting material 
with very good electrical 
properties and good resis- 
tance to PD. Moreover, for a 
resin to be suitable for long- 
term operation, it requires 
high thermal stability with 
low electrical loss at ser- 
vice temperature and at power 
frequency, excellent adhesion 
to mica, high resistance to 
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moisture, chemicals and other 
contaminants, high mechanical 
strength over a range of ser- 
vice temperatures, dimensio- 
nal stability, ability to ope- 
rate at higher temperatures 
and a short cure time at 1500 
C to 1600 C [6] . In typical 
rotating machine insulations, 
mica sheets form a sort of 
sandwich with epoxy resin, 
thus rendering the electrical 
breakdown of such a combined 
system rather difficult [7] . 

One of the problems facing 
the mica/epoxy resin insula- 
tion is the one of electrical 
trees, which may grow and 
eventually bridge the gap be- 
tween the electrodes causing 
thus ultimate failure. Elon- 
gation of electrical trees 
has been experimentally ob- 
served in [8], [9]. In such a 
case, the mica sheets consist 
the harder material and the 
electrical trees propagate 
via the epoxy resin, which is 
the weaker material. Electri- 
cal trees propagate, gener- 
ally speaking, more easily in 
epoxy resin, a fact also con- 
firmed in another paper [10] . 

Normally prior to electri- 
cal treeing, PD take place in 
defects in the machine insu- 
lation. Such defects may come 
about from construction or 
from the stressing of the in- 
sulation. As said, the stres- 
sing in machine insulation 
can be multi-factor stress- 
ing, i.e. the insulation may 



be stressed because of high 
voltage, thermal cycles as 
well as from mechanical load- 
ing [11] . Defects can come 
about as enclosed cavities, 
delaminations, problematic 
interfaces, possible enclosed 
foreign particles etc. The 
consequence of all these is 
local electric field enhance- 
ment, PD activity which sub- 
sequently may result to insu- 
lation damage. Such PD can be 
quite intense in the order of 
1-10 nC [12], [13]. These phe- 
nomena can have a cumulative 
effect and cause aging and 
shortening of the lifetime of 
machine insulation [14], [15]. 

Rotating machine insula- 
tion systems suffer from what 
most of composite insulating 
systems suffer, i.e. the pre- 
sence of interfaces. Mica 
sheets and epoxy resin con- 
sist of a system with multi- 
ple interfaces. Interfaces 
may encourage electric field 
intensifications in the wea- 
ker material and, thus, the 
cause of deteriorating phe- 
nomena. On the other hand, 
thinner (and consequently mo- 
re) mica sheets may delay the 
discharge process in that 
discharges lose energy at the 
interfaces, i.e. a discharge 
having penetrated one layer 
could not enter the next 
layer of material until the 
spot on the interface, cen- 
tred on the channel, had been 
charged to a potential which 
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could produce a field compa- 
rable with that of the chan- 
nel at the level in question 
[16]- 

Breakdown mechanisms in- 
side the insulation may start 
from defects due to excessive 
electrical field. Previous 
work done on epoxy resin sam- 
ples showed that it is possi- 
ble that trees emanate from 
enclosed cavities, causing 
thus conditions for further 
propagation and eventual fai- 
lure [17] . Although such ema- 
nating trees are still put 
into question from some re- 
searchers [18], experimental 
evidence can hardly be re- 
futed. PD cause pits on the 
inner surface of such cavi- 
ties possible and then elec- 
trical trees may ensue. One 
aspect that should be stres- 
sed is that of the applied 
voltage: for meaningful com- 
parisons of data (perhaps with 
several years interval) and/ 
or comparison of data in the 
same laboratory or factory, 
it is essential to use iden- 
tical wave shapes of volta- 
ges [19] . It is something 
that people tend to forget 
but something that comes out 
when PD measurements - com- 
parisons of such measurements 
at different times - have to 
be performed. 

Mica barriers delay tree 
propagation. Depending on the 
dielectric constants of mica 
sheets and of epoxy resin as 



well as on the threshold vol- 
tages, electrical trees may 
take different forms but, in 
general, they seem not to pe- 
netrate the mica barriers 
[20]. Mica barriers may re- 
sult in a major increase in 
breakdown time, this increase 
being depended on both the 
tree growth time and the set- 
in of the failure time [21] . 

The interplay and interde- 
pendence between PD and elec- 
trical treeing has been shown 
before, where in narrow holes 
of short length small PD in a 
rather high number may be pro- 
duced whereas in holes of lar- 
ger diameter and longer length, 
fewer PD but with larger mag- 
nitudes will ensue [22] . In any 
case, trees tend to grow a- 
round the mica barriers [23] . 

A crucial factor determin- 
ing tree propagation along a 
mica barrier is the type of 
chemical bonding between mica 
and epoxy resin. The stronger 
the bonding, the higher the 
resistance to the tree propa- 
gation [24] . Imperfections 
may result from imperfect mi- 
ca sheet overlapping, from 
the creation of cavities in 
parts which are at the edges 
of windings, from not so 
smooth mica sheets or from 
abrupt interruptions of mica 
sheets (because of construc- 
tional faults) [25]. Further- 
more, the layered mica can 
delaminate under thermo-mecha- 
nical stresses and thus cause 
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cavities, which in turn will 
lead to PD [25] . Examples of 
imperfections are given be- 
low. In Fig. 1, a winding in- 
sulation is shown and the 
different radii of mica sheets 
are noted. In Fig. 2, imper- 
fections in winding insula- 
tion are noted. 

Experimental evidence that 
a mica barrier may withstand 
electrical treeing much bet- 
ter than epoxy resin was gi- 
ven in [10], where it was al- 
so noted that the breakdown 
strength of such a combina- 
tion depends on the thickness 
of mica sheets, on the thick- 
ness of epoxy resin, the tem- 
perature, the type of epoxy 
resin as well as the cleanli- 
ness of both materials. 

Generally speaking, the 
time to breakdown is the sum 
of the time from the initial 
PD activity and the creation 
of initial tree channels and 
the growth time of trees to 
the final breakdown. For some 
authors, there is an incuba- 
tion period during which PD 
activity is barely detectable 
and trees grow only slightly, 
then a period of tree expan- 
sion follows and finally a wi- 
dening of the smaller tree 
channels which eventually 
leads to bridging of the elec- 
trodes and the breakdown [9] . 
For others, two stages of tree- 
ing are observed: first, the 
inception period which may be 
for very many cycles, follo- 



wed by a relatively short pe- 
riod of tree growth. PD de- 
tection reveals that a tran- 
sition is accompanied by big 
increase in the PD magnitude, 
which persists until break- 
down [26] . A good account of 
PD measurements and the ef- 
fect of PD in rotating ma- 
chines were given in [27], 
where examples of "good" 
(i.e. relatively free of PD) 
and "bad" (with delamina- 
tions) windings were presen- 
ted. In [27], it was empha- 
sized that for any comparison 
between PD measurements, the 
experimental conditions play 
a predominant role, a state- 
ment that echoes reference 
[19] . Certainly, a condition 
assessment of windings can be 
done by continuously monitor- 
ing the PD activity, taking 
into account that an increase 
in discharge activity occurs 
when the insulation is eroded 
and also that a PD activity 
can manifest itself both as 
internal PD activity and as 
surface discharges [28] . A 
good account of the relation 
between PD and tree structu- 
res was given in [29], where 
it was emphasized that elec- 
tron avalanches, field fluc- 
tuations arising from the 
discharges themselves, local 
variations in permittivity and 
resistance of the insulation 
can play a decisive role for 
the electrical tree propaga- 
tion. Minor variations of 
trapped space charges may 
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lead to preferred directions tions and, thus, may also af- 
for new tree channel forma- feet the tree propagation. 




Fig. 1: Winding insulation. Note the different radii which 
cause the different bending of mica sheets (after [10]) 




Fig. 2: Imperfections of winding insulation at the edges 

(after [10]) 
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It is not to underestimate 
that several years ago, prob- 
lems with very small PD - not 
to say with phenomena below 
the inception voltage - were 
noticed in turbine generator 
armature winding insulations, 
namely that although such win- 
dings withstood standard test- 
ing, they failed short after 
the beginning of their ser- 
vice. Why was that? Possibly 
because extremely small PD, 
not easily detected by con- 
ventional methods, are at 
work and little by little 
could deteriorate the insula- 
tion of such windings [30] . 

Mica sheets offer a good 
protection to breakdown paths 
and/or to electrical treeing. 
This certainly depends on the 
applied voltage. Simulations 
in epoxy resin/mica sheets 
show that there is a voltage 
limit beyond which breakdown 
even of mica sheets is possi- 
ble. From a certain voltage 
value upwards, mica sheets 
are also prone to treeing, as 
the following figures show. 
In Fig. 3, the applied volt- 
age is 28 kV, which creates a 
rather distinguished form of 
treeing in the epoxy resin 
but which cannot penetrate the 
mica sheet. On the contrary, 
in Fig. 4, with another 
higher applied voltage, the 
mica sheet is penetrated by 
the electrical trees and in 
Fig. 5, there is a complete 
failure of the system epoxy 



resin/mica sheet. It is to be 
noted that in Fig. 4 as well 
as in Fig. 5, the trees in 
the epoxy resin are of bush- 
type and in the mica sheet of 
branch-type. This is because 
trees in the weaker material 
are far more numerous and 
thus they are interconnected 
much more densely than in the 
stronger material. The simu- 
lations shown in Figs. 3, 4 
and 5 were performed with the 
method of Cellular Automat 
[31] -[33]. The different ty- 
pes of electrical trees in 
epoxy resin (bush-type) and in 
mica sheet (branch- type) de- 
pend on the applied voltage 
as well as on the type of the 
insulating material [34] . In 
Figs. 4 and 5, the progres- 
sion of treeing towards the 
opposite electrode can well 
be seen. It is indeed a ques- 
tion of time before the tree- 
ing structure reaches the op- 
posite electrode. 

All in all we observe that 
conventional insulating mate- 
rials functioned more or less 
satisfactorily. Both experi- 
ments and simulations indi- 
cated that breakdown paths 
follow the easiest way to the 
other electrode, i.e. through 
the epoxy resin, which is the 
weaker of the two materials. 
Breakdown of the mica sheets 
may be possible but this de- 
pends on their thickness as 
well as on the voltage ap- 
plied . 
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Fig. 3: Propagation of electrical tree in epoxy resin and mica 
sheet. Applied voltage 28 kV, breakdown strength of epoxy re- 
sin is 26 kV/mm, breakdown strength of mica sheet is 35 kV/mm 




Fig. 4: Propagation of electrical tree in epoxy resin and mica 
sheet. Applied voltage 34 kV, breakdown strength of epoxy re- 
sin is 26 kV/mm, breakdown strength of mica sheet is 35 kV/mm 
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Fig. 5: Propagation of electrical tree in epoxy resin and mica 
sheet. Applied voltage 34 kV, breakdown strength of epoxy re- 
sin is 26 kV/mm, breakdown strength of mica sheet is 35 kV/mm 



Is there a way to improve 
even more the electrical per- 
formance of such a composite 
system? The following section 
will concentrate on some of 
the modern materials proposed 
for machine insulation. 

(As an explanatory note to 
this paper: It is evident 
from the present review that 
here we are not interested in 
an estimation of lifetime ex- 
pectancies of the insulation 
in connection with the qual- 
ity of the insulation itself. 
Models of lifetime expectan- 
cies have been developed and 
commented upon elsewhere [35] . 
It should also be emphasized 
that, the present paper is not 
concerned with the technica- 
lities of the measurements of 
PD in rotating machine insu- 



lation either [36] ) . 

Nanocomposite Polymers for 

Rotating Machine Insulation 

A recent paper on nanodi- 
electrics applications has 
pointed out that "three types 
of insulations have been de- 
veloped with great improve- 
ments in the resistance to 
partial discharges: first in 
random-wound wire enamel; se- 
cond in form-wound strand 
enamel; and third in form- 
wound stator bar insulation 
or ground wall insulation" 
[37] . Of these aforementioned 
applications, the one that 
interests us most for the 
time being is the third one, 
i.e. that of stator bar insu- 
lation. The other two types 
of insulation may constitute 
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the subject of yet another 
review paper. In the context 
of the present paper, we will 
concentrate on some off the 
novel insulating materials 
concerning the stator bar in- 
sulation . 

Results reported on 13.8 
kV form-wound hydro generator 
bars made of epoxy resin with 
a high percentage of nanopar- 
ticles of silica, indicated a 
very optimum performance re- 
garding their thermal and me- 
chanical properties. The elec- 
trical performance of such an 
insulation w.r.t. the PD and 
electrical treeing resistance 
was also noted [37], [38]. 
Better resistance to PD and 
electrical treeing implies 
less insulation thickness. 

Before starting an evalua- 
tion of nanocomposite materi- 
als for machine applications, 
let us be clear on a certain 
point: the novel materials 
may offer some new insights 
regarding the endurance of 
the insulation of electrical 
machines, but the condition 
-based maintenance treatment 
is not expected to change 
that dramatically. In other 
words, we are still far from 
new aging models and/or ac- 
celerated laboratory testing 
[40] . In yet other words, sy- 
nergy effects may occur also 
with the new materials and, 
when two or more aging mecha- 
nisms are at work simultane- 
ously, the total degradation 



may not be a single algebraic 
sum of the separate degrada- 
tion outcomes [39], [40]. 

It is rather clear that 
from the two basic components 
that constitute the rotating 
machine insulation, the one 
that is expected to be prone 
for the addition of nanopar- 
ticles, is the epoxy resin. 
Inorganic nitrides and oxides 
(such as, AlN, BN, Si02, ZnO, 
AI 2 O 3 and Ti02 among others) 
may be added in small amounts, 
and if homogeneously disper- 
sed, they may show potential- 
ly better electrical and ther- 
mal properties than the con- 
ventional epoxy resin [41] . 
The need for good dispersion 
of the nanoparticles as well 
as the important role of in- 
terfaces in the nanocomposite 
polymers has been discussed 
elsewhere and there is no 
need for any repetition here 
[42]. Boron Nitride (BN) has 
been used to enhance the di- 
electric properties of the 
groundwall insulation system 
for generators because it has 
a reasonably high resistivity 
( 10 15 ohm. cm) and breakdown 
strength (53 kV/mm) as well 
as a rather small relative 
permittivity (around 4). Such 
inherent properties render BN 
a good nanoparticle material 
to be added to epoxy resin 
[41] . Moreover, the addition 
of BN nanoparticles improved 
the thermal conductivity of 
epoxy resin. In another pub- 
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lication, epoxy resin was mi- 
xed with AI 2 O 3 nanoparticles 
[43] . The breakdown strength 
of such a nanocomposite was 
comparable to the breakdown 
strength of conventional epo- 
xy resin but with a lower 
scatter in values. In [43], 
particular attention was paid 
to the fact that nanoparticle 
agglomerations were dependent 
on the method of preparation 
of the nanocomposite. As is 
well known, nanoparticle agglo- 
merations affect the PD, the 
electrical tree propagation 
and the breakdown strength of 
nanocomposites [44] . Zirconia 
nanoparticles added to epoxy 
resin offer an improvement of 
the breakdown strength as 
well as of the thermal and 
dielectric properties. Addi- 
tion of zirconia nanoparti- 
cles to about 5 wt% gives a 
higher breakdown strength com- 
pared to conventional epoxy 
resin [45]. It must, however, 
be noted that in [45], there 
are in certain parts of the 
paper discrepancies between 
the text and the experimental 
data . 

Fillers in insulation sys- 
tems for rotating machines 
are investigated in [46], a 
most thorough review publica- 
tion. Without going into many 
details about the application 
of nanocomposites for machine 
insulation and the role of 
nanoparticles, that paper 
points out several vital as- 



pects of insulation engineer- 
ing, namely, the effect of 
electrical, mechanical and 
thermal stresses, which leads 
to delaminations and/or void 
formation. Nanoparticles are 
deemed to restrain the charge 
transport processes. As a con- 
sequence, nanocomposites ex- 
hibit less space charge for- 
mation than conventional po- 
lymers or their microparticle 
counterparts. The novel ap- 
proach of [46] is that it in- 
sists on a concept for void 
formation that is based on a 
micro-mechanical approach, i. 
e. void formation can be con- 
sidered as a process by which 
the insulating material ac- 
commodates mechanical energy, 
no matter whether this energy 
is purely mechanical, elec- 
tro-mechanical or thermo-me- 
chanical in nature. In the 
case of nanocomposites, crack 
initiation may be considered 
as having a critical nucle- 
ation size. Crack initiation 
sites may well come from na- 
noparticles aggregates, as was 
noted in [47]. In [47], the 
breakdown strength of epoxy 
resin/clay nanocomposites was 
investigated, the clays being 
Cloisite 20A ( C20A) and Cloi- 
site 30B (C30B) with diffe- 
rent levels of loading. It 
was reported that the break- 
down strength depends on the 
level of loading (around 5wt% 
being the optimum) and also 
on the type of nanocomposite 
structure (intercalated nano- 
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composite or exfoliated nano- 
composite), with the exfolia- 
ted type giving the higher 
breakdown strength. Clay C30B 
gave the better results since 
it is more hydrophilic than 
clay C20B and, consequently, 
it has a greater affinity for 
epoxy resin. Epoxy pre-poly- 
mers can more easily interca- 
late into C30B clay galle- 
ries, increasing thus the in- 
terlayer spacing (more de- 
tails on intercalated and ex- 
foliated structures in [48], 
[49]). In yet another paper, 
it was reported that appro- 
priate nanoparticles, if sui- 
tably added, may enhance the 
breakdown time of conventio- 
nal epoxy resin by a factor 
of ten [50] . 

Nanosized particles were 
also discussed in [51], where 
excellent voltage endurance 
results seemed to be very pro- 
mising. Large percentages of 
silica nanoparticles (up to 
25 wt%) were reported for 
making practical coils in VPI 
(vacuum pressure impregna- 
tion) epoxy resin successful- 
ly [52]. The reported percen- 
tage seems to be excessive in 
view of previous publica- 
tions. Nevertheless, for na- 
nocomposite materials in high 
voltage machinery, the prob- 
lems remain much the same as 
with the more traditional in- 
sulations, namely that the 
heat transfer must be satis- 
factory also for the new 



materials, the mechanical 
strength must be high enough 
and the risk from PD must be 
minimized [51] . Recent work 
on mechanical properties sho- 
wed that, with epoxy resin 
and montmorillonite ( MMT ) clay 
mineral, natural frequency of 
vibration and damping factor 
of the said material increase 
by adding up to 5 wt% of nano 
clay [53] . Recent work also 
indicated that epoxy resin 
with 10 wt% Ti02 nanoparti- 
cles improved greatly the ac 
breakdown strength and the 
time to breakdown [54] . 

The percentage of nanopar- 
ticles to be included in a 
base polymer matrix depends 
on the type of nanoparticles 
as well as on the base poly- 
mer. It has been reported, 
for example, that with epoxy 
resin and layered silicate, 
just small amounts of nano- 
particles are enough for the 
improvement of partial dis- 
charge resistance, whereas in 
other publications, it was 
confirmed that only 2 wt% of 
nanoparticles is sufficient 
to improve the partial dis- 
charge resistance of poly- 
amide/layered silicate nano- 
composites [55] -[57]. Small 
amounts of nanoparticles (3 
wt% of Si02 nanoparticles) 
were also reported to improve 
the glass transition tempera- 
ture of epoxy resin in com- 
parison with the neat epoxy 
resin. This is due to the re- 
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duction of polymer chains mo- 
bility. On the other hand, in 
the same paper, the resisti- 
vity of epoxy resin with 3 wt 
% of Ti02 nanoparticles was 
lower by one order of magni- 
tude with respect to pure 
epoxy resin [58] . Paper [58] 
is a good example of the de- 
pendence of the nanocomposite 
properties on the nature of 
the added nanoparticles to 
epoxy resin. Further evidence 
as to the effect of the per- 
centage wt% of nanoparticles 
in epoxy resin is offered in 

[59] , where POSS (polyhedral 
oligomeric silsesquioxane) na- 
noparticles were added to the 
base material. Tan<5 measure- 
ments as well as thermo- 
gravimetric analysis showed 
that moderate percentages in 
wt% (between 1% and 4%) of- 
fered the best results. The 
importance of nanoparticle 
percentage and functionaliza- 
tion was also emphasized in 

[60] , where poly (butylene 
te- rephthalate) based poly- 
mers containing alumina nano- 
particles were investigated, 
as alternative to epoxy re- 
sin. It was reported that an 
optimum nanoparticle percent- 
age exists for giving a low- 
ering of the permittivity of 
the resulting nanocomposite 
as well as a lowering of 
tan5. The lowering of the 
aforementioned parameters can 
probably be ascribed to the 
restriction of polymer chain 
movement by nanoparticles due 



to the modified molecular 
structure and chain dynamics, 
which cause a strong surface 
interaction between the nano- 
particle and the polymer ma- 
trix [61] . 

Treeing effects in nano- 
composite epoxy resin propa- 
gate through the base mate- 
rial and do not go through 
the nanoparticles [44], [62]. 
Erosion depth was found to be 
minimal for a combination of 
micro- and nano- particles 

[62] . The desirable result of 
having good thermal conduc- 
tivity and low dielectric con- 
stant is more difficult to 
obtain. In [60], it was shown 
that epoxy resin with h- or 
c- boron nitride nanoparti- 
cles presents higher thermal 
conductivity at the expense 
of a higher dielectric con- 
stant, whereas epoxy resin 
with silica nanoparticles has 
a much lower dielectric con- 
stant but with a far lower 
thermal conductivity. The re- 
ported lower thermal conduc- 
tivity of epoxy resin with 
silica nanoparticles, howe- 
ver, was contradicted in [52], 

[63] . In [52], it was men- 
tioned that nanosized Si02 
particles act as barriers to 
the treeing phenomena and hin- 
der propagation. Moreover, the 
mechanical and thermal pro- 
perties are improved signifi- 
cantly, thus giving a promi- 
sing new insulation system 
with less thickness and bet- 
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ter heat transfer. The dis- 
agreement between [62] and 
[52] may be due to the dif- 
ferent processing methods as 
well as to the different size 
and/or shape of the nanopar- 
ticles . 

Differentiation between 
short-time breakdown and long 
-term failure in nanocompo- 
sites was reported in [64] . 
As the authors pointed out, 
short-term breakdown proper- 
ties depend on the applied 
voltage waveforms as well as 
the bonded region of the 
nanoparticles, whereas for 
long-term aging and failure, 
the transitional region and 
the cohesive energy density 
(CED) of the polymer matrix 
play the dominant role. They 
remarked that the percentage 
weight of nanoparticles to be 
included in a polymer matrix 
depends on the matrix itself, 
the chemical nature of the 
nanoparticles, their func- 
tionalization, their size and 
their bonding to the polymer 
matrix. For example, Ag nano- 
particles of about 20 nm in 
size mixed with epoxy resin 
at about only 0.05 wt% may 
improve the short-term break- 
down strength by about 30% 
w.r.t. the pure epoxy resin. 
Regarding nanoparticle con- 
tent and the polymer matrix, 
it was shown that with 1 wt%, 
epoxy resin nanocomposite has 
a better long-term electrical 
aging resistance than its 



polyethylene counterpart. The 
authors of [64] also remarked 
that PD resistance improves 
as the size of nanoparticles 
decreases because the proba- 
bility of electron collision 
with nanoparticles increases 
leading the electron trans- 
port to become harder. The 
latter statement agrees with 
simulation results presented 
in [42] . 

Rotating machine insula- 
tion will be better served 
with nanomaterials, if such 
materials include nanometric 
layered silica nanoparticles, 
since the latter offer better 
PD resistance and improved 
mechanical properties. As the 
design field nowadays for 
conventional machine insula- 
tion is limited to only about 
3 kV/mm [65], silica nanopar- 
ticles may help to increase 
the design field [66] . Such 
layered silica nanoparticles 
present a barrier behavior, 
rendering them interesting for 
applications . 

Possible Charging Phenome- 
na Below Inception Voltage 

The present paper did not 
deal with either the techni- 
calities of PD measurements 
in rotating machine insula- 
tion or the modeling of life- 
time of such insulation under 
a variety of simultaneous 
stresses (electrical, thermal 
and mechanical) [67] -[69]. 
The literature on such topics 
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is very rich to be dealt with 
in the present work. This pa- 
per did not deal either with 
the possibility of charging 
phenomena below the so-called 
inception voltage. Relevant 
work done in the past re- 
vealed that it is possible to 
have sudden failures in insu- 
lation systems (including 
those of rotating machines), 
even though the equipment 
passed the relevant specifi- 
cation tests [30], [70] - [76] . 
More recent research on the 
topic of charging phenomena 
below inception voltage indi- 
cated that in both base epoxy 
resin [77], and in epoxy re- 
sin with Ti02 nanoparticles 
and microparticles, charging 
phenomena were observed [78] . 
It has, however, to be noted 
that in the case of epoxy re- 
sin with nanoparticles and 
microparticles, charging phe- 
nomena below inception were 
rather sporadic. This may be 
due to the bonding strength 
between fillers and matrices, 
the interfiller space or ma- 
trix volume surrounded by 
neighboring fillers and to 
the morphology in the inter- 
filler space [79] . The whole 
subject of possible charging 
effects below the so-called 
inception voltage cannot be 
dealt with in the present pa- 
per. It is, however, a sub- 
ject which unjustly does not 
attract much attention from 
the insulation community. The 
authors intend to come back 



to this subject, possibly 
with another review paper con- 
centrated on this subject on- 
ly. Certainly, for the treat- 
ment of this question, impor- 
tant publications such as 
[80] - [83] , must be taken into 
account . 

Further Developments 

It is understandable that 
there may be alternatives to 
the epoxy resin as insulating 
material for rotating machi- 
nes. Such alternatives may be 
silicone based, resin rich in- 
sulation materials due to 
their thermal stability, fle- 
xibility, anti-vibration and 
very good electrical proper- 
ties. It would be interesting 
to see admixtures of such ma- 
terials with nanoparticles, 
as a further exploration for 
possible applications [84] . 
Furthermore, more fundamental 
research has to be performed 
regarding the combined stres- 
ses on nanocomposite poly- 
mers. Since the various stres- 
ses (electrical, thermal, me- 
chanical etc.) are applied 
not sequentially but combined 
[85] -[88], experimental work 
has to be done in this re- 
spect. In [63], detailed steps 
for future work have been pro- 
posed, such as thermal aging 
and classification tests at 
different temperatures, elec- 
trical aging - voltage en- 
durance tests at various le- 
vels resulting in lifetime 
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curves, multifactor aging as 
well as thermo-mechanical ben- 
ding endurance in parallel 
with electrical stress tests. 
On a more general basis, op- 
timization of nanocomposite 
material fabrication methods, 
a better understanding of in- 
terfaces and possible combi- 
nations of micro- and nano- 
composites may be research 
fields in the future. The 
need for nanocomposites ha- 
ving high breakdown strength, 
low thermal expansion coeffi- 
cient, high thermal conduc- 
tivity, satisfactory long- 
term aging and good with- 
standing capability to multi- 
stressing, will be even more 
pronounced in the coming 
years [89] . The variety of 
nanoparticle sizes and types, 
the variety of polymer matri- 
ces as well as the variety of 
processing methods, leaves us 
with the hope that optimal 
combinations w.r.t. the elec- 
trical, mechanical and ther- 
mal properties, may be found 
for the benefit of the elec- 
trical machines industry. 

Conclusion 

In this paper a review was 
performed for both traditio- 
nal and modern insulating ma- 
terials for rotating machine 
insulation. Traditional insu- 
lation mainly consists of mi- 
ca sheets and epoxy resin. 



with the former being the 
stronger of the two materi- 
als. Electrical trees tend to 
propagate through the epoxy 
resin and have greater diffi- 
culty in breaking through the 
mica sheets. Nanocomposites 
on the other hand offer gene- 
rally better insulating prop- 
erties. The nanoparticles that 
are dispersed in the polymer 
matrix tend to act as ex- 
tremely small barriers, pre- 
venting thus the propagation 
and growth of electrical 
trees. The performance of the 
nanocomposite polymers de- 
pends on a variety of parame- 
ters, such as, for example, 
the type of polymer matrix, 
the type of nanoparticles, 
their functionalization and 
their size. 
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RadPat4W - FLOSS for MS Windows or Wine Linux 
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Abstract 

This paper briefly highlights the features of the software 
tool [RadPat4W], named after Radiation Patterns for Windows, 
that is based on an alternative exposition of fundamental An- 
tenna Theory. This stand-alone application is compatible with 
the [Wine] environment of Linux and is part of a freeware 
suite, which is under active development for many years. Nev- 
ertheless, the [RadPat4W] source code has been now released 
as FLOSS Free Libre Open Source Software and thus it may be 
freely used, copied, modified or redistributed, individually 
or cooperatively, by the interested user to suit her/his per- 
sonal needs for reliable antenna applications, from the sim- 
plest to the more complex. 



Keywords 
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Introduction 

Useful software has to work 
exactly as someone wants, so 
the authors' group decided to 
develop its own mini-Suite of 
software tools for antenna 
applications [1] . This pro- 
ject started in the middle of 
90s, when the PC with www ac- 
cess became power enough to 
cover the increased require- 
ments of antenna analysis and 
design, as well as of their 
results presentation and dis- 
tribution. Since then, the 
development of the mini-Suite 
has been orientated towards 



the personal needs of the in- 
dividual user or of the inde- 
pendent member of a small, 
open, loosely connected group, 
like the authors' one, who is 
interested in antenna educa- 
tion, research and engineer- 
ing, i.e., a student, an edu- 
cator, a researcher, a pro- 
fessional engineer or a radio 
amateur. Such a user has 
enough bibliographical resour- 
ces provided by the Open Ac- 
cess movement, but only li- 
mited technical resources for 
construction and measurement. 
The mini-Suite is intended 
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then for the informed user 
who at least can construct an 
experimental thin-wire an- 
tenna model and at most has 
access to a VNA Vector Net- 
work Analyzer to test this 
model - by the way, nowadays, 
the cost of a certified re- 
furbished VNA is just a small 
percentage of its new price. 
For that reason, the mini- 
Suite specifically includes 
the stand-alone application 
[RadPat4W] . 

[RadPat4W] 

The active development of 
this tool attempts to bridge 
the increasing gap of today 
approximate simulation tech- 
niques, which dominate an- 
tenna applications, to clas- 
sic exact analysis methods, 
which concern the demanding 
user who wants to know what 
s/he is really doing with 
these marvellous antenna sim- 
ulators. To achieve this goal 
and facilitate the study of 
antenna application results, 
either approximated or exact, 
[RadPat4W] computes and/or 
plots the antenna geometry, 
its characteristics, as well 
as 2D main-plane cuts of its 
radiation pattern and 3D Vir- 
tual Reality objects for its 
geometry and pattern. Cur- 
rently, the tool uses by de- 
fault: (1) working formulas 
produced by the analysis 
method of the authors' alter- 
native exposition of funda- 



mental Antenna Theory [2] 
that is quickly but rigo- 
rously results in the most 
general complex vector ex- 
pression for the radiation 
pattern of any thin-wire an- 
tenna, and (2) numerical re- 
sults from approximation tech- 
niques based on the Moment 
Method implemented by the two 
antenna simulators [DA] and 
[RichWire], which are inclu- 
ded in the mini-Suite [1] . 
Finally, to support the seri- 
ous user to judge the re- 
sults, the current beta ver- 
sion of [RadPat4W] incorpo- 
rates the superposition on 
the plotted results of scien- 
tific VNA measurements with 
systematic errors first esti- 
mated by the authors in 2008 
[3], a process that is now 
accomplished semi-manually 
using a combination of other 
separate mini-Suite tools. 

Besides [RadPat4W], which 
was always distributed 
through the internet as free- 
ware, other non-commercial sof- 
ware, less related to [Rad- 
Pat4W] and from developers 
with a diverse knowledge of 
Antenna Theory, is distribu- 
ted under various terms of 
use. These ware kinds of the 
free have been exhaustively 
examined by the members of 
the USENET group [alt. comp, 
freeware] with the purpose to 
warn the candidate user about 
the actual content of the 
corresponding licenses. How- 
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ever, to the best of authors' 
knowledge, it seems that the- 
re is still no completely 
FLOSS for reliable antenna 
applications. Therefore, with 
the aim to further encourage 
the independent user to tweak 
[RadPat4W] according to her/ 
his personal needs or even to 
be involved in this modern, 
most promising, cooperative 
activity of the FLOSS move- 
ment, the authors decided 
lately to release the entire 
source code of [RadPat4W] un- 
der the approved, by the OSI 
Open Source Initiative, MIT 
License . 

The source code, now in 
version 4.4 with help in ver- 
sion 1.0, is developed from 
scratch, without using any 
other code, in MS Visual Ba- 
sic 6 SP6 for 32-bit MS Win- 
dows and the executable, 
which is also compatible with 
the [Wine] environment of Li- 
nux, needs about 8.5 MB of 
free hard disk space for its 
installation. The application 
usability has been multiple 
checked during antenna cour- 
ses and theses elaboration, 
thus its source code is in a 
mature state for a long time 
now although, from time to 
time, new features are added 
to it. The code is available 
for download from authors' 
group website "http://www.ant 
ennas.gr/floss" or GoogleCode 
repository at "http://code.go 
ogle . com/p/rga/" . 



The features of [RadPat4W] 
can be divided according to 
their functionality in two 
groups: (1) pattern computa- 
tion and plotting using work- 
ing formulas from Antenna 
Theory, and (2) pattern plot- 
ting using numerical results 
from antenna simulators. To 
exemplify these features by 
examples, a number of antenna 
education, research and engi- 
neering applications, from 
the simplest to the more com- 
plex, are presented in the 
following . 

Working Formulas 

Fig. 1 shows the software 
application form of [Rad- 
Pat4W] for the three main- 
plane cuts of the E-norma- 
lized radiation pattern of a 
linear, center-fed, standing- 
wave dipole, which is paral- 
lel to z-axis and has a 
length of 2.35A, where the 
length is the only one input 
parameter with values in the 
range [0.001, 10]A. To over- 
come the practical constrains 
of the limited number of 
screen pixels that obscures 
the detailed view of zero-E 
directions, which determine 
the radiation pattern lobes, 
a useful feature has been in- 
troduced in all application 
forms that is the magnifica- 
tion of the pattern up to six 
times. Each magnified diagram 
shows the zero-E directions 
by magenta colored radial li- 
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nes on the screen. The [Max- 
Zero-D] button opens the win- 
dow of Fig. 2 in which the 
computed directions of maxima 
and zeros, as well as, the 



directivity and the maximum 
value of pre-normalized E ra- 
diation pattern, are shown 
[ 2 ]- 



\ Radiation Patterns 




Fig. 1: [One Dipole on z] : A magnified main-plane cut 




Fig. 2: [One Dipole on z] : 
Maximum, zeros and Directivity 



Fig. 3 shows the three 
main-plane cuts of radiation 
pattern for the same dipole 
but in the space direction 
described by the input data 
of its unit directional vec- 
tor: (0.000, 0.707, -0.707). 

The [Zero] button opens the 
window of Fig. 4, where the 
computed directions of zeros 
on the three main-plane cuts, 
as well as on a plane that 
contains the dipole axis, are 
shown . 
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Fig. 3: [One Dipole in Space]: Three different main-plane cuts 




Fig. 4: [One Dipole in Space]: Zero pattern directions 



Fig. 5 shows the input 
data for a uniform linear ar- 
ray of exactly parallel lin- 
ear standing-wave dipoles, in 
two frames for the generator 
dipole with complex vector 
pattern G and the array of 
isotropic sources with com- 
plex number pattern A, re- 
spectively. The first frame 



defines the length and the 
direction of the generator or 
reference dipole. The second 
frame defines the geometrical 
and electrical input data for 
the array: the number of iso- 
tropic point sources, their 
constant phase difference in 
degrees, their constant equi- 
distance per wavelength, and 
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the unit directional vector 
of array axis. The shown val- 
ues are for an array of 3 
dipoles, 2A each, in the di- 
rection of +y axis, with 
phase difference of -144°, 
0.4A apart, on the array axis 
direction of +x. The [Show 
Graphics] button opens the 
window of Fig. 6, in which 
the dipole array complex vec- 
tor pattern E = AG, i.e., the 
Principle of Radiation Pat- 
terns Multiplication, is shown 
in absolute 3D form: the nor- 
malized norm pattern ||E|| 
results as product of the 
normalized absolute pattern 
| A | by the normalized norm 
pattern ||G|| and by a non- 
shown constant spherical pat- 
tern | A | max | | G | | max / | | E | | max > 1 
[2] . In the [Directivity] 
frame of the window, the di- 
rectivities of the array Da, 



of the generator Dg, and of 
the dipole array D are shown. 
The [Max-Zero of Array Fac- 
tor] button opens the window 
of Fig. 7, in which the di- 
rections of zero-A and |A| max 
and the values of |A| max , || 



and 



are shown. 



The buttons [A-3D], [G-3D] 

and [E-3D] produce the re- 
spective 3D Virtual Reality 
radiation patterns, which are 
shown in Fig. 8 as three 
screen captures of the free 
Platinum Worldview VRML vie- 
wer plug-in for MS Internet 
Explorer. By the way, the 
contemporary free VRML Cor- 
tona viewer is available for 
a number of other web brow- 
sers too, under MS Windows, 
while under Linux the most 
appropriate add-on for ice- 
weasel (Mozilla Firefox) is 
the FreeWRL VRML viewer. 




Fig. 5: [Array Data]: Uniform Linear Array input 
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Fig. 6: [Array Graphics]: 9 main-plane pattern cuts, 2 zeros 
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Fig. 7: [Max-Zero of Array Factor]: Computed results 




Fig. 8: Virtual 3D Principle of Radiation Pattern 
Multiplication 
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Antenna Simulators 

The second group of [Rad- 
Pat4W] features concerns its 
ability to plot antenna geo- 
metry and patterns from the 
clear text data files such 
those of the two antenna si- 
mulators of the mini-Suite: 
[DA] and [RichWire] . The sour- 
ce code of [DA], now in ver- 
sion 1.0.8, has been written 
in Compaq Visual Fortran 6.1 
as Quick-Win 32-bit applica- 
tion for MS Windows, and the 
executable, which is also com- 
patible with the [Wine] envi- 
ronment of Linux, can be eas- 
ily installed on a PC with 
free hard disk space of only 
about -500 KB. 

In essence, this applica- 
tion is a restricted varia- 
tion of [RichWire], which is 
a fully analyzed, corrected 
and redeveloped edition of 
the original Moment Method 
thin-wire computer program by 
J . H . Richmond, available in 
the public domain by NASA 
since 2005 [4] . [DA] is used 
for antenna simulation by 
half-wave dipoles, with just 
one active. The program re- 
quires an input data file to 
derive three output data fi- 
les. All these data files are 
used by [RadPat4W] . In Fig. 
9, the simplest data for only 
one dipole in space are 
shown. [RichWire] data files 
are similar. The usability of 



both antenna simulators has 
been also multiple checked. 
The simulators are available 
as freeware from the men- 
tioned repositories. 

In Fig. 10, the [RadPat4W] 
engineering application for a 
commercial VHF Yagi-Uda an- 
tenna is shown. The [Antenna] 
button reads the antenna geo- 
metry [RichWire] data input 
file. The [geo.wrl] button 
produces the 3D Virtual Rea- 
lity antenna geometry and si- 
multaneously opens the [GL Vie- 
wer for Mathematica] for an 
immediate view [5]. In Fig. 
11, the drawn results produ- 
ced by [RadPat4W] are shown 
for an educational applica- 
tion of a flat airplane mo- 
deled with non-overlapping A/2 
dipoles in [DA] [6]. Fig. 12 
illustrates the drawn results 
produced by the currently be- 
ta version of [RadPat4W] for 
a research application of a 
constructed Hentenna model, 
simulation designed with [Rich 
-Wire] and measured with a 
VNA system [7] . 

[RadPat4W] Development Plans 

Scheduled expansions of 
[RadPat4W] include the fol- 
lowing facilities, which are 
already available in other 
mini-Suite tools: (1) choice 
of other plane- or conical- 
cuts, (2) key-in of any exact 
analysis working formula E(0, 



FUNKTECHNIKPLUS # JOURNAL 



Vl-40 



ISSUE 5 - YEAR 2 




ANTENNA RADIATION PATTERNS: RADPAT4W - FLOSS FOR MS WINDOWS OR WINE 



cp), (3) selection of the Z,- , 
0-, and cp- plotting step, (4) 
automation of fine Cartesian 
pattern plotting, (5) super- 
position of VNA measurements 
on 3D Virtual Reality pat- 
terns, such as that in the 
left part of Fig. 13 [8], and 
(6) superposition of VNA mea- 
surements with their diffe- 
rential error cloud [3] on 2D 



plots, such that in the right 
part of Fig . 13 [9] . 

Any other expansion of the 
freely available code is of 
course welcomed. In authors' 
group, there are no plans in 
the near future to upgrade 
the mini-Suite to an inte- 
grated environment for its 
tools . 




Fig. 9: [Da]: Input and output data for the simplest case 
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Radiation Patterns - [Radiation Pattern By Data] 
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Fig. 10: [RichWire] : A commercial VHF Yagi-Uda antenna 




Fig. 11: [RadPat4W] : Results for a flat airplane model 
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Fig. 12: [RadPat4W] : Design, Construction and Measurement 




Fig. 13: 3D VNA measurements and their 2D error cloud 
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Antennas Research Group, Austria - Hellas [2, 3] 

Abstract 



This paper introduces the FLOSS Free Libre Open Source 
Software [VEMSA3D] , a contraction of "Visual Electromagnetic 
Simulator for 3D Antennas", which are geometrically modeled, 
either exactly or approximately, as thin wire polygonal stru- 
ctures; presents its GUI Graphical User Interface capabili- 
ties, in interactive mode and/or in handling suitable formed 
antenna data files; demonstrates the effectiveness of its use 
in a number of practical antenna applications, with direct 
comparison to experimental measurements and other freeware 
results; and provides the inexperienced user with a specific 
list of instructions to successfully build the given source 
code by using only freely available IDE Integrated Develop- 
ment Environment tools-including a cross-platform one. The 
unrestricted access to source code, beyond the ability for 
immediate software improvement, offers to independent users 
and volunteer groups an expandable, in any way, visual an- 
tenna simulator, for a genuine research and development work 
in the field of antennas, adaptable to their needs. 

Keywords 



FLOSS, antennas, modeling, simulation 



Introduction 

A lot of amazing visual EM 
software simulators, both com- 
mercial and freeware, exist 
for many years now. However, 
to the best of authors' know- 
ledge, none of these simula- 
tors has been ever released 
under an Open Source license- 
paid, free gratis or free li- 
bre. This situation, as it 
was discussed recently, pre- 



cludes the independent users 
and volunteer groups with 
limited resources, like the 
authors' nonprofit group, from 
the genuine, state-of-the- 
art, scientific research, sin- 
ce, neither the code improve- 
ment, nor its adaptation to 
specific needs, is possible 
[ 1 ], [ 2 ]. 

For that reason, the au- 
thors decided to develop and 
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release under the GNU public 
license their own Visual Elec- 
tromagnetic Simulator for 3D 
Antennas [VEMSA3D] , although 
it certainly is a less well- 
equipped software application 
for the moment-but still a 
fully expandable one. That re- 
lease is now perfectly per- 
missible, because in 2005 
NASA released in the public 
domain the FORTRAN source co- 
de of the well-known MoM Me- 
thod of Moments Thin-Wire Com- 
puter Program, by J. H. Ri- 
chmond [3], [4]. This is ex- 
actly the code on which the 
authors' group based its free- 
ware simulators: [RichWire], 
a CLI Command Line Interface, 
and [DA], a MS Quick-Win FOR- 
TRAN derivative of it. These 
two simulators are under un- 
interrupted revision, impro- 
vement, expansion, and rede- 
velopment, since a long time 
ago [1], [2]. Therefore, the 
[Rich-Wire] FORTRAN code was 
translated-entirely, line-by- 
line, without using any paid 
or free translator-to C++, to 
form the core of scientific 
EM computations in [VEMSA3D] . 

On the other hand, the au- 
thors' group requirements for 
scientific accuracy in visual 
representation of the produ- 
ced EM simulation antenna re- 
sults from its simulators, we- 
re already enforced the soft- 
ware development of the Vir- 
tual Antennas, that is the 
Virtual Antennas laboratories. 



in VRML [5], [6], the FLOSS 

application [RadPat4W] for an- 
tenna radiation pattern pre- 
sentation, in MS VB6 [2], 
[7], as well as, the recently 
developed visual antenna ap- 
plication for the Wolfram De- 
monstrations Project in Ma- 
thematica [8] . The visualiza- 
tion ideas implemented in the 
aforementioned graphics ap- 
plications were also expres- 
sed-from the scratch-in C++, 
to use the Open Source cross- 
platform [wxWidgets] library 
with OpenGL and form the core 
of scientific EM graphics in 
[VEMSA3D] [9] . 

The authors, having taken 
into account that no famil- 
iarization with software use 
is possible without getting 
hands-on experience, restrict 
themselves to a brief discus- 
sion of the current [VEMSA3D] 
characteristics. The code, the 
antenna applications data, as 
well as, any other informa- 
tion, referenced open bugs or 
future code releases, will be 
always available in authors' 
group repositories, at 
"www. antennas . gr /floss" 
and in GoogleCode website, at 
"http : //code . google . comp/p/rg 
a/" . 

GUI Interactive Mode of 

Operation 

It is assumed that the 
user has already some experi- 
ence in the sketching of a 
polygonal wire outline model, 
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for an antenna under consid- 
eration, consisting of num- 
bered wire segments and no- 
des, with their 3D space co- 
ordinates, including the po- 
sitions of any antenna cir- 
cuit elements, that is volt- 
age generators and lumped 
loads, and s/he wants then to 
key-in these model and cir- 
cuit data into [VEMSA3D] us- 
ing its GUI Graphical User 
Interface in interactive mo- 
de . 

The GUI main window is 
shown in Fig. 1, while Fig. 2 
shows the menu items along 
with their available submenu 
options numbered from 1 to 7. 

The GUI itself is divided 
in three panels named [An- 
tenna in space (3D)], [Anenna 
Elements], and [Data-sets]. 



The function of each panel is 
briefly described in the fol- 
lowing . 

The [Antenna in space 
(3D)] viewing panel is used 
to project all the generated 
3D and 2D graphics. By de- 
fault, the application starts 
with a simple linear dipole 
antenna in space, as it is 
shown in Fig. 1 from a view- 
point on the diagonal of tri- 
hedral angle OXYZ, where the 
usual letters associated with 
the Cartesian axes do not ex- 
ist. Instead, a one-to-one 
correspondence implicitly ex- 
ists between (X, Y, Z) axes 
and (R, G, B) colors [5]. The 
3D image of any antenna can 
be manipulated through the 
[Study] options, as shown in 
Fig. 2.3. 




Fig. 1: GUI: The main window 
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Study Graph Simulation 
Tool Bar 2 

' Status Bar 
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• Antenna Elements Panel 



Status Window 
Results Window 
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ZOX View 



Options Help 
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Help 

Viewing Options... 0 

Dataset Simulation Options... 



Options MIS 



About. .. 2 1 



Simulation Options Help 
• Antenna in space (3D) 4 

Antenna segments in space (3D) 

E Radiation Pattern in space (3D) 

E Radiation Pattern on 3 main planes (2D) 

E Radiation Pattern on xOy plane (2D) 

E Radiation Pattern on yOz plane (2D) 
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E phi Radiation Pattern in space (3D) 
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Fig. 2: Main window: The unfolded menu items 



The [Antenna Elements] pa- 
nel contains four tabs, cor- 
responding to model and cir- 
cuit data, in four self-ex- 
plained building blocks. In- 
stead of attempting to de- 
scribe the block data in gen- 
eral the presentation conti- 
nues with the default data: 

[Nodes] : 2, with coordi- 
nates 1: (0.2, 0.1, 0.25), 2: 
(0.2, 0.1, -0.25), in meters 
[m] or in terms of wavelength 
A [wl]. 

[Wires] : 1, starting at 
node: 1, ending at node: 2, 
that is a directed segment. 

[Segments] is either: (a) 
a positive integer, (b) a 
zero or (c) a number between 
0.00 and 0.25, to respec- 
tively define the division of 
this wire in segments: (a) in 
the indicated number of seg- 
ments, (b) in a calculated 
number of segments of length 
no longer than 0.05A (that is 
10 segments, in this case), 



(c) in segments with length 
no more to that number. No- 
tably, no segment can be 
longer than 0.25A [3]. 

[Generators] : 1, connected 
in the wire:l, at a distance 
of 0.5 times this wire 
length, from its starting no- 
de, with rms value of 1^0° 
[V] . There can be only one 
generator on each wire. 

[Loads]: Has the same stru- 
cture as [Generators] . Resis- 
tors, inductors and capaci- 
tors are inputted in [Q], [H] 
and [F], respectively. It is 
empty in this case, since no 
lamped loads exist. 

Notably, the GUI has been 
developed in such a way that 
efficiently supports the in- 
teractive handling of the an- 
tenna modeling by the user on 
the screen. Thus, the selec- 
tion of an element in any of 
these four tabs simultane- 
ously highlights this element 
on the viewing [Antenna in 
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space (3D)] panel-if the [Se- 
lect] option under [Study] 
(Fig. 2.3) has been already 
chosen by the user, and this 
is the default selection in 
any case-and conversely: the 
selection of any element on 
this viewing [Antenna in spa- 
ce (3D)] panel highlights the 
element in its table. Conse- 
quently, a full control on 
the antenna modeling struc- 
ture is achieved. 

The [Datasets] panel may 
contain independent sets of 
model and circuit data for va- 
rious antennas. Such a Data- 
set has its own [Dataset Si- 
mulation Options], which are 
chosen through submenu of 
Fig. 2.6, as the default se- 
lection of them is shown on 
the left part of Fig. 3. 

A [Wire conductivity] of -1 



corresponds to perfect wire 
conductivity, while the [Spa- 
ce conductivity] of 0 and the 
[Space dielectric constant] 
of 1 defines the Free Space 
EM environment. The [Integra- 
tion steps] accepts a posi- 
tive number for the steps of 
approximated integration used 
in the MoM impedance computa- 
tions-a zero means closed- 
form integration. 

Notably, the [Maximum seg- 
ment length] affects the di- 
vision of wires: the smaller 
this number is, the more seg- 
ments will be considered. It 
is currently known that there 
is an implied, still program- 
matically unimproved, priori- 
ty of the number of segments 
defined in [Wires] tab, over 
this selection-that is a bug, 
which may crash the applica- 
tion . 



Dataset Simulation Opti 


ons X 


Operating frequency : 


300 


[MHz] 


Wire radius : 


0.001 


[m] 


Wire conductivity : 


1 _1 


[MS/m] 


Speed of light : 


3e4008 


[m/sec] 


Space conductivity : 


0 


[S/m] 


Space dielectric constant : 


1 


]B 


Maximum segment length : 


0.05 


Pambda] 


Integration steps : 


0 


]0 


3D angle step : 


5 


[deg.] 


2D angle step : 


1 


[deg.] 



OK ~j | Cancel | | Default ] 



Status 



20:20:35.453 


Simulation started. 


20:20:35.453 


Running Dataset : Dataset 


20:20:35.687 


SEG started. 


20:20:35.703 


SEG finished. 


20:20:35.703 


SORT started. 


20:20:35.750 


SORT finished. 


20:20:35.750 


SGANT started. 


20:25:32.718 


SGANT finished. 


20:25:32.718 


GANT1 started. 


21:03:11.906 


GANT1 finished. 


21:03:11.906 


PATTERNS started. 


21:03:22.500 


PATTERNS finished. 


21:03:22.500 


Simulation finished. 


21:03:22.500 


Total time: 00:42:47.047 



Close 



Fig. 3: Windows: [Dataset Simulation Options] - [Status] 
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The EM simulation of an an- 
tenna starts by selecting [Si- 
mulate Current Dataset], un- 
der [Simulation] (Fig. 2.5), 
while its progress is shown 
in the [Status] window, such 
as the one at the right of 
Fig. 3, which concerns the 
default dipole. Notably, mul- 
tiple Datasets can be simu- 
lated, one after the other, 
through the [Simulate Marked 
Datasets] shown in Fig. 2.5. 

The course of the simula- 
tion may be roughly described 
in five steps: 

(1) Wires are automatically 
segmented, using the provided 
parameters, to produce the fi- 
nal model structure of points 
and segments 

(2) Model structure is ana- 
lyzed and adjacent segments 
are combined to form dipole 
current modes with sinusoidal 



distribution 

(3) Impedances are calculated 
and the MoM algebraic system 
of equations is formed 

(4) The system of equations 
is solved and the segment cur- 
rents are calculated 

(5) 3D and 2D radiation pat- 
terns are calculated, as well 
as other useful results like 
those shown in Fig. 4. 

In Fig. 5, the [Viewing 
Options] under the [Options] 
(Fig. 2.6) is given with all 
of the available options for 
the graphics. Through [Graph] 
menu item (Fig. 2.4), a vari- 
ety of 3D and 2D plots, which 
include either the normalized 
electric far field E radia- 
tion pattern or the relative 
radiation intensity U in dB, 
as well as their 0 and cp 
parts, can be illustrated. 




Fig. 4: Window: [Results] - Table: [Power] 
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The [Antenna segments in 
space (3D)] illustrates the 
final antenna structure with 
points, segments and segment 
currents in amplitude and 
phase, as shown in Fig. 6, 
where the resulting current 
distribution on default di- 
pole it seems to be sinu- 
soidal, that is as it was ex- 
pected to be. The well-known 



2D intersections of the 3D 
radiation patterns by the 
three main planes yOz, xOy, 
zOx, as well as these 3D pat- 
terns themselves are also 
shown for the default dipole. 

Notably, modeled antennas 
can be imported and exported 
in [RichWire] data file for- 
mat (Fig 2.1) . 




Fig. 5: [Viewing Options] 




Fig. 6: GUI: Default dipole results 



Build under MS Windows XP 

To build [VESMA3D] MS Win- 
dows executables, the use of 



IDE CodeBlocks with GCC/Mingw 
is suggested [10] : 
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Set below as X:\ the system disk, e.g. X:\ -> C:\ 

1 Download (~20MB) , Install (~110MB) : 

http : //sourceforge . net/pro jects/codeblocks/f iles/Binaries/ 
[View all files] [Binaries] [8.02] 
codeblocks-8 . 02mingw-setup . exe 
Select: [Full], [Run], Compiler: [GCC/Mingw] 

2 Download wxWidgets (~20MB) , Extract (~120MB) : 

http : //sourceforge . net/pro jects/wxwindows/f iles/2 . 8 . 11/ 
[View all files] wxWidgets-2 . 8 . 11 . zip 
wxWidgets-2 . 8 . 11 . zip-> X:\ [Yes to All] 

3 IDE CodeBlocks, Open: 

X : \wxWidgets-2 . 8 . ll\include\wx\msw\setupO . h 
Replace, at line 1006: 

#def ine wxUSE_GLCANVAS 0 -> 1 
Repeat it, to file: setup. h 

4 At system variable [path] : 

[My Computer] [Properties] [Advanced] 

[Environment Variables] [System Variables] 

Path [Edit ] [Variable value], Add: 

" ; X : \Program Files\CodeBlocks\MinGW\bin; " 

5 Command Prompt : >cd\ 

>cd X : \wxWidgets-2 . 8 . ll\build\msw 
>mingw32-make -f makefile. gcc MONOLITHIC=0 
S HARED =0 UNICODES USE_OPENGL=l 
BUILD=release 

Re-Command, with option: 

BUILD=release -> BUILD=debug 

6 Download (~100KB) , Extract (~700KB) : 
http : //code . google . com/p/rga/ downloads/ 

VEMSA3D_source_l . zip -> X : \ 

7 IDE CodeBlocks, [File] [Open] : 

X : \VEMSA3D\build\win-cb\VEMSA3D . workspace 

* If system disk is not C:\ then: 

(#) [Project] [Build options. . .] [VEMSA3D] 

* [Search directories] Correct in each of: 

* [Compiler] , [Linker] , [Resource compiler] 

* C:\ -> X : \ [Yes] 

* Repeat from (#) , for [Release] 

* Repeat from (#) , for [Debug] 

[Build] [Select target] [Release] [Build] 

[Build] [Select target] [Debug] [Build] 
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This process results in 
[VEMSA3D] executables of ~35 
MB Debug and of ~5MB Release 
versions, which run under NT4, 
W2K and WXP, at least. 

Alternatively, the use of 
the freeware IDE MS Visual 
Studio Express with C++ 2008 
Compiler is also suggested. 
For that, first download MS 
Visual Studio 2008 Express 
iso-image from Microsoft web- 
site, burn it into a CD-R, 
and install it. Then, down- 
load from our repositories 
the file [VEMSA3Dfiles-4win- 
cb.zip], extract it, and fol- 
low the included setup in- 
structions. This process re- 
sults in [VEMSA3D] executa- 
bles of -5MB Debug and of 
~2MB Release versions, which 
run under W2K and WXP, at 
least, but definitely do not 
run under NT4 . 

Practical Antenna Applica- 
tions and Results 

This section presents the 
implementation of [VEMSA3D] to 
produce eight antenna models, 
from simple to more compli- 
cated structures. The input 
data were imported through 
[Import Geometry] of the [Fi- 
le] menu (Fig. 2.1). 

Fig. 7 shows an array of 2 
dipoles for operation at the 
frequency of 1111 MHz, dis- 
tanced by 0.85A, constructed 
by bare copper wire of 1 mm 
(0.0037A) radius, and mea- 
sured. In the same figure, 



the 3D radiation intensity 
pattern and its 2D main plane 
cuts are shown [11] . The con- 
tinuous line represents [VEM- 
SA3D] results, the dashed li- 
ne, analytically produced pat- 
terns, and the dots, measure- 
ments made using the authors' 
group VNA system [12], [13]. 

Fig. 8 illustrates the re- 
sults for a constructed and 
measured improved Hentenna mo- 
del, at 1110 MHz, with height 
of A/2, width of A/6, and 
with the active element at a 
distance 7A/12 from its bot- 
tom [14] . 

The third antenna consists 
of a A/4 monopole (5.83 cm) 
at 1286 MHz over a circular 
counterpoise of 14 cm radius. 
The monopole was constructed 
by 2 mm diameter copper wire 
and the counterpoise has been 
printed on a circuit board of 
29 cm x 29 cm, as a circle 
with four radials of 2 mm 
width [15]. In Fig. 9, the 
left 3D and 2D patterns are 
for vertical polarization, 
while the right ones for ho- 
rizontal. For a better repre- 
sentation, the patterns have 
been normalized. 

A somehow more complicated 
antenna is shown in Fig. 10. 
The left 3D and 2D patterns 
correspond to the vertical po- 
larization of a monopole over 
circular counterpoise with 16 
radials, while the right co- 
lumn patterns are the corre- 
sponding 3D and 2D patterns 
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for the same antenna, but with 
an additional grounded disc 
below it, and without elec- 
trical connection to the coun- 
terpoise antenna, at a dis- 
tance of 0.9 cm. The grounded 



disc was constructed from thin 
copper sheet and was modeled 
in [VEMSA3D] with 5 concen- 
tric circles and 64 radials 
[15]. 




Fig. 7: Array of 2 dipoles: 73 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 8: Hentenna: 109 dipole current modes 
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Fig. 9: Monopole over a counterpoise with 4 radials: 
167 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 10: Monopole over a counterpoise with 16 radials: 
419 dipole current modes. 

Monopole over a counterpoise with 16 radials over 
a grounded disc: 2594 dipole current modes 
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In Fig. 11, the results of 
[VEMSA3D] for a commercial 
UHF antenna are shown [16] . 
The antenna model is presen- 
ted in detail and separately, 
for the active element sys- 
tem, as well for the wire- 
frame reflector. The patterns 
correspond to the center fre- 
quency of its operation at 
650 MHz. 

Fig. 12 illustrates the 
model of a small jet airplane 
from the well-known freeware 
[4NEC2X] antenna simulator 
[17] . Simulation was carried 
out at the frequency of 10 
MHz with the same number of 
points and segments in both 
simulators. There is a good 
agreement between them for 



the radiation intensity pat- 
terns although some deviation 
in input impedance and direc- 
tivity is observed. 

Fig. 13 shows the results 
for a horn antenna with di- 
mensions proposed by K. Pitra 
and Z. Raida [18]. A small 
bow-tie feeder with a trian- 
gular perimeter of 0.47A and 
flare angle 37.50°, is con- 
sidered. The antenna was ini- 
tially simulated at the fre- 
quency of 40 GHz in both 
[4NEC2X] and [VEMSA3D] simu- 
lators, with 1205 segments to 
be consistent with the re- 
stricted maximum number of 
1500 segments of [4NEC2X] . 
There is a good agreement be- 
tween the produced results. 




Fig. 11: Typical commercial 
723 dipole current 



TV UHF antenna: 
modes . 
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[VEMSA3D] 



Fig. 12: Jet airplane modeled in [4NEC2X] and [VEMSA3D] : 
391 dipole current modes 



Finally, the most compli- 
cated antenna model, which is 
presented, corresponds to the 
same horn antenna at the same 
frequency of 40 GHz, with a 
dense wire-frame consisting 
of 3266 points and 4362 seg- 
ments. The model is shown in 
Fig. 14. The total number of 
dipole current modes is 5458, 



as shown in Fig. 15. The pro- 
cess time ranges from ~40 min 
in an AMD Phenom X2 550 3.11 
GHz CPU to -220 min, in an 
Intel Pentium 4 1.7 GHz CPU. 
Fig. 16 illustrates the re- 
sulting radiation patterns 
and Fig 17 the current ampli- 
tude and phase. 
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Ji 

w 

[4NEC2X] 




[VEMSA3D] 




Impedance | 


30 -i 


0.06 


S.W.R. 50 j 


1.67 




Frequency 


1 40000 


Mhz 


Wavelength 


| 7.e-3 


mtr 


Norm-All : 


: 6.53 dBi 






Resistance [Ohm] 

4 . 673E+001 

Reactance [Ohm] 

-3 . 929E+000 

VSWR [] 

1.1115 

Directivity 
6. 702E+000 [dBi] 



Fig. 13: Horn antenna with bow-tie feeder modeled in 
[4NEC2X] and [VEMSA3D] : 1370 dipole current modes 




Fig. 14: Horn antenna with bow-tie feeder: 4362 segments 
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Points Segments Dipole Current Modes ; input Sources Power 





Point 1 


Point 2 


Point 3 


Segment A 


Segment B 


Amplitude [A] 


Phase [deg.] 


5455 


1121 


3263 


1122 


4319 


4320 


1 . 053E-002 


-170.85 


5456 


1121 


3264 


1123 


4321 


4322 


1 . 042E-002 


-172.78 




1122 


3265 


1123 


4323 


4324 


1 . 615E-004 


-127.73 


^5458 


) 1121 


3266 


1118 


4325 


4326 


2 . 602E-002 


12.29 



Points Segments Dipole Current Modes 



Total Input Power (Real) 


2 . 543E-002 [W] 


Total Input Power (Imag) 


-5 . 540E-003 [W] 


Total Radiated Power 


2 . 543E-002 [W] 


Total Dissipated Power 


0 . 000E+000 [W] 


Power Efficiency 


100% 


Power Gain 


5 . 719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 


Directivity 


5.719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 



Fig. 15: Horn antenna with bow-tie feeder: 
5458 dipole current modes 




Fig. 16: Horn antenna with bow-tie feeder: 
Radiation intensity 
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Fig. 17: Horn antenna with bow-tie feeder: 
Current Amplitude and Phase 



Conclusion 

The first stable version 
of a visual EM simulator has 
been developed and released 
as FLOSS. The computational 
results of its use were found 
to be in good agreement with 
experimental measurements as 
well with the comparable free- 



ware simulator [4NEC2X] . Sin- 
ce the number of possible mo- 
difications and additions to 
the attributes of [VEMSA3D] 
seems to be endless, no at- 
tempt will be made to suggest 
a particular direction for its 
future development: any con- 
tribution from the antenna 
community is very welcomed. 
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Abstract 

The aim of this paper is the study of various parameters 
affecting the lifetime of transformer oil in transformers of 
150/20 kV. Fifty (50) samples of oil were taken from such 
transformers in the major Athens area, Greece. The parame- 
ters investigated - according to international standards - 
were breakdown strength, oil color, humidity, interfacial 
tension and tan5. Thermal and mechanical stresses have as 
result the oxidation of transformer oil and the deteriora- 
tion of its insulating properties. Humidity and foreign par- 
ticles also consist factors contributing to the lowering of 
the breakdown strength of transformer oil. In most examined 
samples the breakdown strength and tan5 were satisfactory. 
It is shown that the slightest contact with the atmospheric 
air may affect humidity. A color index of 3 does not neces- 
sarily mean that the oil is bad. Foreign particle presence 
combined with humidity may decrease the interfacial tension. 
Generally speaking, the 50 investigated transformers were in 
a satisfactory state and none of them was required to have 
its oil replaced. A main conclusion of this work is that we 
should not base our judgment about the oil quality on only 
one or two parameters but on a combination of more parame- 
ters . 



SATURDAY 3 JANUARY 2015 



Vl—27 



FUNKTECHNIKPLUS # JOURNAL 




M.G. DANIKAS, S. GEORGHIOU, I. LIAPIS, A.B.B.ABD. GHANI 



Keywords 

Transformer oil, partial discharges, diagnostic tools, 
breakdown voltage, dielectric strength 



Introduction 

Transformer oil is a very 
important component of a 
transformer. It must have 
good thermal and insulating 
properties [1] . The insulat- 
ing oil is used for providing 
insulation between the live 
parts of the transformer and 
the grounded parts as well as 
for carrying out the heat 
from the transformer to the 
atmosphere [2] . Faults in 
transformers are rare (l%-2% 
per year in power transmis- 
sion and distribution sys- 
tems), but they can be very 
costly in terms of economic 
and technical consequences. 
Faults can due, among others, 
to lightning and switching 
overvoltages, insulation fai- 
lure, humidity, foreign par- 
ticles and bad maintenance. 
The latter is a significant 
factor for the lifetime of a 
transformer [3], [4]. The main 
factors which may accelerate 
the ageing of the oil are hu- 
midity, temperature varia- 
tions, oxidation and foreign 
particles. Various factors af- 
fecting the breakdown strength 
of transformer oils have been 
discussed and analyzed [5], 
[6] . There is no single mea- 
surement which can deliver 
enough information as to the 



ageing and/or deterioration of 
transformer oil, mainly be- 
cause transformer insulation 
is a dynamic system, in which 
e.g. humidity may migrate 
from the oil to paper insula- 
tion and from paper insula- 
tion back to the oil [3], 
[7] . With this in mind, a va- 
riety of diagnostic methods 
were employed in order to stu- 
dy the state of fifty trans- 
formers of 150/20 kV. The who- 
le work was carried out with 
the aid of the Public Power 
Corporation (PPC) Transformer 
Division in Athens, Greece. 
The fifty investigated trans- 
formers were from the major 
Athens area. 

Diagnostic Methods 

Warning signs about the 
state of a transformer are, 
among others, a big increase 
of partial discharges (» 
2500 pC), a visible deterio- 
ration because of foreign me- 
tallic and carbon particles, 
the presence of humidity in 
the solid insulation about 3- 
4% and the presence of slud- 
ge. The latter is the last 
visible state of deteriora- 
tion. Experience indicates 
that the breakdown behavior 
and breakdown voltage are de- 
termined much more from the 
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above mentioned factors than 
that of pure insulating li- 
quids [4], [8], [9] . 

Several diagnostic methods 
were used in order to see the 
quality of the transformers 
in question. The characteri- 
zation of the oil color (DIN 
51517 - ASTM 155) was per- 
formed through a device (chro- 
mometer) including standard 
glass disks and two glass 
jars with lid. The control of 



breakdown voltage was mea- 
sured by a typical Foster 
test cell, according to IEC 
156/95 (Fig. 1). 

The control of humidity in 
the oil was measured by a 
Metrohm - 684 KF Coulometer, 
which consisted of a glass 
container with a stirrer ti- 
tration in which the reagent 
from container storage is ad- 
ded . 




Fig. 1: The test cell for breakdown voltage measurements 
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The device is fully auto- 
mated and once the experi- 
menter gives the settings, it 
measures the moisture content 
of the oil. The measurements 
were performed according to 
IEC 814. The control of in- 
terfacial tension (ASTM D971 
- 91) was performed via a 

tensimeter, which gives the 
value in dynes per centimeter 
in a direct reading. 

The device that performed 
measurements of tan5 and of 



resistivity, is the BAUR-DTL 
fully automated device for 
measuring dielectric losses 
of oils. Such a system has a 
fully automated process for 
measuring dielectric loss, re- 
lative dielectric constant 
and resistivity (Fig. 2). The 
measurements were performed 
according to IEC 247. The 
density of oil was performed 
according to DIN 51517, with 
the aid of a pipe of 250 ml, 
an electronic thermometer and 
a glass cylinder. 




Fig. 2: Device for measurement of tan5 and resistivity 
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It is true that no single 
diagnostic method can give 
full information as to the 
state of the transformer oil. 
The aforementioned methods 
may give a better picture of 
its state. 

Results 

The sampling was done ac- 
cording to specification ASTM 
D 923. Sampling should take 
place in clean conditions (ab- 
sence of humidity and pollu- 
tion), suitable glass vessels 



should be used and the latter 
should be kept clean and her- 
metically closed. Every sam- 
ple should be kept away from 
light according to VDE 0370/ 
9.61. Every sample should be 
taken while the oil is warm. 

In Tab. 1 a classification 
of values of the various in- 
vestigated parameters of in- 
sulating oil is given. An oil 
can be classified as good, 
acceptable or bad according 
to Tab . 1 [10] . 



Tab. 1: Classification of insulating oils 



Oil Parameters 


Good 


Acceptable 


Poor 


Color 


< 2 


- 


> 2 


Breakdown Voltage (kV) 


> 40 


30 - 40 


< 30 


Humidity [ppm] 


< 10 


10 - 25 


> 25 


Interfacial Tension (mN/m) 


> 28 


22 - 28 


< 22 


tan5 


< 0.1 


0.1 - 0.5 


> 0.5 


Resistivity (p) (GQ-m) 


> 3 


0.2 - 3.0 


< 0.2 



The sampling of oil from 
the fifty transformers 150/20 
kV was done with due care and 
according to the standard 
practice. In Figures 3-7, the 
results of the measurements 
are shown, regarding tan5, co- 
lor, humidity and interfacial 



tension respectively. Green 
color symbolizes the good sam- 
ples, yellow color the ac- 
ceptable samples whereas the 
red color shows the bad sam- 
ples. From Figs. 3 and 4, it 
is clear that the breakdown 
voltage and tan5 values of 
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most of the samples are very 
good. This is due to syste- 
matic control of the oil and 
the good maintenance. Figs. 3 
and 4 indicate that these two 
factors, which are related to 
ageing and oxidation, are 
relatively stable. 

In Fig. 5, humidity is in 
relatively acceptable levels. 
Only a small percentage of 
the transformers (4%) seems to 
have high humidity. Most of 
the samples are within the 
limits prescribed by the in- 
ternational standards. It 
must be emphasized that the 
humidity level is a parameter 
which changes easily, since 
the slightest contact of the 
oil with the atmosphere may 
change its characteristics. 
In Figs. 6 and 7 the results 
regarding the color as well 
as the interfacial tension 
are shown. 

As time passes by, the oxi- 
dation products change the 
oil color. Most of the sam- 
ples had a rather acceptable 
color. Even a color index of 
3 does not necessarily con- 
sist an objective indication 
of the oil quality. For this 
reason, color measurements 
should in fact be accompanied 
by other parameter measure- 
ments. In fact, although 36% 
of the samples showed a ra- 
ther dark color (Fig. 6) other 
parameter measurements indi- 
cated that these samples were 
good or acceptable. The pre- 



sence of foreign particles in 
combination with humidity may 
reduce the interfacial ten- 
sion of the oil. In the in- 
vestigated samples, a per- 
centage of 28% (Fig. 7) is 
characterized as poor. This, 
however, is not particularly 
annoying, if we take into ac- 
count for these samples also 
the other parameter measure- 
ments . 

The density of the inves- 
tigated oil samples was mea- 
sured in the generally ac- 
ceptable values, i.e. between 
0.85 and 0.92 gr/ml (with the 
lowest recorded being 0.85 
gr/ml, whereas the highest was 
0.91 gr/ml). Although the oil 
density does not consist per 
se an individual characteris- 
tic of the examined sample, 
its increase may imply an in- 
crease of degradation by- 
products. In the context of 
the present work the oil den- 
sity was used for the calcu- 
lation of the interfacial 
tension [11] , [12] . 

It can be said that, in 
general, the state of the in- 
vestigated transformer oil 
samples was more or less sa- 
tisfactory. In a few cases, 
there is a need of further 
filtering and possibly a se- 
cond sampling it should be 
carried out. Whereas no trans- 
former functioned with a par- 
ticularly bad oil, it is true 
that, with transformer age- 
ing, the oil suffers from so- 
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lid impurities, free and dis- 
solved water particles and 
dissolved air. Frequent sam- 
pling is necessary in order 
to ensure the good function- 
ing of such transformers. The 
results reported here are in 
line with those published be- 
fore [11] . Although a statis- 
tical approach of the whole 
subject is desirable (i.e. to 
try to correlate the data 
collected here as well as 
from previous published work 
with the pre-history of each 
individual transformer), it 
is difficult to be realized 
since the Greek Electricity 
System has transformers from 
a variety of manufacturers. 



This inhomogeneity of suppli- 
ers certainly renders the re- 
lation between the data col- 
lected with the pre-history 
of transformers very diffi- 
cult. It also should be noted 
that with the term "pre-his- 
tory" we mean the detailed 
registration of all faults, 
faulty conditions, lightning 
strokes, switching overvolt- 
ages etc. which have occurred 
in a transformer. Previous 
work done in this direction 
was only partially successful 
since at that time, pre-his- 
tory of transformers was only 
related to one parameter, 
namely that of breakdown 
strength [13], [14], [15]. 



2% 




Fig. 3: Graph of breakdown voltage results 
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Fig. 5: Graph of humidity results 
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It can be said that this 
piece of work here does not 
consist per se an original 
piece of research. This paper 
does not claim to have ex- 
plored new inroads regarding 
the mechanisms of breakdown 
of dielectric liquids, as for 
example in [16], or novel di- 
rections regarding new diag- 
nostic methods, as for exam- 
ple in [17]. The whole pur- 
pose of this paper is, how- 
ever, to show that monitoring 
work is necessary - not to 
say essential - for the cor- 
rect maintenance of trans- 
former oil in bigger trans- 
formers. An attempt to also 
relate the present state of 
transformer oil in the above 
mentioned transformers with 
their previous history will 
follow. 

Some Further Remarks 

One may ask which from the 
parameters investigated, are 
the most important. It is a 
difficult question to answer. 
However, if one has to choose 
between the aforementioned pa- 
rameters, he would most pro- 
bably select two of them, na- 
mely humidity and breakdown 
strength. If one looks care- 
fully Fig. 3 (referring to 
breakdown voltage results) 
and Fig. 5 (regarding humid- 
ity results), one may see 
that 98% of the investigated 
transformers had good oil 
breakdown strength and 96% of 



the investigated transformers 
had good or acceptable levels 
of humidity. This means that 
Fig. 3 and Fig. 5 are in more 
or less good agreement. Hu- 
midity plays a critical role, 
since it can contribute to a 
dramatic lowering of break- 
down strength, as was also in- 
dicated in some older but 
nevertheless classical publi- 
cations [18], [19], [20]. Ge- 
nerally speaking, inclusions 
of humidity more than about 

10 ppm (at normal tempera- 
ture) cause a lowering of 
breakdown strength. A low oil 
breakdown strength may imply 
that there are foreign parti- 
cles and/or admixtures in the 
oil. On the other hand, a 
high breakdown strength does 
not necessarily mean that the 

011 is good. It may be possi- 
ble that the quantity of fo- 
reign particles may not be 
sufficiently large, so that 
it can influence the break- 
down voltage [21] . 

The interfacial tension, 
although it gives an idea of 
the concentration of oxida- 
tion by-products in a trans- 
former oil, is not necessa- 
rily an indicator for defi- 
nite conclusions about the 
oil under investigation. This 
is because in warmer periods 
of the year, the oil tempera- 
ture increases and humidity 
may affect the oil more than 
in cooler periods [22] . 

Change of the oil color 
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may mean the existence of by- 
products or the presence of 
foreign particles. Although 
such a change may imply a 
certain degree of pollution 
of the oil under investiga- 
tion, the color by itself may 
not be considered as a very 
reliable indicator of the oil 
quality, as was shown in this 
paper and was also reported 
in [23]. 

Tan5 results (Fig. 4) 
match very well with Fig. 3 
results of breakdown 
strength. Tan5 changes as the 
oil degrades. Although tan5 
measurements cannot be taken 
as a sole criterion of oil 
quality, its results match ex- 
tremely well with the break- 
down strength data, in the 
context of the present paper. 
The low value of tan5 depends 
on the nature of the oil as 
well as on its processing 
[22], [24], 

A last remark should be 
made concerning the monito- 
ring of both transformers of 
150/20 kV and of distribution 
transformers (both kinds of 



transformers in the major 
Athens area) : it seems that 
the former have a larger per- 
centage of good/acceptable oil 
than the latter [25] . To va- 
lidate, however, this point, 
more work has to be done with 
transformers of both kinds 
from the major Athens area. 

Conclusion 

Insulating oil samples 
from transformers of 150/20 
kV have been investigated. In 
the context of this work, se- 
veral parameters - through 
the appropriate diagnostic 
techniques - affecting the 
state and lifetime of trans- 
former oil have been studied. 
No single parameter can fully 
describe the state of the oil 
of a transformer. The variety 
of parameters investigated he- 
re may give a more complete 
picture. In the context of 
this work, the majority of 
the investigated transformer 
oil samples were found to be 
good or acceptable. This 
points out to the further con- 
tinuing sampling at regular 
intervals . 
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Abstract 



Partial discharges (PD) may cause damage to high-voltage 
(h.-v.) electrical insulation and eventually breakdown. It 
is known, however, that sudden breakdowns occurred in indus- 
trial insulations after only a few months in service, al- 
though they had passed the suggested international specifi- 
cations tests. In this paper, we investigate the possibility 
of damage of an h.-v. insulation even below the inception 
voltage, giving the historical background which led to cer- 
tain thoughts. References to previous work are given and a 
differential equation is proposed regarding the possibility 
of having charging phenomena below the inception voltage. 



Keywords 
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Introduction 

PD may cause significant 
damage to the h.-v. insulation 
[1], [2]. PD may start from 

asperities on the electrodes, 
from enclosed voids, from 
fissures and/or from enclosed 
foreign particles [3] . Nume- 
rous publications tackled the 



voltage insulation, cavities, 
rents 

question of the relation be- 
tween PD parameters and insu- 
lation damage [4]-[7]. 

Although a voluminous 
amount of research has been 
performed on the questions of 
insulation damage and PD, no 
particular attention was gi- 
ven to possible charging 
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events below inception volt- 
age. However, such events we- 
re held responsible (at least 
partly) for the failure of 
insulation equipment back in 
the fifties and the sixties 
of the last century [8] -[11]. 
Such incidents were reported, 
regarding big electrical ma- 
chines as well as high volt- 
age switches that failed ab- 
ruptly, although they passed 
all tests the international 
technical committees were sug- 
gesting [8] . Which were these 
incidents? In the 1950's one 
major US manufacturer produ- 
ced and sold large rotating 
machines with cast polymer 
insulation. At that time all 
available diagnostic techni- 
ques suggested that there 
should be no problem at all 
from low leakage currents. 
Surprisingly enough, these ro- 
tating machines failed while 
in service in a few months. 
Another occasion of impres- 
sive failure was registered 
when another US manufacturer 
producing 15 kV class vacuum 
switches, reported that, al- 
though all normal tests at 
that time were done, unex- 
plained erosion of cast epoxy 
insulation occurred in these 
switches, only after a few 
years in service. Such were 
the incidents that led to the 
exploration of the below in- 
ception regime of charging 
phenomena [12] . Scientific 
evidence of "something" hap- 
pening below inception was 



provided by Filippini [13], 
who suggested that treeing 
propagation followed growth 
along current paths. Further- 
more, Brancato [14] indicated 
the possibility of relating 
current below inception with 
temperature rise. In his own 
words, "perhaps the most ex- 
citing but speculative pro- 
ject is the development of 
techniques to determine tem- 
perature rise of void sur- 
faces due to partial dis- 
charge. ... electrical aging 
may in reality be a chemical 
reaction in voids induced by 
the deposition of energy. 
This deposition is expended in 
the raising of the tempera- 
ture of the walls of the 
void, in creation of ozone and 
possibly nascent oxygen which 
tend to oxidize the surface 
chemical reactions which in 
the presence of water compo- 
nents form nitric or nitrous 
acids which deteriorate the 
void surfaces and finally me- 
chanical erosion by the re- 
sulting ion stream impinging 
on the void". In the same 
publication, Brancato also 
pointed out that measuring or 
sensing techniques may not be 
adequate for all classes of 
dielectric materials or to 
all aging environments and/or 
the inadequate sensitivity of 
some instruments to detect 
changes. Recent papers point 
out to the slow decline of 
inception voltage with time 
[15]. Is this also an indica- 
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tion of charging phenomena be- 
low the so-called inception? 
Would such phenomena play a 
role regarding the decrea- 
se/decline of inception volt- 
age? 

Hints about possible da- 
maging events below inception 
voltage were discussed in a 
brief report by E. Brancato 
back in 1991 [16] . Previous 
work on air gaps with a non- 
uniform electrode arrangement 
suggested that such charging 
events were possible [17]- 
[20] . To the best of our 
knowledge, until now there 
are indications of possible 
charging effects below incep- 
tion voltage but a definite 
proof is missing. 

It is the purpose of the 
present paper to offer the 
historical background of the 
approach regarding charging 
events below the inception 
voltage and to explore a bit 
further the events below in- 
ception. An equation, which 
was proposed many years ago 
by Bruning, is discussed. 

On the Problem of Inception 

Voltage - Historical Back- 
ground 

Inception voltage is cal- 
led the lowest voltage under 
a.c. conditions, at which re- 
petitive PD of a specified 
magnitude in successive cy- 
cles are observed, as the ap- 
plied voltage increases [21] . 
Above the inception level, an 



insulation may deteriorate de- 
pending on the magnitude of 
the imposed voltage. Things 
become more blurred when an 
insulation is subjected to an 
applied voltage which is lo- 
wer than the extinction volt- 
age (according to the seminal 
publication by Kelen [21], the 
extinction voltage is the 
voltage at which discharges 
of a specified magnitude will 
recur when an alternating volt- 
age, which exceeds the dis- 
charge-inception voltage, is 
reduced) . 

In his important paper, 
Bruning and colleagues [9], 
pointed out that unexpected 
insulation failures occurred 
despite the fact that the in- 
sulation was operating below 
the inception voltage. In 
this doctoral thesis, Bruning 
[8] pointed out that there 
were problems in tackling and 
understanding the basics of 
current -voltage relationships 
for new insulating materials 
and/or new devices. For exam- 
ple, the design failure of 
epoxy insulators, which used 
the application of criteria 
for porcelain discharge on- 
set, did not sufficiently cha- 
racterized low current per- 
formance with contamination. 
On another occasion, there was 
a lack of fundamental cur- 
rent vs. voltage analytical 
model for the design of a new 
type of transformer windings. 
In many cases of transformer 
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failures, discharge testing 
did not indicate low level 
pre-discharge activities with 
operation. Problems also aro- 
se with underground cables, 
where aging effects, exacer- 
bated by the presence of wa- 
ter, may give higher dis- 
charge currents which are be- 
low the level which can be 
satisfactorily measured. Bru- 
ning speculated that low le- 
vel discharge currents cause 
a destructive local tempera- 
ture effect. Problems related 
to the above questions were 
revealed in the USA with elec- 
trical machines, cables and 
high-voltage switches, which 
failed unexpectedly although 
all of them passed the re- 
quired tests prescribed by 
the international technical 
committees . 

All the above point to the 
fact that there is a need for 
a fresh approach for this pro- 
blem, an approach which may 
give a possible solution to 
charging effects below incep- 
tion. Efforts have been un- 
dertaken in Bruning's seminal 
publications [8], [9], where 
questions such as the follow- 
ing were put: 

a) what is the insulation 
leakage current-voltage cha- 
racteristic? 

b) what is the effect of 
such a leakage current? 

Efforts were also under- 
taken in some publications. 



albeit with a different insu- 
lating material and a much 
simpler electrode arrangement 
[17] -[20]. In those papers, 
indications were presented 
that effects below the so-cal- 
led inception level may exist. 
Moreover, experiments perfor- 
med in polyethylene samples 
with enclosed cavities, indi- 
cated that at relatively low 
voltages sharp current wave- 
forms were detected. This in- 
dicates sudden streamer PD 
mechanism in enclosed cavi- 
ties. In previous years, com- 
ments regarding the relation 
between such sudden PD wave- 
forms and events at or below 
the inception level were of- 
fered. It was speculated that 
there may be sudden bursts 
related to local rising of 
temperature in a cavity [22] - 
[26], 

Possible Relation Between 

PD Events at Inception Volt- 
age and Charging Events 

below Inception 

Bruning and colleagues [8], 
[9], [27] proposed an equa- 

tion describing the quasi- 
steady state for current flow 

V2v + W*Vs/s = 0 (1) 

where, V is the gradient and 
V 2 is the Laplace operator 
expressed as a function of 
the location, V is the local 
voltage and s the local con- 
ductivity of the fluid (in 
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our case air). They also pro- 
posed that it would be plau- 
sible to state that the local 
air temperature is related to 
the local power dissipation 
which is leakage current ti- 
mes the local voltage gra- 
dient. This being true, the 
local electrical conductivity 
is related to the local tem- 
perature. The question was as 
to whether such an equation 
was leading to a burst or 
sudden PD pulse even below 
inception. Would this be pos- 
sible? 

In fact the above equation 
was characterized as the 
"thermal model" in Bruning's 
Ph. D. Thesis [8]. It sug- 
gests a thermal model of 
gaseous conduction which un- 
der certain conditions of li- 
mited diffusion rate of ion- 
ized species, low radiant en- 
ergy loss, and low thermal 
conductance, a diffuse cur- 
rent flow may generate a dif- 
fuse conduction process. Si- 
milar pulseless regimes were 
suggested in [28]. Further- 
more, the same direction of 
research, without however be- 
ing so specific, was pointed 
out in [29], [30], where pseu- 
doglow regimes were observed. 
It may well be that the cur- 
rent not indicated on a con- 
ventional PD detector is a 
current below the detection 
level . 

More in detail, if one as- 
sumes that in the gas there 



is a dissipative current flow, 
using classical notation. 



w* 


J = 0 


(2) 


J = 


sE 


(3) 


V*sE = 0, but E = W 


(4) 


and 


consequently. 




V2V 


+ W*Vs/s = 0 


(5) 



The above equation, in the 
case of a constant conducti- 
vity s, reduces to the well 
known Laplace/Poisson equa- 
tion. However, for a gaseous 
conduction process, where the 
there is a partial variable 
ionization as J varies with 
time and position, we have 
finally Eq . (5) . 

As was noted in a previous 
paper [27], solution of the 
distributed - perhaps pulse- 
less current flow - then pro- 
ceeds from ionization deter- 
mined from both Saha's equa- 
tion and the heat balance de- 
termining the local tempera- 
ture. The complexity of the 
transport processes have not 
permitted ab initio calcula- 
tions. However, approxima- 
tions indicate a variation of 
current density and field 
strengths, which indicate the 
possibility of current below 
the detection level. Whether 
the continuous thermal condu- 
ction model arises from true 
continuous conduction or from 
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pulses too short to be de- 
tectable by conventional PD 
detectors, it is a question 
in need of an answer. It is 
interesting to note, however, 
that sudden bursts of pulses 
at or below the so-called in- 
ception voltage were observed 
recently, not only in conven- 
tional polymers but also - 
albeit less frequently - in 
nanocomposites [31], [32]. 

Research by Bruning and 
colleagues [9], [10] indicated 
that damage observed below 
inception was very similar to 
the damage above inception. 
In these papers, indications 
were offered that chemical 
changes below inception were 
similar to those noted above 
inception. Such indications 
gave strong ground for sug- 
gesting that "something" goes 
on below the inception volt- 
age. Chemical changes that 
indicate chemical deteriora- 
tion below inception imply 
the existence of charging phe- 
nomena that may render shor- 
ter the lifetime of the insu- 
lation. Charging effects be- 
low the inception voltage may 
imply that there is no volt- 
age below which no deteriora- 
tion takes place. More pre- 
cisely, it is known that the 
lifetime L (time to failure) 
is given by the formula 

L = c(V - Vo) -k (6) 

where, V is the applied volt- 



age, Vo is the voltage below 
which no deterioration takes 
place and k is a constant, k, 
Vo, and c are constants de- 
pending on the insulation ma- 
terial (it should be noted 
that original form of Eq . (6), 
referred to in [33], was 

L = (A - a)k(v - aA)-k 

where, L is life in minutes, 
A is the 1 min strength, a is 
a material constant, k is an- 
other material constant and V 
is the applied voltage). How- 
ever, if one assumes that the- 
re is damage below the incep- 
tion voltage, then 

Vo = 0 (7) 

and Eq . (6) becomes 

L = cV-k (8) 

As was noted in [9], "most 
equipment designers use this 
empirical relation without 
having settled the fundamen- 
tal question as to whether Vo 
= 0 or not, since empirical 
experiments in reasonable ti- 
me periods cannot distinguish 
between the two forms", i.e. 
that of Eq . (6) from that of 
Eq . (8). Ambiguities regard- 
ing Eq . (6) still remain to 

this day, especially if one 
looks at insulation lifetime 
models [34] . The question of 
whether V0 = 0 is also rele- 
vant to the sensitivity of PD 
detector as well as to the 
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thickness of the examined in- 
sulation, deserves to be fur- 
ther explored. It was re- 
ported before, that scaling 
laws exist, i.e. a thicker 
insulation requires a more 
sensitive PD detector [35], 
[36] . Needless to say that 
the sensitivity of a PD de- 
tector determines the incep- 
tion voltage of either an in- 
sulating material or a full 
insulation system, or in the 
words of one of the most dis- 
tinguished researchers of the 
good old generation, E. Bran- 
cato, "... electric stresses 
in the absence of internal 
discharges can cause changes 
in material properties. Some 
ascribe the changes to elec- 
trochemical reactions, while 
others express the suspicion 
that the observed changes are 
really caused by partial dis- 
charges but the corona de- 
tecting system is too insen- 
sitive" [37] . 

From the above it is indi- 
cated that, although there is 
not yet a universal name to 
it, "something" must occur 
below the inception voltage. 
Whether this "something" can 
be referred to as "charging 
phenomenon" or can be mani- 
fested as "signal", is not yet 
clear. Moreover, if one sees 
the more practical questions 
and if one tackles the prob- 
lems of the role of antioxi- 
dants in cable insulation. 



one may say that antioxidants 
are incorporated to prevent 
premature degradation during 
extrusion, but the role of 
the residual antioxidant on 
insulation response to aging 
stresses has not yet been 
quantified, especially in 
voltages below the inception 
level [38] . 

Conclusion 

The present paper is an 
introduction to the problem 
of possible charging effects 
below the so-called inception 
voltage. The historical back- 
ground of an approach - dif- 
ferent from the usual ap- 
proaches - is given as well 
as the indications that char- 
ging phenomena below incep- 
tion may exist. Below incep- 
tion sudden pulses were ob- 
served with an electrode ar- 
rangement of point-plane with 
air as insulating material. 
It is speculated that minute 
abnormalities on the cavity 
surface may act as emission 
sites and thus provoke small 
charging phenomena not easily 
detected by normal conven- 
tional PD detectors. An equa- 
tion regarding the current 
flow in a cavity was given 
and commented upon. In a fu- 
ture publication, a solution 
of the said differential equa- 
tion will be given together 
with appropriate comments. 
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Abstract 

A planar Delta-Cross shaped loop antenna is proposed. The 
presentation includes results from an analytic study of the 
radiation pattern, as well as, from a simulation study for 
both the electric and radiation characteristics: input imped- 
ance, standing wave ratio, radiation pattern and directivity. 
The loop was initially shaped as 4 non-overlapping equilat- 
eral triangles on a plane, symmetrically oriented around a 
common vertex, center fed at one triangle base. In order to 
improve the antenna characteristics, its shape was then modi- 
fied by equally changing the triangle base angle-while keep- 
ing the loop length constant and equal to 2 wavelengths WL. 
In this way, the final antenna loop was shaped with 4 isosce- 
les triangles. The analytical and simulated results for the 
radiation pattern were found to be in good agreement. Fur- 
thermore, a comparison with antenna's dipole counterpart 
characteristics showed a much better performance of the pro- 
posed antenna. 



Keywords 

Cross loop antenna, delta elements, analysis, simulation, 
improvement 



Introduction 

The cross loop antenna con- 
sisting of four delta ele- 
ments was examined as a pro- 
totype antenna during prepa- 
ration of an EECE diploma 
thesis [1] . The main avail- 
able tools for its study were 
the antenna theory [2], the 
[RadPat4W] computer program 



for antenna patterns [3] the 
[RICHWIRE] simulation program 
[4] and the mini-Suite of 
software tools [5] . Its sim- 
ple plane figure of double 
symmetry with respect to two 
axes which means easy geomet- 
rical representation for the 
theoretical consideration of 
its radiation pattern as well 
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as for the simulation and 
easy construction, was the 
basic reason for this pro- 
posal. The initial antenna 
consists of four equilateral 
triangles with side length 
A/6 and almost coincident 
their four vertices. Thus, 
the total length of the an- 
tenna is 2A. The feed source 
was at the middle of the base 
of one of the triangles as it 
is shown in Fig. 1. 




48° angle was selected as the 
improved Delta-Cross Shaped 
Loop antenna. 

Analysis 

The initial antenna was 
analytical studied consider- 
ing standing waves, that is, 
sinusoidal, current distribu- 
tion in a parallel wire 
transmission line. Since the 
perimeter of each equilateral 
triangle is A/2, the total 
antenna length is equal to 2A 
and it corresponds to a prop- 
erly formed short circuited 
transmission line. However, 
the study was carried out in 
a 5A/2 piece of open cir- 
cuited transmission line, 
formed properly as in Fig. 2, 
since the available formula 
of equation (1) concerns open 
circuited two parallel wired 
transmission lines 

i(£) = i sin((3(h - |£|)) (1) 



Fig. 1: Planar Delta-Cross 
Shaped Loop antenna 

Then, the antenna was mod- 
ified only with respect to 
the base angles of the trian- 
gles, keeping constant the 
triangle perimeter equal to 
A/2, resulting a cross an- 
tenna of four isosceles tri- 
angles. The angle was varied 
between 48° and 78°, in order 
to improve both of its elec- 
trical and radiative charac- 
teristics and the one with 



with 



5 

-h<f<+h where h = — A . 

4 

The above current distribu- 
tion (1) was used to evaluate 
the radiation pattern of the 
antenna through the relations 



iP R k r . 

E = e r PF 







PF = 



J ±(^’)e ;L|3 ' r, d t 

<A 



( 2 ) 

( 3 ) 
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According to the basic 
standing wave theory: i) the 
direction of current changes 
per A/2 segment, ii) between 
two successive standing wave 
nodes the current phase re- 
mains the same and iii) if a 
source is between two succes- 
sive standing waves the cur- 
rent direction does not chan- 
ges. These principles were 
applied at the properly ben- 
ded 2-wire transmission line 
of Fig 2, keeping in mind 
that the Kirchhoff laws must 
be simultaneously satisfied. 




Fig. 2: 5A/2 properly bended 
2-wire transmission line 



The double arrows indicate 
the direction of it product 

for every one segment, and 
the dotted circles at the al- 
most coincident vertices of 
the four equilateral trian- 
gles correspond to a zero in 
the current distribution. 

The two parallel A/4 seg- 
ments are almost coincident 
with opposite currents, so 
they are mutually canceled as 
for the radiation and they 
are excluded from further stu- 
dy. Thus, the remaining seg- 
ments form the desired an- 
tenna, with total length 2A 
and a current maximum at the 
source, as it is shown in 
Fig. 3. 

For every one segment 0- 
12, of Fig. 3, we had to de- 
termine four quantities: the 

length of its starting point 
f A , the length of its end 
point t T , the position vector 
of its center R k , as if they 
were part of the correspond- 
ing line of 2A length, and 
its unit direction vector t , 

I 

with its direction to be in 
all segments from t k to t T . 
The planar antenna was ar- 
ranged on yOz as it is shown 
in Fig. 4, along with all the 
R k position vectors as they 
were determined by the rela- 
tion 

R k =0A v -tf A f v (4) 

V V I 
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and Tab. 1 contains all the 
above mentioned geometrical 
quantities . 



i m 




i 

Fig. 3: 13 segments 



Then through (l)-(3), that 
is, after the evaluation of 
the definite integral for the 
Pattern Factor PF [2], [6] 

the radiation pattern of all 
the 13 segments was formu- 
lated and according to the 
law of superposition the to- 
tal radiation pattern is 
given by 



: ve = e =t e 0 + £ (p $ = 




which, after a lot of sym- 
bolic manipulation where the 
real parts are mutually elim- 
inated, leads to a more sim- 
plified expression 




and by the symmetry of the 
geometry radiation patterns 
of 0th, 12th and 6th wire are 
combined, as well as the five 
(5) couples: 1, 7-2, 8-3, 
9-4, 10 and 5, 11 as it is 
obvious from Tab. 1. Thus, 
the total radiation pattern is 



E Eg 12 6 + E-l 7 + E 2 g + 



+ ^3,9 + E 4 ^ 10 + E 511 



Simulation 

The simulation process was 
carried out in terms of [RICH- 
WIRE] program. The center fre- 
quency was 1111 [MHz] where 
the wavelength is A = 0.27 [m] 
and the wire radius was 1 
[mm] . We studied the antenna 
characteristics both electric 
and electromagnetic (SWR, In- 
put Impedance, Directivity, 
Radiation Pattern) in the 
range [600, 1300] MHz around 
the center frequency with a 
10 [MHz] step. 
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Fig. 4: 13 position vectors 
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Tab. 1: Geometrical Data 



# 






R k 


P R k r 








y 

i 


z 

1 




0 


-A 


11 A 

12 A 


-A 


H a 
12 A 


V3 TI 

- 2 nsin 0 sincp+— — cosG 


1 


-H a 


-H a 




- 3,g x 


3^ sinGsinco 3 V3 ti q 


12 


12 


8 A 


. O _L II U J _L 1 1 VI/ m V/U Jw 

4 4 


2 


-H a 


-H a 


- 33 a 


& 


3 A 3 ^ °i nfi c i nrn+ p n ° G 


12 A 


12 A 


8 A 


4 4 


3 






12 


A 

2 


-HZsin 0 sincp-ncos 0 


4 






H a 
~ 8 _A 


A 

8 


H— sin 0 sin(p+Hos 0 
4 4 


5 




A 

12 


A 

8 


H 


- ^sin 0 sincp-HZ COS 0 


6 


A 

12 


A 

12 


0 


H a 
“ 12 A 


-""cose 

6 


7 


A 

12 


& 


A 

8 


H a 
8 ~ A 


^sin 0 sincp-HZ COS 0 
4 4 


8 


& 


& 




A 

8 


- H^sin 0 sincp+Ilcos 0 


9 


& 


h* 




A 

2 


sin 0 sincp ttcos 0 


10 




H a 


3 3 a 


£> 


3 3 ^ 0 - p|Qo-j nrn-|- 3^ p n ° fl 


12 A 


12 A 


8 A 


^ o _LI 1 U o _LI 1 v|^ T ^ OUoO 


11 


H a 


H a 


4* 


- 3 H a 


3 n o -j i nrn 3 Vs" TT ~ 


12 A 


12 


8 A 


_ o _L 1 1 U o _L 1 1 Vl/ U U o U 

4 4 


12 


11 A 

12 A 


A 


A 


H a 
12 A 


_ ... \'3tt . 

2 nsin 0 sincp+— — cos 0 
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# 


t 

1 




^0 

1 






y 

i 


Z 

1 








0 


-i 


0 


- sinBsincp 


- cosBsincp 


-coscp 


1 


1 

2 


Vs 

2 


^-sinBsincp- cosB 


^-cosBsincp+V^-sinB 


1 

— coscp 


2 


Vs 

2 


1 

2 


- -^-sinBsincp+^-cosB 


- -^-cosBsincp- ^-sinB 


Vs 

-—coscp 


3 


0 


-1 


- cosB 


sinB 


0 


4 


V3 

2 


1 

2 


V^-sinBsincp+^-cosB 


V^-cosBsincp-^-sinB 


Vs 

—coscp 


5 


1 

~2 


Vs 

2 


- ^ sinBsincp - V^cosB 


- ^-cosBsincp+V^-sinB 


1 

--coscp 


6 


1 


0 


sinBsincp 


cosBsincp 


coscp 


7 


1 

~2 


Vs 

2 


- ^-sinBsincp+V^-cosB 


- ^-cosBsincp sinB 


1 

--coscp 


8 


Vs 

2 


1 

“2 


^-sinBsincp - ^ cosB 


V^-cosBsincp+^-sinB 


Vs 

—coscp 


9 


0 


1 


cosB 


- sinB 


0 


10 


Vs 

2 


1 

“ 2 


- ^-sinBsincp- ^ cosB 


- 4^-cosBsincp+^-sinB 


Vs 

-—coscp 


11 


1 

2 


Vs 

2 


^-sinBsincp+^-cosB 


^-cosBsincp- sinB 


1 

— coscp 


12 


-1 


0 


- sinBsincp 


- cosBsincp 


-coscp 
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The 


total 


number of 


96 


segments for 


the whole 


an- 


tenna. 


where 


all sides 


had 


equal number 


of points. 


was 



selected for the simulation, 
after the investigation for 
13, 24, 48, 72, 96, 120, 144, 
168, 192, 216, 240 segments. 
The above mentioned charac- 
teristics remained almost 
constant after this number of 
96 segments. An appropriate 
computer program was devel- 
oped in Fortran to create the 
input file with antenna geom- 
etry for [RICHWIRE] . 

A very good agreement be- 
tween simulated and analyti- 



cal produced radiation pat- 
terns is shown in Fig. 5 for 
the three main planes in 1111 
[MHz], although there is a 
point for discussion here. A 
closer examination, after 7 
years from the initial work 
[1], of what was really com- 
pared in Fig. 5, proved that 
the analytical patterns illus- 
trated corresponds to rather 
an upper bound of the abso- 
lute radiation pattern since 
it results from the sum of 
the 13 absolute radiation 
patterns and not from the 
total absolute radiation pat- 
tern . 



— Analytical Richwire 




Fig. 5: Radiation patterns in 1111 [MHz] 



Furthermore, the simula- 
tion model, as it is shown in 
Fig. 6, consists inevitable 
not in four equal delta ele- 
ments but in four truncated 
ones, as care must be taken 

a) to avoid the overlapping 
of the wires in a future con- 
struction and 

b) to reduce the strong cou- 



pling caused by the proximity 
effect . 

Thus, the minimum permit- 
ted square region with side 
(2a+a/10), where "a" is the 
used wire radius, was removed 
from the center intersection 
point of the four delta ele- 
ments, while in the same 
time, in order to achieve the 
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connection of the neighboring 
sides, the base angles of ev- 
ery element are slightly 
greater that 60°, and the to- 
tal length of such formed el- 
ement was slightly greater 
than the initial A/2. 

Z Z 

i i 





Fig. 6: Simulation models 

Extremely high values of 
Standing Wave Ratio [SWR] in 
1111 [MHz] for the three 
characteristic impedances 50, 
75 and 300 [Q] as it is shown 



in Fig. 7, made clear the 
need for further research on 
the antenna's geometry. 

The antenna was modified 
only with respect to the base 
angles for constant delta 
perimeter, resulting a Delta- 
cross Loop antenna of four 
isosceles triangles. A com- 
puter program was also devel- 
oped in Fortran producing the 
corresponding input data for 
[RICHWIRE] . 

The investigation was per- 
formed from the value of 48° 
for the base angles to 78° in 
step of 1° at the frequency 
of 1111 [MHz]. The best per- 
formance was noticed between 
48° and 52° and the model, 
shown in Fig. 8, with 48° 
base angles was selected. 
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Fig. 8: Improved antenna 

In Figs. 9-13 SWR, Input 
Impedance [Zinp] as Input Re- 
sistance Rinp and Input Reac- 
tance Xinp in [Q] and Direc- 
tivity D in [dB] of the im- 
proved antenna together with 
that of the initial one and 
of 2A dipole, are shown re- 
spectively. Fig. 12 shows a 
properly selected window of 
the total frequency range for 
the Input Impedance in order 
to clarify the improvement of 
the modified antenna. In 
Figs. 14 - 15 a comparison of 
the radiation patterns of the 
three antennas is given at 
five frequencies in the three 
main planes zOx, xOy and yOz . 
In Fig. 16 the 3D analytical 
radiation pattern at 1111 
[MHz] and in Fig. 17 the cor- 
responding 3D patterns in all 
frequencies are shown. 

Conclusion 

Although, the agreement of 
results between simulation 
and analysis was very good, 



the characteristics of the 
initial antenna were not sat- 
isfactory. An improvement 
process was needed and the 
research was here limited to 
the change of the base angles 
only to the rather small 
range of -12°, +18° around 

the initial 60° base angle. 
Thus, the improved antenna 
had the best performance 
about 100 MHz lower of the 
desired frequency of 1111 
[MHz], as it is shown in 
Figs. 9-13. The SWR for the 
50 [Q] characteristic imped- 
ance has its best value 
closer to 1111 [MHz] while 
the initial antenna and the 
2A dipole of equivalent 
length had a relative small 
value of SWR at about 800 
[MHz] (Fig. 9). It is also 
obvious that the Input Imped- 
ance as well as Directivity 
is relatively stable and 
slightly better for the im- 
proved antenna and much bet- 
ter from the initial one 
(Figs. 10 - 13). Correspond- 
ing values for SWR, Input Im- 
pedance and Directivity for 
the three antennas are given 
in Tab. 2. 

The radiation pattern of 
the improved antenna is clo- 
ser to that of a centered fed 
sinusoidal A/2 dipole, avoid- 
ing the dispersion of radia- 
tion in multiple directions 
of the rather long 2A dipole 
antenna (Figs. 14 - 15, 17). 
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Tab. 1: SWR, D, Zinp 





SWR 


D [dB] 


Zinp [Q] 


1111 [MHz] 


60 0 


43.7 


2.49 


40 - i 287 


2A 


26.1 


4.09 


577 - i 647 


o 

00 


6.55 


3.10 


70 + i 127 


800 [MHz] 


60° 


1.44 


2.61 


36 + i 5 


2A 


2.55 


3.21 


88 - i 52 


o 

00 


42.7 


2.30 


24 - i 219 


1010 [MHz] 


60° 


29.6 


2.85 


1297 + i 485 


2A 


20.4 


3.74 


993 + i 168 


o 

00 


1.26 


2.90 


40 + i 3 



An extensive investigation 
concerning the analytical de- 
termination of the antenna 
radiation pattern and its 
normalized absolute values is 
in the immediate future plans 
of the authors in order to 
achieve the clarification of 
its relation with the one 
produced by simulation. 

A wider range of the base 
angle variation or the varia- 
tion of the total length of 
each delta element are some 
design ideas for a future re- 
search work. Nevertheless, 
the final step is, without 
doubt, the construction and 
measurement of the proposed 
antenna together with the 
comparison of the correspond- 
ing results. 
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Frequency [MHz] 



. 10: Rinp for improved-48 0 , 2A and initial-60 0 antenna 




Frequency [MHz] 

Fig. 11: Xinp for improved-48 0 , 2A and initial-60 0 antenna 
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800 1000 1200 800 1000 1200 
Frequency [MHz] Frequency [MHz] 

Fig. 12: Rinp, Xinp for the three antennas in [800, 1200] MHz 




600 700 800 900 1000 1100 1200 1300 

Frequency [MHz] 



Fig. 13: Directivity of the three antennas 
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— 48° — 2X 60“ 






Fig. 15: Radiation patterns at 1111 and 1300 [MHz] 



z 

60 ° 1 1 

1111 [MHZ ] 




Fig. 16: 3D Analytical Radiation Pattern 
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1010 [MHz] 



llll[MHz] 



1300 [MHz] 



Fig. 17: 3D Simulated Radiation Patterns 
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www. mfa . org/collect ion s/object /plate -153852 

The simple meandros is a symbol for eternal immortality. 

"Warrior with a phiale", ca. 480-460 BC, 

by Berliner Maler, 

Museo Archeologico Regionale "Antonio Salinas" di Palermo: 
" http : //commons .wikimedia. or g/wiki/File : War rior_MAR_Palermo_ 

NI2134.jpg" 
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